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Abstract 


This document details the methodology and use of the CAM- 
RAD. Mod 1/HIRES codes, which were developed at NASA Langley Research 
Center for the prediction of helicopter harmonic and Blade- Vortex Interac- 
tion (BVI) noise. CAMRAD.Modi is a substantially modified version of 
the performance/t rim /wake code C AM RAD. High resolution blade load- 
ing is determined in post-processing by HIRES and an associated indicia! 
aerodynamics code. Extensive capabilities of importance to noise prediction 
accuracy are documented, including a multi-core tip vortex roll-up wake 
model, higher harmonic and individual blade control, tunnel and fuselage 
correction input, diagnostic blade motion input, and interfaces for acoustic 
and CFD aerodynamic codes. Modifications and new code capabilities are 
documented with examples. A user's job preparation guide and listings of 
variables and namelists are given. 




Chapter 1 

Introduction 


1.1 Introduction 

With growing noise restrictions being imposed on rotorcraft, a means to 
accurately and efficiently predict noise generated by a wide variety of ro- 
torcraft configurations is needed. Many of the existing rotorcraft computer 
codes that are available are intended to calculate only rotorcraft perfor- 
mance quantities. The calculation requirements of the rotor system in a 
performance analysis often involves only the lowest frequency loading re- 
sults. For example, in level steady flight, a simple performance analysis 
might only require knowledge of the mean rotor thrust and drag. To pre- 
dict such quantities, a high resolution loading calculation on the rotor is 
not necessary. But, for noise calculations, a detailed, high resolution ra- 
dial and azimuthal loading solution is needed in order to accurately define 
events such as Blade- Vortex Interaction (BVI) noise. A computer code sys- 
tem designed to fill this requirement is presented in this documentation. 
This computer code system uses, as a “base" code, the original 1980 version 
of the Comprehensive Analytical Model for Rotorcraft Aerodynamics and 
Dynamics (CAMRAD) (Ref. [1] and [2]). This document is intended to sup- 
plement the original documentation of the CAMRAD code, not to replace 
it. It is assumed that the reader is already familiar with the original version 
of CAMRAD. Since this document enumerates the changes that have been 
made to the original CAMRAD code to create the code system that is now 
collectively known as “CAMRAD. Mod 1/HIRES” (Ref. [3]), it is intended 
to document the modifications and to be a reference for the new coding. 
Presented are the three major parts of the code: CAMRAD.Modl, HIRES, 
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and the Indicia] Post- Processor. It also contains an updated Users' Manual 
that lists ail variable inputs to the code system. 


1.2 Organization of Documentation 

This chapter provides an introduction, a system overview, an outline of 
the documentation, and a discussion of the sample cases to be used in this 
document. Chapter 2 deals with the modifications made to the low resolu- 
tion part of the original version of CAMRAD to obtain CAMRAD.Modl. 
Within each chapter, a section discusses each set of modifications. Where 
applicable, the first subsection of each section is an introduction to discuss 
the motivation for the modification and the methods used. Subsequent sub- 
sections discuss details of the modification including the actual code changes. 
When applicable, the last subsection of each section discusses how the mod- 
ification is related to the high resolution modifications. 

Chapter 3 discusses the changes made to include a high resolution post- 
processor, known as HIRES. This coding is part of the CAMRAD.Modl 
code, but is executed after the trim loop of the low resolution calculations. 

Chapter 4 deals with the code known as the Indicia] Post- Processor 
(IPP). This code is a standalone code that incorporates many aspects of the 
works of T.S. Reddoes and Cordon Leishmann with regard to empirical use 
of indicia! aerodynamic functions. Many of their formulations are directly 
applicable in the code, but others were modified such that they could be 
cast into a form compatible with the CAMRAD.Modl system. 

Chapter 5 is the Users’ Manual and describes namelist inputs, informa- 
tion on codes needed to prepare input data, and other general user-related 
information for CAMRAD.Modl. This chapter is intended to be a supple- 
ment to Reference [2]. 

1.3 System Overview 

A code system has been developed to expand the capabilities of previ- 
ous rotorcraft performance and noise codes. As a “base” code, the original 
version of CAMRAD was chosen. The original CAMRAD version is capa- 
ble of performing comprehensive rotorcraft calculations such as performance 
and low resolution loading calculations, for various rotorcraft configurations 
including a single rotor in a wind tunnel, a conventional helicopter, a tan- 
dem rotorcraft, a coaxial rotorcraft, and a tiltrotor. The original analysis 
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is divided into a several parts. First, a “Trim” analysis determines the ro- 
torcraft configuration orientation, control settings, etc.) required to 

match a specified flight condition. Second, a “Flutter" analysis linearizes 
the rotorcraft equations of motion about the trimmed configuration and 
determines eigenvalues, etc. Third, the “Transient” analysis determines a 
rotorcraft non-equilibrium response to a particular input such as a gust. The 
majority of the work in this document is related to the Trim analysis and 
post- processing of the Trim analysis results. Since the focus of the work in 
the CAMRAD.Mod 1 effort is on the Trim solution and post -processing of the 
Trim results, no effort has been made to update the Flutter and Transient 
analyses; as such, their usage in CAMRAD.Mod! is neither recommended 
nor are they discussed further in this document. 

Since most of the following document pertains to the Trim analysis and 
new follow-on procedures, a brief review of the trim process used in CAM- 
RAD is in order. It is assumed here that the reader is somewhat familiar 
with the CAMRAD prediction capability and the details of the rotor trim 
methods used t herein. In Figure 1.1, the box labeled CAMRAD.Mod 1, that 
includes the Trim, Transient, and Flutter analyses, represents the low res- 
olution portion of CAMRAD.Mod 1. This set of analyses, discussed above, 
are analogous to the original CAMRAD code. Also shown are several in- 
put paths and an additional output path into HIRKS. These paths will be 
discussed later. Figure 1.2 is an expansion of the Trim box of Figure 1.1. 
As the trim procedure is being carried out, all processes except the current 
process are held constant, as in the original CAMRAD code. 

For example, first, the wake and the wake influence coefficients are de- 
termined for a fixed configuration. That is, the influence coefficients are 
determined for fixed blade motion, fixed circulation, fixed blade controls, 
etc. Once these wake influence coefficients are known, they are held fixed 
and the next stage proceeds. After the wake influence coefficients are known, 
the necessary control settings are determined to match the target flight con- 
dition. A modified Newton-Raphson technique is used to increment the 
controls to determine a guess at the actual control settings required to meet 
the target flight condition. Then, with the wake influence coefficients fixed, 
the blade controls fixed, and blade motion fixed, the circulation distribu- 
tion is determined. With the just calculated circulation distribution, the 
blade motion is recalculated. This motion/circulation iteration continues 
until successive iteration differences for both have converged to below a 
prescribed tolerance. With the motion/circulation iterations converged, the 
next guess at a control setting is made, and the process is continued until all 
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Figure 1.1: Flowchart of CAMRAD.Modl and HIRES with associated in- 
puts. 
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iterations are converged and the target flight condition is met. This entire 
trim process may be successively repeated for several (or the same) wake 
models: a uniform inflow model, prescribed wake model, free wake model, 
and rollup wake model. 

Figure 1.1 also illustrates a flow chart of the input requirements for the 
CAMRAD.Modl code system. First, airfoil characteristics (he., lift, drag, 
and moment coefficients) are normally available in a standard “CN1 airfoil 
table”. These airfoil tables are used by the airfoil preparation program, here- 
after denoted as AIRFOIL, to generate a binary airfoil file, represented by 
"airfoil, tab” in the figure. Alternatively, airfoil characteristics may be gener- 
ated using namelist inputs to AIRFOIL. Second, BLOCKDATA information 
for the rotorcraft is prepared by the input preparation program, hereafter 
denoted as BLOCKFILE. This information is also normally converted to a 
binary input file, labeled "inputfile.bin” in the figure, for use in the analysis. 
In addition to the BLOCKDATA, namelist inputs are used to set specific 
run conditions such as RPM, advance ratio, etc. These namelists are located 
in the script file or command file used to run the analysis. Also, other files 
may be input to the analysis for use in the tunnel/fuselage correction model, 
denoted "TF corrs” in the figure. 

Several new output paths emerging from CAMRAD.Modl were intro- 
duced to predict such events as B\T and system noise. Figure 1.3 again 
shows the CAMRAD.Modl box (a box containing a trim loop, a transient 
loop, and a flutter loop), and the HIRES box with their associated outputs. 
As seen in this figure, several new branches have been made out from the 
end of the trim loop; none of these branches existed for the original version 
of CAMRAD. The first of these branches, the Trim-FPRBVI branch, along 
with the possible return branch FPRBVI-Trim, is called the CFD interface 
and is discussed later in this document. This branch uses the low resolution 
wake and blade position information and applies an external CFD code, in 
this case, FPRBVI (Ref. [4]), to calculate the high resolution loading uti- 
lized by the rotorcraft noise code, WOPWOP (Ref. [5]), to predict BVI 
noise. A return path from FPRBVI to the trim loop is possible in an open- 
loop manner, to take advantage of the loading calculated by the CFD code 
in the calculation of rotorcraft trim. Along the same branch, once the tone 
noise has been calculated by WOPWOP, the ROTONET (Ref. [6]) code 
system could be applied to compute propagated noise. 

In the second branch (Trim-ROTONET), the flight condition, the blade 
position, and the low resolution loading information is made available for use 
in the systems noise code ROTONET to predict tone noise and propagation 
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H ROTONET 


Figure 1.3: CAMRAD.Modl and HIRES output paths to other codes. 









effects. In the third branch (Trim- WOPWOP), the flight condition, the 
blade position, and low resolution loading information is made available 
for use in the rotor tone noise prediction code WOPWOP. The WOPWOP 
results could then be used in ROTONET to calculate propagated noise. 
Both of these branches are discussed later in the document. Note in the 
figure that both of these paths are surrounded by dashed lines. This is done 
to indicate that, even though these paths exist in the code, their use is not 
recommended since they do not include higher harmonic loads and their use 
may produce misleading results. 

The next branch involves the extensions of CAMRAD to include a high 
resolution wake and/or loading calculation known as HIRES. Figure 1.4 
provides a brief introduction to the solution procedures used in this portion 
of the code; a more detailed discussion is provided later in the document. 
After the trim solution has been obtained, the far wake influence coefficients 
are obtained using a high resolution reconstruction of the blade position and 
wake position. With the far wake influence coefficients known, airloads due 
to the far wake effects may be calculated. To account for the near wake 
effects, two choices are possible. One choice is a near wake lattice model 
used to calculate influence coefficients of a near wake lattice, followed by 
an airload analysis. Although this method exists in the code, is has not 
been exercised thoroughly or validated. The other choice is an Judicial 
Post- Processor (IPP) code that accounts for the near w r ake effects using 
indicia] aerodynamic functions. Once the airloads are known, they may be 
used in the rotorcraft tone noise code WOPWOP, optionally followed by the 
ROTONET code to account for noise propagation effects. 

1.4 Sample Case 

Throughout Chapter 2 of this documentation, a sample case will be used 
to illustrate modifications to CAMRAD. Modi. Examples of results from 
HIRES and the IPP can be found in the literature (Ref. [3]) and thus are 
not presented here. In general, when comparisons are being made between 
this sample case and the same case showing the modification, the upper plot 
of a given figure is the result from the sample case. The lower plot is the 
same sample case including the particular modification in question. The 
sample case is a model BO- 105 hingless rotor in a wind tunnel. Some of the 
properties used in the sample case are listed in Table 1.1. 



Far Wake Influence 
Coefficients 



Figure 1.4: The major computational loops within 


Far Wake 
Calculations 


Near Wake 
Calculations 


HIRES. 



Table 1.1: Sample Case Information 


radius 

2.0 meters 

chord 

0.121 meters (rectangular planform) 

number of blades 

4 

flap hinge location 

(none) 

lag hinge location 

(none) 

sweep of quarter-chord 

0.0 degrees 

airfoil section 

NACA 23012 

precone 

0.0 degrees 

nominal advance ratio 

0.15 

nominal RPM 

1041.0 

nominal shaft tilt 

5.3 degrees (aft tilt) 

nominal Ct/signia 

0.05007 



Chapter 2 

CAMRAD.Modl 


In this chapter, modifications made to the low resolution part of CAM- 
RAD are presented. There are a number of sections, each describing in 
detail the specific modification and/or enhancements made to CAMRAD. 
Where appropriate original CAMRAD predictions are compared to CAM- 
RAD.Modl (modified CAMRAD) predictions. 

2.1 Changes to Free Wake Azimuthal Resolution 

2.1.1 Introduction 

There are two vortex wake models in CAMRAD to determine vortex 
geometry. These models are the rigid (or prescribed) wake model and the 
free wake model. One of these wake models is used during the wake in- 
fluence coefficients calculation (shown in Figure 1.2) to determine the tip 
vortex geometry. As for the free wake geometry model, CAMRAD.Modl 
relies on the Scully Free Wake method (Ref. [7]) as does the original version 
of CAMRAD. However, for CAMRAD.Modl, a higher resolution free wake 
analysis is desired and the smallest possible azimuth step size was changed 
from 15 to 10 . Changes were made mostly bv redimensioning arrays to 
allow 10 azimuth steps in the free wake geometry calculations. In the early 
development of CAMRAD.Modl, under certain circumstances, the free wake 
geometry calculations gave an error messages indicating too many "transi- 
tion points-', which caused the program to stop. (For details on transition 
points, see Ref. [7]). In order to allow CAMRAD.Modl a better chance 
of completing the free wake portion of the program without stopping, the 
number of allowable transition points was doubled from 16 to 32. One of the 
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input parameters that the free wake geometry calculation uses from other 
parts of CAMRAD. Modi is the maximum bound circulation. Modifications 
were made to include choices of different circulation options to be used in 
the free wake geometry analysis. The options available are as follows: the 
original maximum bound circulation, the maximum positive bound circu- 
lation, the maximum negative bound circulation, the maximum outboard 
bound circulation, the maximum inboard bound circulation, and the “large 
core" circulation (discussed in the Section 2.17). The choice of option is 
dependent on the problem being explored. 

2.1.2 Sample Case Discussion 

The effect of changing the azimuthal resolution from 15 to 10 for the 
sample rotor in a descent condition is shown in Figures 2.1 (a) and 2.1 
(b). In this figure, contours of local lift coefficient, Ci, are shown over 
the rotor disk. These predictions were made using the original maximum 
bound circulation option in the free wake model. Figures 2.2 (a) and 2.2 (b) 
show the lift coefficient as a function of span for the 15 and 10 azimuth 
step cases for several azimuth locations. The major difference in this case 
compared to original CAMRAD is the azimuth resolution used for the free 
wake analysis. Some small differences are seen; but no systematic study of 
these small differences has been made. However, the primary purpose of 
the modification is to provide higher resolution wake and blade geometry to 
HIRES reconstruction than would be possible with the original CAMRAD 
free wake resolution. 

2.1.3 Code Modifications 

The following changes were made in the free wake subroutines: 

1. Array dimensions were changed as follows in common blocks SQCAL, 


SSPLOT, ? 

tiid SGAM: 


(97) 

becomes 

(145) 

(3,25) 

becomes 

(3,37) 

(6,25) 

becomes 

(6,37) 

(3,25,97) 

becomes 

(3,37,145) 

(25) 

becomes 

(37) 

(6,25,16) 

becomes 

(6,37,16) 
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Figure 2.2: The local lift coefficient as a function 
of higher azimuthal resolution in the free wake n 

14 


2. The array KTR(6,37,16) in the common block SQCAL, already modi- 
fied bv the redimensioning in (1) above, was changed to KTR(6,37,32). 
This change was made in subroutines DCALC, GE0MF1, GEOMF2, 
NWCAL, VSCAL, and WQCAL. 

3. In subroutine NWCAL, the line: 

IF ( KM ( I , J) ,GE. 16) GOTO 490 
was changed to: 

IF (KM(I,J) ,GE. 32) GOTO 490 

4. In the subroutines CHEKR1 and CHEKR2, the line: 

IF (LEVEL .EQ. 2) .AND. (MPSI , GT. 24) GOTO 21 

was changed to: 

IF (LEVEL .EQ. 2) .AND. (MPSI .GT. 36) GOTO 21 


5. A change was made in QCVL to handle instances when the variables 
AL and/or BL are < 0. The following lines were added: 

IF ((AL .GT. 0.) .AND. (BL .GT. 0.)) THEN 
Q = . . . (original line from code) 

ELSE 
Q = 0. 

ENDIF 

6. Two FORMAT statements were modified to allow proper output of 
parameters during use of the DEBUG variable. FORMAT statement 
number 2 was changed such that the variable, H, is output with the 
“13” instead of the “12” format, and the following “3X” was changed 
to “2X” to retain the same field width. Also, FORMAT statement 
number 4 was modified so that the “13” format is used instead of the 
“12” format. These FORMAT changes were made in both GEOMF1 
and GEOMF2. 
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7. Options for different maximum circulations used in the free wake calcu- 
lations are controlled by input parameters OPMXFWG, OPROLLU, 
and IFWLGC in namelist NLTRIM (see Chapter 5). These parame- 
ters were added at the location where, originally, the maximum bound 
circulation was stored in an array for use in the free wake geometry 
calculations. 

2.1.4 Extensions to High Resolution 

The free wake geometry is determined by calculations in the low resolu- 
tion portion of CAMRAD.Modl; whereas, in HIRES, the wake geometry is 
interpolated as needed from the low resolution information. This interpola- 
tion is applied between known, low resolution wake endpoints, as discussed 
in Chapter 3. As such, the modifications discussed here are automatically 
included in the HIRES reconstruction procedure without further coding con- 
siderations. 

2.2 Modification to Allow 90 Degrees of Near 
Wake 

2.2.1 Introduction 

When calculating the near wake portion of the wake influence coefficients 
in the trim solution (see Figure 1.2), the original version of CAMRAD placed 
a limit on the extent of wake behind the reference blade that could be des- 
ignated as near wake. This limit was a function of the number of azimut hal 
and radial resolution being used and was due only to array sizes in the 
code. In order to have the capability to test the effects of an extended near 
wake, a modification was made to the code to allow up to 90 of near wake. 
Originally if using the maximum number of radial stations (MRA — 30), 
combined with 36 azimuth steps (10 steps), one was limited to 20 of near 
wake (KNW =2). Modifications were made so that under these conditions, 
one could use up to 90 of near wake (KNW = 9). Figure 2.3 illustrates the 
modification. The near wake is represented symbolically by a circular arc, 
whereas in the code, the near wake is a vortex lattice model. 
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(a) original C'AMRAI) prediction for a near wake extended 20°, 
MRA=30, KNW=2. Ac = 10° 



(b) C'AMRAD.Modl prediction for a near wake extended 90°, 
MRA=30, KNW=9, A«/> = 10° 

f igure 2.3: Extent of the near wake ilust rated for the original C'AMRAI) 
and C'AMRAD.Modl. 
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2.2.2 Sample Case Discussion 

Figures 2.4 (a) and 2.4 (b) show the contours of lift coefficient using 
20 of near wake and 90 of near wake, respectively. Figures 2.5 (a) and 
2.5 (b) show the same information plotted at a several azimuth stations. 
Though there are not large effects apparent in these plots, there may be 
occasions where the extent of the near wake might become an issue; this is 
an engineering choice that is left to the user. 

2.2.3 Code Modifications 

The common blocks were changed as follows: 

In subroutine FILE.]: 

/WKC1CM/ WKC1(7) ,01(135000) ,CNW1 (29600) 

/WKC2CM/ WKC2(7) ,02(135000) ,CNW2(29600) 

was changed to : 

/WKC1CM/ WKC1(7) ,01(252720) ,CNW1 (972000) 

/WKC2CM/ WKC2(7) ,02(252720) ,CNW2(972000) 

In subroutines VINDCAL1, V1NDCAL2, WKC1INT, WKC2INT, and C'FD- 
WAKE: 


/WKINT/ CINT(3, 72000) ,CNWINT(3, 30000) 

was changed to: 

/WKINT/ CINT(3, 72000) , CNWINT(3 , 1 10000) 

In subroutines VVAKECl, WAKEC2, WAKECl, WAKEN2, VVKCTINT, and 
WKC2INT: 

/WKC1CM/ . . . ,CNW(3, 97200) 

/WKC2CM/ . . . ,CNW(3, 97200) 

was changed to: 


/WKC1CM/ ... ,CNW(3, 324000) 
/WKC2CM/ ... ,CNW(3, 324000) 


18 





(a) ( ) for 10 wake resolution with nea 



(b) Ci for 10 wake resolution with nea 


Figure 2.5: The local lift coefficient as a 
of changing the near wake extent. 
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In subroutines CHEKR1 and CHEKR‘2: 


NWMAX = MRG*MRL*MPSI*MAX0(2 ,KNW+1) 

IF (NWMAX .GT. 97200) GOTO 41 

was changed to: 

NWMAX = MRG*MRL*MPSI*MAXO (10 ,KNW+1) 

IF (NWMAX .GT. 324000) GOTO 41 

With these changes, and using the maximum values of MRA=30, 
MRG=30, MRL=30, and MPSI— 36, the near wake may be extended to 
90 (i.e., KNW=9 in namelist NLWAKE). 

2.2*4 Extensions to High Resolution 

The only high resolution variable affected by this change is the array 
size of the variable CNWINT. This array has been dimensioned to be com- 
patible with the low resolution 90 near wake modification. Other common 
blocks (WKOlCM and WKC2CM) listed in the high resolution subroutines 
WKC1INT, WKC2INT, VINDCAL1, and VINDCAL2, have been dimen- 
sioned so that their size is consistent throughout the code. 

2.3 Modifications for Higher Harmonic Control 
(HHC) of Pitch 

2.3.1 Introduction 

In the past, there has been considerable interest in the concept of Higher 
Harmonic Control (HHC 1 ) as a means to modify the certain aspects of heli- 
copter behaviors, such as vibratory loads and acoustic signatures (ref [8]). 
Recently there have been major experimental programs conducted to study 
the use of HHC to reduce the vibratory loads of forward flight conditions 
and to reduce the BVl noise levels for descent flight conditions (Ref. [9]). 
In an attempt to predict or systematically study the effects of HHC 1 , options 
were added to CAMRAD.Modl to include a fixed, user prescribed HHC. 

In a typical wind tunnel trim case, the blade pitch in the CAM- 
RAD.Modl trim loop is adjusted at the collective (0/rev) and cyclic (1/rev) 
levels until a trimmed solution is obtained. It is desirable to include in the 
trim solution, an additional blade pitch that represents HHC 1 . The HHC 1 is 
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a fixed, open loop quantity added to the existing rigid pitch motion cal- 
culated in the “motion iteration” (see Figure 1.2). Therefore, the pitch of 
the blade will include control inputs (0/rev and 1/rev), blade elastic torsion 
(due to sources such as pure torsion, pitch-bending coupling, etc .), and a 
prescribed HHC. CAMRAD.Modl has two modifications for HHO. These 
two modifications vary only in their generality and inputs. 

2.3.2 HHC Pitch Equations 

CAMRAD.Modl now calculates the blade pitch input from the following 
equation: 


0(t) = #o + 0i c cos(i/’) + 0 U sin(0) + Ohhc(V') (2.1) 

where 0 O is the collective pitch, ffu- is the lateral cyclic pitch, 0 ls is the 
longitudinal cyclic pitch, and 6hhc is the new HHC input. 

Two HHC modifications are in CAMRAD.Modl. The first HHC mod- 
ification (input via namelist NLHHC), made early in the CAMRAD.Modl 
development process, provides only the capability to model 3/rev, 4/rev, 


and Tj/rev HHC via the following equation: 

= term 1 + term‘2 + term 3 (2.2) 

terml = cos( Ir - A<p co a) (2.3) 

Urrn'2 = 0[ at cos (40 - 4 <p lat ) sin(v') (2.4) 

termi = $ lon cos(4i/’ - 4 <p ton ) cos (0) (2.5) 


where 0 C olh&lat<0lon,<!>coll,<i>lati and <t>lon are input values in degrees. Figure 
2.6 illustrates each of the 3 terms in the above equations. In each plot, the 
values of Ocoll , <Piai- an< l 0/on- have been set to zero to demonstrate each 
term; these quantities serve only to phase-shift the waveform. The upper 
plot in Figure 2.6 shows terml in the Equation 2.3. It, illustrates a prescribed 
1 amplitude, 4/rev HHC pitch input. The lower left plot shows tervi'2 in 
Equation 2.4 (solid line) along with the sine wave “envelope” defined by 
the term. The lower right, plot illustrates termS in Equation 2.5 (solid line) 
along with the cosine “envelope” defined by the term. Again, both terms are 
plotted for a 1 amplitude HHC input and the Oi,,i and <i>i on terms serve to 
phase shift each waveform. It can be shown that the above equations can be 
used to generate a pure 3/rev, 4/rev, or 5/rev HHC pitch signal. Equation 
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Figure 2.6: HHC :i Terms" 
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2.2 can also produce “wavelets” that are comprised of 3/rev, 4/rev, and 
5/rev HHC components. 

In order to have a more flexible HHC input, and to allow the possibil- 
ity of a form of Individual Blade Control (IBC), the second, more recent 
modification (input via namelist NLHHC2), uses a truncated Fourier series 
to represent the HHC 1 pitch. The HHC equation for this modification is as 
follows: 


12 

^HHc(t) — &hhc\ o + cos (n0) + B n sin (nw)) (2.0) 

n— 1 

where A n ,B n , and 0hh(\o are input values in degrees. The current max- 
imum allowable number of IIHC input harmonics is twelve. With Equa- 
tion 2.0, any HHC 1 waveform may be input approximately by a twelve term 
Fourier series. It should be noted that zeroth and first harmonics are usually 
redundant inputs since the rotor in CAMRAD.Modl, in a typical wind tun- 
nel scenario, is trimmed by adjusting the collective and first harmonics of 
pitch. If this is the case, the trimmed collective and cyclic pitch values will 
merely compensate for these input HHC values. They have been included 
here for completeness and are normally always equal to zero. 

Also now included in the CAMRAD.Modl code are additional motion 
terms due to HHC 1 pitch rate and HHC 1 pitch acceleration as calculated bv 
the following equations: 


0HHC = X(l + DPS) (2.7) 

0HHC = A(1 + DPS) 2 + X(DDPS) (2.8) 

where DPS and DDPS are internal CAMRAD.Modl quantities used to 
account for hub and shaft motion, and A is the azimuthal derivative of either 
Equation 2.2 or 2.(5 and A is the second azimuthal derivative of Equation 2.2 
or 2.(5. depending on the HHC model being used. 

2.3.3 Sample Case Discussion 

Figures 2.7 (a) and 2.7 (b) show the 10 azimuth case without HHC 1 
and with the inclusion of a 4/rev, 1 amplitude HHC pitch which is a pure 
cosine wave starting at t = 0 . Figures 2.8 (a) and 2.8 (b) show the same 
information plotted radially at several azimuth locations. In the HHC case, 
a 4/rev pattern can be seen in the loads due to the 4/rev pitch input (Figure 
2.7 (b)). Since the HHC 1 is included in the blade motion, the 4/rev loading 
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is not necessarily in phase with the HHC pitch input. This can be seen 
in Figure 2.7 (b); the 4/rev loading is not a pure cosine wave starting at 
= 0 as is the HHC input. This case demonstrates that HHC can greatly 
impact the loads for a given rotor. 

2.3.4 Code Modifications 

In subroutines INPTR1 and INPTR2, changes were made to read in both 
HHC' option parameters via namelist NLHHC for the first modification and 
NLHHC2 for the second modification. The parameters in both namelists 
are initialized to zero and are converted internally to radians after being 
read. The input values are thus in degrees. The first set of modification 
parameters are saved in t he common block HHCl for rotor- 1 and in common 
block HHC -3 for rotor-2. The second set of modification HHC parameters are 
saved in common block HHC2 for rotor-1 and in common block HHCl for 
rotor-2. The HHC equations are programmed in the subroutines HHCTRM 1 
and HHCTRM2 for rotor-1 and rotor-2, respectively. The pitch, pitch rate, 
and pitch acceleration terms due to HHC are included in the blade motion 
subroutines MOTNB1 and 1V10TNB2 for rotor- 1 and rotor-2, respectively. 

2.3.5 Extensions to High Resolution 

Since the blade motion in HIRES is obtained from the low resolution 
blade motion determined in CAMRAD.Modl, the HHC modifications dis- 
cussed here are automatically included in the HIRES reconstruction. 

2.4 Modifications for Aerodynamic Sweep Effects 

2.4.1 Introduction 

It is widely known that, in fixed-wing aircraft, swept wings have advan- 
tages in reducing the compressibility effects of high speed flight (ref [10]). 
In the case of rotorcraft, the blade tips are traveling at high subsonic Mach 
numbers, and thus encounter compressibility effects. Many rotorcraft man- 
ufacturers are or have been using some form of swept tip rotor blade design. 
The original version of CAMRAD did not apply any models to account for 
planform sweep. To study aerodynamic effects of mildly swept planforms, 
a modification was made to CAMRAD.Modl to model these effects. This 
modification provides a means by which airfoil characteristics (i.e., lift, 
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of HHC inputs. 
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drag, and moment coefficients) calculated in the “circulation iteration” (see 
Figure 1.2) for swept, yawed flow conditions can be analytically related to 
unswept, unyawed conditions. Since tabulated yawed, swept airfoil charac- 
teristics are not normally available, this analytical relation facilitates use 
of available 2-D tabulated airfoil characteristics. Note that these modifica- 
tions are only included in the aerodynamics of the rotor, and therefore are 
included in the aerodynamic forcing functions for the rotor; no modifications 
have been made to alter the blade dynamics to account for effects of swept 
plan forms (that is, assumptions such as a straight elastic axis, etc. are still 
in place). 

2.4.2 Angle of Attack Correction 

Since airfoil data, is tabulated for 2-D unswept sections, it is convenient to 
relate section properties for a yawed, sw r ept planform to those of an unyawed, 
unswept section so that these tabulated tables may be used for yawed, swept 
blade sections. C AMRAD and CAMRAD.Modl already account for yawed, 
unswept, flow effects as follows: 


GH 


C/;2di(an, A It) 
cos 2 (A) 


(2.9) 


C/(n) 


( d;2d( a t2, A/f) 

cos(A) 


( m ( O ) — ( rn ,2 d ( &t 1 ? A// ) 


( 2 . 10 ) 

( 2 . 11 ) 


an = n cos ( A ) cos(A) 


( 2 . 12 ) 


0*2 = ocos(A) 


(2.13) 


w here A is the yaw angle between the flow' and a section perpendicular to 
the span wise reference line of the blade, O/j and a %2 are angles of attack to 
be used in the airfoil table interpolation, A/* is the Mach number to be used 
in t he airfoil table interpolation, a is the calculated angle of attack for the 
2-D section, and quantities with the subscript 2 d are values found by the 
airfoil table interpolation. The airfoil tables are interpolated to determine 
the CiadiCdM, and C m ;zd at the angle of attack, a n (or a t2 for the drag) 
and at the Mach number, M t • Once these 2-D values are known, the above 
equations relate the unyawed, unsw r ept values to the desired yawed, unswept 
values to be used in the analysis. Note that the cosines in the denominators 
of the lift equation arise from dynamic pressure differences between the 
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(b) Swept planform of sweep angle T and flow yaw angle A. 
Figure 2.9: Non-swept and swept, blade planforms. 


yawed and the unyawed flows. The single denominator of the C\i equation 
results from the assumption that the total drag is in the yawed flow direction 
(a cosine factor in the numerator has canceled out one of the cosine factors 
in the denominator). The cosine factors in the angles of attack, a t \ and 0 * 2 , 
used in the airfoil table interpolation come from two sources. One cosine 
comes from the angle of attack difference between the yawed, unswept and 
unyawed, unswept flows. This is the origin of the cosine in the angle of 
attack in the 072 equation above. The second cosine in the a t 1 equation 
above arises from the "'swept wing equivalence assumption'' (Ref. [11]). 
This cosine accounts for the difference in lift curve slopes between yawed, 
unswept and unyawed, unswept flows. In the context of Ref. [1] and [11], a 
"swept wing” refers to the entire, straight blade, being angled (skewed) to 
the freest ream (Figure 2.9 (a)). 

To model swept planform effects on the aerodynamics, the "swept, wing 
equivalence assumption v is modified to include the local sweep angle effect 
on the angle of attack and Mach number used in the table interpolation. 
This effect does not affect the dynamic pressure portion of the equations 2.9 
- 2.13 above (/.e.. the denominators remain unchanged) since the modifi- 
cation merely relates the yawed, swept properties to the unyawed, unswept 
properties. 

The angle of attack modification is made solely to relate the yawed, 
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swept flow to an equivalent unyawed, unswept flow. This technique is used 
so that the airfoil tables for a 2-D airfoil section (perpendicular to the ref- 
erence span) may be retained and such that the 2-D airfoil tables remain 
independent of the yaw and/or sweep angles. Without this modification, 
t he airfoil tables would necessarily be a function of not only angle of attack 
and Mach number, but also vaw angle and sweep angle. Figure 2.9 (b) illus- 
trates the sweep of a planform and the meaning if yaw angle versus sweep 
angle. The sweep angle actually required for the analysis is the sweep of the 
quarter-chord line if the section is tapered. 

To include modification for the angle of attack, the CAMRAD.Modl 
Equations 2.9 - 2.13 shown above are changed to the following: 


c, 

, . Ct. 2d (a tl , M t ) 

« = 21 W 

(2.14) 


cos z ( A ) 

Cj 

t x Mi) 

(o) = 

(2.15) 

cos(A) 

( m 

•M 

c 

II 

2 

(2.16) 

On : 

- o cos(A + T) cos(A) 

(2.17) 


a t2 = acos(A) 

(2.18) 


where T is the sweep angle of the quarter chord line of the blade with 
respect to the reference span line. The reference span line has the same 
definition as in the original CAMRAD version. Note that due to the "'swept 
wing equivalence assumption”, only one of the cosine terms in the equations 
above is affected. Next, the calculation of the Mach number to use in the 
table interpolation is discussed. 

2.4.3 Mach Number Correction 

As discussed previously, the airfoil tables are interpolated to find a value 
of C'tjd, Cdjd, and C m ,2d at a particular angle of attack and Mach number. 
Then, the equations above are applied to calculate th e C/, Q, and C m used in 
the CAMRAD.Modl analysis. In the previous subsection, the modification 
to the angle of attack used in the table interpolation was discussed. This 
subsection discusses the Mach number modification. 

It is well known from swept wing analysis that there is a compressibility 
relief due to local sweep of a planform. The Mach number modification 
is therefore cast in the form of a compressibility relief term referenced to 
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the original velocity vector at the section. Swept wing theory, along with 
a high aspect ratio assumption, implies that the correct Mach number to 
use in aerodynamic calculations is the Mach number perpendicular to the 
quarter chord line. Thus the Mach number to be used in the airfoil table 
interpolation of 2-D loading must be modified to account for such a relief. 
For the case of no sweep, the Mach number calculation, is as follows (this is 
the form used by the original version of CAMRAD): 


Mt = Ml 




2 \ u$) 


M r 


«oo 


(2.19) 


where M f is the Mach number to be used in the table interpolation, A / _[_ is 
the Mach number in a plane perpendicular to the straight blade (be., the 
straight elastic axis), U p is the velocity perpendicular to the hub plane at the 
current section. I is the velocity parallel to the hub plane at the current sec- 
tion, (ir^ is the sj>eed of sound as calculated internally in CAMRAD.Modl, 
and M corr is a user input constant (input as a function of blade span) that 
could be used to account for any desired constant compressibility relief. For 
a straight, unswept planform, this equation is consistent with the swept wing 
analysis. In order to account for a swept planform, however, this equation 
must be modified. The modification involves calculating the total Mach 
number from the value of M_i as follows: 


M total = ~~~r ( 2 . 20 ) 

cos A 

then, calculating the Mach number normal to the quarter chord (see Ref. 
[10]), as follows: 


M n = Mtotul cos (A + T) (2.21) 

Combining these two equations, the Mach number to be used in the 2-D 
airfoil table lookup of properties is: 


M t = 


M_l cos ( A + T) 
cos(A) 


( 2 . 22 ) 


2.4.4 Sample Case Discussion 

Figure 2.10 shows the lift coefficient contours for 2 cases. Figure 2.10 (a) 
is the 10 azimuth step case with no sweep and Figure 2.10 (b) is the same 
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case except that 30 of aft sweep outboard of r/R = 0.81 is included. Figure 
2.11 illustrates the same information, plotted radially at several azimuth 
locations. For this particular case, the sweep correction model has a very 
small effect. 

2.4.5 Code Modifications 

In subroutines INPTR1 and INPTR2, the following were added: 

REAL SWPL0(30) , SWPHI (100) 

NAMELIST /NLSWP/ SWPLO 

COMMON /SWPCM1/ SWPLO ,SWPHI (added to INPTR1) 

COMMON /SWPCM2/ SWPLO, SWPHI (added to INPTR2) 

A read of the namelist NLSWP was added to INPTR1 and INPTR2 after the 
read of namelist NLBED (to be discussed in a later section). Input variables 
are converted to radians after input and are saved in the common blocks SW- 
PCM1 for rotor-1 and in SWPCM2 for rotor-2. All SWPLO and SWPHI 
quantities are initialized to zero before reading namelist NLSWP. Once the 
low resolution sweep quantities are input, the subroutines INITHR1 and 
1NITHR2 interpolate these quantities to the required high resolution quan- 
tities for use in HIRES. In subroutines AEROS1 and AEROS2, changes are 
as follows: 

add the lines: 

REAL SWPL0(30) ,SWPHI(100) 

COMMON /SWPCM1/ SWPLO, SWPHI (added to AER0S1) 

COMMON /SWPCM2/ SWPLO, SWPHI (added to AER0S2) 

change the lines: 

AEL = ADL+COSLSQ 
AED = ADD*C0SL 
AEM = ADM*COSLSQ 

to the following: 

YAWANG = ACOS(COSL) 

C0SL2 = COS (YAWANG + SWPLO(IR)) 

C0SL3 = C0SL*C0SL2 
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IF (COSL .NE. 0.) THEN 
ML = ML*C0SL2/C0SL 
MD = MD*C0SL2/C0SL 
MM = MM+C0SL2/C0SL 
ENDIF 

AEL = ADL*C0SL3 
AED = ADD+COSL 
AEM = ADM*C0SL3 

2.4.6 Extensions to High Resolution 

The low resolution sweep modifications discussed in this section are also 
applied in HIRES. No additional user input is required for t his to occur. The 
low resolution inputs from the variable SWPLO are internally interpolated 
to the high resolution radial stations input by the user in the array RAEINT 
of namelist NLH1RES. 


2.5 Modification to Airfoil Tables 

2.5.1 Introduction 

The number of and size of the airfoil tables input to CAMRAD.Modl is 
limited. In order to input more C81 airfoil tables and/or more angles/Mach 
numbers per table, several common blocks and several IF statements were 
changed. In the airfoil table preparation program AIRFOIL (see Figure 
1.2), the 081 airfoil tables are read and converted to a “OAMRAD airfoil 
table" (“airfoil. tab” in Figure 1.2) format. This format allows for efficient 
interpolation of the airfoil aerodynamic information during CAMRAD.Modl 
execution. For some rotorcraft, multiple tables need to be read and used. 
This modification effectively increases the number of and size of the input 
airfoil tables. 

2.5.2 Code Changes 

In AIRFOIL and in the subroutine AEROT of AIRFOIL, the common 
block TABLES was changed such that the dimensions of the variables CLT, 
CDT, and CMT were increased from 5000 to 10000. Also, the IF statement: 


IF (NA(NAB)*NM(NMB) *NRB .GT. 5000) ICHECK = 1 



was changed to: 


IF (NA(NAB) *NM(NMB) *NRB .GT. 10000) ICHECK = 1 

in the airfoil preparation program. 

In CAMRAD.Modl, similar changes were made. Common blocks 
A1TABL and A2TABL in the subroutines AEROTl, AEROT2, AET1INT, 
AET2INT, FILER, INPTA1, and INPTA2, were changed such that the di- 
mensions of the variables CLT, CDT, and GMT were increased from 5000 
to 10000. Also, the IF statement: 

IF (NMAX .GT. 5000) GOTO 12 

was changed to: 

IF (NMAX .GT. 10000) GOTO 12 
in subroutines INPTA1 and INPTA2. 

2.5.3 Extensions to High Resolution 

Since the common blocks in subroutines AET1INT and AET21NT were 
changed, no other user intervention is needed for application to the high 
resolution part of the code (HIRES). 

2.6 Modification to Motion Convergence 

2.6.1 Introduction 

When CAMRAD.Modl fails to converge to a trimmed condition, many 
times, an inner loop is the cause of convergence failure. For example, if 
the circulation loop (see Figure 1.2) diverges, most likely the trim loop will 
also diverge. To assist in circulation loop convergence, a lag (relaxation) 
factor is employed in CAMRAD.Modl. However, in some instances, the 
trim divergence is caused by motion loop divergence. In the original version 
of CAMRAD, there is no relaxation factor in the motion loop to assist 
convergence. To help motion convergence in these situations, a relaxation 
factor was added inside the motion loop. This relaxation factor was added 
to the rotor forcing function in order to make the trim convergence more 
robust. The relaxation factor is a user specified factor input to linearly lag 
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the rotor forcing function between successive motion iterations. The form 
of the relaxation is as follows: 


F{k, t) = F old { A\ V’) Ml- FACTM) + F(Ar, VO * (FACTM) (2.23) 

where / is the forcing function, F 0 i d is the forcing function from the previous 
revolution, FACTM is the user specified motion relaxation factor, k is the 
mode shape index, and V T is the current azimuth location. This equation 
is utilized in subroutines INRTMl and INRTM2 for rotor- 1 and rotor-2, 
respectively. The same relaxation factor is used for both rotors. For FACTM 
= 1.0. Equation 2.23 produces the same result as the original CAMRAD 
motion iteration. 

2.6.2 Code Modifications 

The relaxation factor has been added to the namelist NLTRIM, which 
is read by the subroutine, INPTN. The default is FACTM = 1.0, which 
produces no relaxation in the forcing function, as was implemented in the 
original version of CAMRAD. A common block, FORCCM, was added to 
subroutines IN IT HI, INITR2, INPTN, INRTMl, 1NRTM2, and PRNT. The 
common block variables are OLDF1 (16,36), OLDF2(lG,36), and FACTM. 
Coding was added to subroutines INITR1 and INITR2 to initialize the vec- 
tors OLDF1 and OLDF2 to zero. Coding was added to subroutine INPTN 
to include FACTM in the namelist NLTRIM. Coding was added to the sub- 
routine PRNT to include a listing of the value of FACTM in the “INPUT 
DATA" section of the printed output, if the section is requested. 

2.6.3 Extensions to High Resolution 

Since HIRES does not re-trim the rotor, or update blade motion in any 
way, this modification has no effect in HIRES. 

2.7 ROTONET/WOPWOP Interface 

2.7.1 Introduction 

System noise predictions are frequently used to determine the effects of 
design changes in parametric studies. Since many configurations are eval- 
uated, short computer run times are essential. However, accuracy is also 



needed, which necessitates high-quality airloads. As a part of the NASA 
Langley rotor noise prediction efforts, a method was developed for connect- 
ing the airloads calculations of CAMRAD.Modl to the ROTONET rotor- 
craft systems noise code (Ref. [6]). 

This method includes several new subroutines and input parameters to 
CAMRAD.Modl. Execution of CAMRAD.Modl may yield two ASCII text 
files for use in ROTONET. The first file, named ROTPARAM.DAT, con- 
tains standard ANOPP control statements for defining various PARAME- 
TER inputs to the ROTONET functional modules. The second file, named 
ROTABLES.DAT, has table members which provide rotor aerodynamic and 
dynamic information (normally computed by ROTONET modules LRP, 
RWCL RIN, RRD, and RLD) in the correct form for use by the ROTONET 
source noise modules LRN, RTN, and RBN. The first file is intended to 
be "cut-and-pasted” into a ROTONET input deck which executes LRN, 
RTN, and/or RBN. The second file is a self-contained ROTONET input 
job, and when input to ROTONET, will UNLOAD the table members into 
the CAMROT.WRK library file. This library file is then LOADed into 
the ROTONET input file which executes source noise modules, thereby the 
source noise modules in ROTONET can then utilize the airloads calculated 
by CAMRAD.Modl. 

The CAMRAD.Modl notation for two-rotor vehicles is used; that is, 
"rotor-1" and "rotor-2”. For conventional helicopters, the main rotor is 
rotor- 1 and the tail rotor is rotor-2. For tandem helicopters, the forward 
rotor is rotor- 1 and the rear rotor is rotor-2. For side-by-side rotors, such as 
tiltrotors, rotor- 1 is the starboard rotor and rotor-2 is port rotor. Variables 
with "Rl” and "R2” in the names are for rotor- 1 and rotor-2, respectively. 

In addition to t he ROTONET information, 4 files, two for each rotor, are 
output for use in the rotor tone noise code WOPWOP. The first file for each 
rotor, named WOPWOP-Rl.DAT and WOPWOP-R2.DAT (rotor-1 and 
rotor-2, respectively) contain the WOPWOP input namelist, INPUT. The 
second file for each rotor, named W0PF0RCE-R1.DAT and WOPFORCE- 
R2.DAT (rotor-1 and rotor-2, respectively) contain the vertical and inplane 
sectional forces at the CAMRAD.Modl radial and azimuthal locations. 

2.7.2 Code Modifications 

A variable, NOISFL, was added to the NLCASE namelist in 
CAMRAD.Modl as a switch to turn on/off the output of the RO- 
TONET/ WOPWOP information. If NOISFL = 0, then no information 
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for these programs is output. If NOISFL = 1, then information is out- 
put. NOISFL was also added to the common block CASECM in the main 
program, CAMRAD, and in the following subroutines: FILEE, FILER, 
FILEY, FLUT, INPTN, INPTO, PRNTC, PRNTJ, ROTNET, STAB. 
STABD, STABE, TRAN, and TRIM. If NOISFL = 1, the subroutine 
TRIM calls the new subroutine ROTNET to calculate and output the RO- 
TONET/WOPWOP information. In addition to ROTNET, three new sub- 
routines were added for use by ROTNET: RMTN1, RMTN2, and II AVAR. 
To draw on information calculated already in the main part of the code, a 
new common block. RTNCM, was added to these subroutines: PERFR1, 
PERFR2, PRNTC, and ROTNET. 

2.7.3 Known Caveats 

The use of these modifications is strongly NOT recommended as they 
do not include higher harmonic loads and may produce misleading results. 
These modifications as such have not been exercised nor have they been 
fully tested. 

2.7.4 Extensions to High Resolution 

Since these modifications are intended to output low resolution infor- 
mation for use in other codes, these modifications have no bearing on the 
HIRES portion of the code. 


2.8 CFD (FPRBVI) Interface 

2.8. 1 Introduction 

Use of most CFD codes for rotor problems requires a priori knowledge 
of aircraft trim, rotor dynamics, and wake aerodynamics. The standard 
method for obtaining these quantities is to use CAMRAD. Modi to perform 
the usual trim and performance calculations, and then output quantities 
ready for use in isolated-blade CFD codes such as FPRBVI. A common 
method for transferring the aerodynamic environment calculations to CFD 
codes is through a “partial” angle of attack table. This partial angle sums 
the effects of all blade motions and fluid velocities, in a lifting-line form, less 
the effect of the reference blade's own near wake computed explicitly in the 
CFD analysis within the computational domain (hereafter denoted by the 
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term “CFD box”; see Figure 2.12). The C'FD code reads this table, then 
uses it to modify the velocity field through which the blade travels. If mod- 
eling BVI events, detailed information about the vortex wake is required by 
the CFD (FPRBVI) analysis. A method to calculate the required vortex 
wake information is implemented using a non-rotating “BVI box” (see Fig- 
ure 2.13) that surrounds the rotor (discussed later). Another method for 
modeling the aerodynamic environment in the CFD analysis is to account 
for blade dynamic motions (including pitch inputs) and wake-induced veloc- 
ities separately. The careful use of these effects allows for a more accurate 
calculation, including such effects as pitch rate. Upon computing the air- 
loads for one revolution or more, the CFD code may output lift and moment 
coefficients. The file containing these coefficients can then be used by CAM- 
RAD. Modi to modify the airloads distribution (lift coefficient only) used in 
computing the aircraft trim, wake response, and rotor dynamic response. 

This work is based in part on work performed by industry on contract 
to NASA Langley. Numerous updates, corrections, and features have been 
added bv NASA Langley to improve the quality and quantity of information 
provided to FPRBVI (see Ref. [4]). Several of these features include (1) a 
new blade-wake coordinate transform for rigid blade cases, (2) a new elastic 
blade motion interface to out put elastic motion information for the FPR BVI 
analysis. (3) a new direct blade motion modeling interface to FPRB\ I- (4) 
a new interface to include a vortex rollup model, a vortex multi-core model 
and (5) a new interface to pass tip and secondary vortex trajectories, multi- 
core core properties, vortex strengths, and vortex locations relative to the 
C'FD and BVI boxes. 

2.8.2 Code Modifications 

The main routine, CAMRAD, calls the subroutine INPTN which reads 
the namelist NLTRIM. In the NLTRIM namelist, the variable OPREAD(2) 
is used as a switch to enable reading of the CFD input namelist, NLCFD, 
after the read of namelist NLROLL. At present, the CFD interface is only 
applicable to rotor-1. If OPREAD(2) = 2, CAMRAD.Modl expects to 
read NLCFD after NLROLL. Also, if OPREAD(2) = 2, the subroutine 
INPTCFD reads the namelist NLCFD which contains the variables OPCFD, 
OPBVI, PHICFD, RDB(6). BDB(6), and OPMOTN. These variables are 
listed and described in Chapter 5. 

Once the variables are read via the NLCFD namelist, they are stored 
in the common block CFDDATA. This common block has been added 
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Figure 2.12: CFD Box definitions used in CAMRAD.Modl and FPRBVI. 
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Figure 2.13: BYI Box definitions used in CAMRAD.Modl and FPRBVI 


starboard 
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to the executive program, CAMRAD, and to the subroutines AEROF1, 
AERBED1, BVIBOX, CFDAERO, CFDBOX, CFDWAKE, GEOMBVI, 
INPTCFD, ROTNET, TRIM, and WAKECl. 

If the variable OPCFD = 1, subroutine TRIM calls the subroutine 
WAKECl at the end of the trim loop, such that it will recalculate the influ- 
ence coefficients excluding all wake elements not in the CFD box. In order to 
test if a particular wake segment is in or out of the CFD box, WAKECl has 
been modified such that calls to the subroutines VTXL, VTXL2, and VTXS 
include arguments to enable or disable testing based on the type of segment. 
For example, the CFD box only extends around the reference blade; there- 
fore, elements from other blades need not be tested for being in or out of 
the CFD box. (That is, the CFD box includes only wake elements in the 
near wake of the reference blade since the near wake is included implicitly in 
the CFD code.) Both subroutines, VTXL and VTXL‘2, call the subroutine 
CFDBOX, which does the actual vortex segment testing for the inclusion 
in the CFD box. Upon return from WAKECT, subroutine TRIM calls the 
subroutine, CFDWAKE, which calculates a velocity parallel and perpendic- 
ular to the hub plane for each blade section using the recalculated influence 
coefficients. This newly calculated velocity is subtracted from the previ- 
ously stored "dull 1 ' velocities parallel and perpendicular to the hub plane at 
the blade section to yield the “partial velocities”. These new partial veloc- 
ities are used to calculate the required “partial angles of attack”, a p , as a 
function of radius and azimuth. These a p (r, ^) are written to a file called 
ALPHAP.DAT to be read by an external CFD code such as FPRBV1. 

Including only the partial angles in the CFD code is appropriate when 
there are not BVI events. This is because in “non-BVT” flight conditions, 
the vortices are sufficiently far from the blades that usage of only the partial 
angles is sufficient. However, in flight conditions were there are significant. 
BVI events, it is necessary to model the blade vortex passages in a more ac- 
curate manner. For this purpose, the full potential rotor code, FPRBVI, has 
a method for directly computing the down wash at the blade due to tip vortex 
segments generated bv an external free wake model. CAMRAD. Modi has 
an option to add the tip vortex segment information to the ALPHAP.DAT 
file in the form of a wake table containing vortex segment endpoints and 
strength. In addition to the option OPCFD = 1, if the variable OPBVI 
= 1, an additional test is done on each tip vortex segment (via a call to 
the subroutine BVIBOX from the subroutine CFDBOX), in order to deter- 
mine whether or not that vortex segment is in the BVI box. If the segment 
is inside the BVI box, the velocity contribution of the segment is also re- 
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moved from the “partial” inflow as is done if a segment were in the CFD 
box. The velocity contribution is removed to avoid “double-counting” of 
the vortex influence since the velocity due to the vortex will be included 
in the CFD code. This tip vortex wake table includes all vortex end point 
locations, strengths of each endpoint, and a flag used by FPRBVT to deter- 
mine which elements to use in its own BV1 calculations. In addition, the 
tip vortex wake information (including the secondary vortex) is written to 
a file called “ALLWAKE.DAT”, which tabulates the tip vortex trajectory, 
core properties, strength and location relative to the CFD and BVI boxes. 

If OPMOTN — 0, the tip vortex wake table is written out with segment 
endpoint position defined relative to the flapped blade position to account 
for the effect of flap displacement on blade-vortex miss distance. This is 
accomplished by subroutine GEOMBVI. If OPMOTN = 1, the tip vortex 
wake table is written out with segment endpoint position defined relative to 
the unflapped blade posit ion. The effect of flap displacement on blade-vortex 
miss distance must then be modeled directly in the FPRBVI calculations 
using the information in the MOTION.DAT file. 

To allow direct modeling of blade motion, both rigid and elastic, in 
subsequent CFD calculations, rigid blade motion harmonics and elastic cor- 
rections are written to the MOTION.DAT file. This file contains the rigid 
flap, rigid lag, and rigid pitch motion harmonic coefficients as well as the 
additional blade elastic flap and pitch deflections needed at each azimuth 
and radial station to reconstruct the blade position. (At present, the elastic 
lag is not included in the output of this file.) The total blade flap deflection 
fitotal a f a particular azimuth and V’ can be reconstructed from the rigid flap 
deflection, fi r (ij'), plus an elastic correction tc). 

fitotal{r, </’) = A-(^) + v ') (2.24) 

where 

fir (VO = fio + Ac cos t + fils sin 0 (2.25) 

If OPMOTN = 1, the partial angle-of-attack table written to the AL- 
PIIAP.DAT file is replaced with a table of the wake-induced partial inflow 
at the blade, as a function of radius and azimuth. FPRBVI can then use 
this wake induced velocity, the shaft angle, the blade rigid pitch inputs, 
and the rigid flap motion to reconstruct the aerodynamic environment ex- 
perienced by the reference blade. This method of modeling the flow with 
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FPRBVI allows modeling of phenomena such as pitch rate effects that are 
not contained in the original method using a p information. 

In addition to the output file ALPHAP.DAT, the subroutine CFDWAKE 
creates a file named CAMAERO.DAT. This file is one of two files needed 
if the option to rerun CAMRAD.Modl using the externally generated lift 
coefficients. The other file needed is CFDAERO.DAT which is generated by 
the CFD analysis. The option to rerun CAMRAD.Modl using a combina- 
tion of externally generated lift coefficients and internally calculated values 
for parts of the blade not included in the external calculation, is begun by 
choosing the option OPCFD = 2 for the rerun. This variable choice forces 
CAMRAD.Modl to read the files CAMAERO.DAT and CFDAERO.DAT 
via a call to the subroutine, CFDAERO. Both files contain lift coefficients 
as a function of radial and azimuthal location. These data are stored in t he 
arrays CLOLD and CLEXT, respectively, which are in turn stored in the 
common block EXTAERO. The common block EXTAERO has been added 
to the subroutines AEROF1, AERBEDl, and CFDAERO and contains the 
arrays CLOLD(30,36) and CLEXT(30,36). If OPCFD — 2, during the trim 
process, the lift coefficient is replaced by the formula: 

Cl{r, V’) = Of r, V’) - </’) + C/,, T/frn „/(r, V’) (2.20) 

where C7(r, t ? ) i* s the currently calculated lift coefficient, C f / t0 / r /(/\ 0) is 
the lift coefficient from the previous CAMRAD.Modl execution, and 
Cl, external i r ^ V') is the lift coefficient from the external CFD analysis. At 
the end of the run, the CFD code is rerun, if desired, and the loop repeated 
in a. open loop manner as the user deems necessary. 

2.8.3 Known Caveats 

It should be noted that when using the Trailed Wake Algorithm (TWA) 
in the Indicial Aerodynamics, partial angles computed will be erroneous 
since the TWA, as implemented, can not truncate the near wake to exclude 
vortex elements inside the CFD box. Thus, it is not recommend at this time 
to execute the indicial aerodynamics option and the CFD interface option 
together. 

2.8.4 Extensions to High Resolution 

These modifications are intended to output low resolution information 
for use in the external CFD codes. Therefore, they have no bearing on the 


45 



HIRES portion of the code. 


2.9 Modification for Tunnel/ Fuselage Corrections 

2.9.1 Introduction 

A modification was made to CAMRAD.Modl to include a tunnel and/or 
fuselage correction model in the aerodynamic calculations. Actually, the 
information provided to the code through this modification does not inher- 
ently assume a fuselage or tunnel wall; it could be any body introducing a 
velocity field near the rotor. The correction involves including an additional 
velocity distribution and an additional wake distortion due to the influence 
of wind tunnel walls or due to a fuselage body (or due to any object that 
produces a steady velocity distribution and wake distortion for the rotor in 
question). The additional velocity distribution is superimposed on the wake 
induced velocity at the rotor in the circulation iteration (see Figure 1.2), 
which is used to calculate airloads. The additional w T ake distortion is added 
to the wake distortion used to calculate the wake influence coefficients in 
the trim iteration (see Figure 1.2). 

The additional velocity over the disk and additional wake distortion are 
fixed, user input quantities read from two files for each rotor being used. 
The first file contains the additional velocity for all radial and azimuthal 
stations and one additional velocity at the hub. The additional velocity 
values (read from the file) should be non-dimensionalized by the rotor tip 
speed. Depending on the user's method of running CAMRAD.Modl, the 
input velocity may be a “total additional velocity” or a “delta additional 
velocity". If the delta form is used for an isolated rotor in a wind tunnel, for 
example, the mean velocity correction, externally calculated and removed 
from the velocities in the velocity file, may be included as an alteration to the 
advance ratio and to the shaft tilt. The velocity file would therefore consist 
of an incremental velocity distribution at the rotor disk, due to a tunnel 
and/or fuselage presence, excluding the mean additional velocity. The other 
method, the total additional velocity method, involves not modifying the 
flight condition ( z.e., the advance ratio and the shaft tilt), but instead, 
including the mean velocity correction in the additional velocity distribution. 
Traditionally, the former method is used since, in a typical wind tunnel test 
of a rotor, a shaft tilt correction is used to offset wind tunnel wall effects. 

The second input file for each rotor contains an additional wake geometry 
distortion. This additional distortion is used in the code as a correction to 
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the wake geomet ry to account for the effects of tunnel walls and/or a fuselage 
body. The additional distortion is input as an incremental distortion at each 
wake endpoint location. Since the tip vortex wake endpoints are identified by 
a .11 azimuth index and an age index inside CAMRAD.Modl, the additional 
wake distortion vector is also identified, and input, in the same manner. The 
additional wake distortion vector is calculated externally by multiplying a 
local additional velocity at the current wake endpoint location by the time 
needed for the rotor to advance through one time (azimuth) step. Since this 
provides an incremental distortion at the wake endpoint location, the total 
additional wake distortion a wake endpoint is calculated internally as the 
sum of the incremental distortions along a path from the creation time of 
the vortex to the current age of the vortex. 

The wake geometry of a. particular tip vortex segment is determined in 
CAMRAD.Modl using either the rigid wake or the free wake method. In the 
calculation of the influence coefficients, the total additional wake distortion 
vector, as discussed previously, is added to the current tip vortex position 
vector of the wake geometry. This vector is then used in the influence coef- 
ficients calculation lor the particular vortex segment in question. Since the 
tip vortex position vector is calculated as needed, not stored, the total ad- 
ditional wake distortion vector is calculated as needed also. For this reason, 
this procedure is duplicated in the wake geometry subroutines (GEOMP1 
and GE0MP2) and in the CFD interface BYT wake geometry calculation 
s u b r o u t i n e ( G EO M B V I ) . 

2.9.2 Equations 

CAMRAD.Modl calculates the position vector for a wake endpoint as 
needed. When the position vector is calculated, an additional w r ake distor- 
tion vector is also computed. The additional distortion vector is added to 
the current wake endpoint according to the location of the wake endpoint in 
the (t\<3>) coordinate system, where 0 is the azimuth at which the endpoint 
was released from the blade and </> is the current age of the endpoint under 
consideration. Therefore the wake location f wa ^ is found as follows: 

1'wake(ty\4 > ) — 1* wake, be fore ( 0 ) T ^ r wake ( 0 ? 0 ) ( 2 . 2 /) 

where r wa ^^ f .j orf {%\ <p) is the vortex endpoint location vector before the 
tunnel/fuselage correction and Ar0 a .kr(0, <f>) is the additional distortion due 
to the current correction. 
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At a radial station, r^, and an azimuth station, ip , during the calculation 
of wake induced velocity in subroutines WAKEN1 and WAKEN‘2, an addi- 
tional velocity is superimposed on the previously calculated wake induced 
velocity at that point: 

Vind 4 otai{n, t) - VinAn, t) + f) (2.28) 

where Al7(n>, V ? ) is the additional velocity from the current correction. As 
discussed earlier, this velocity could be an incremental (i.e., not includ- 
ing the average additional velocity) or a total velocity (i.e., including the 
average additional velocity), depending on the user's method of running 
CAMRAD.Modl. 

2.9.3 Code Modifications 

A parameter w r as added to CAMRAD.Modl to control usage of the tun- 
nel/fuselage correction model. The parameter OPWFCOR was added to 
the namelist NLHIRES (discussed in a later chapter) and is saved in the 
common blocks INTAZ and INTAZ2, for rotor-1 and rotor-2, respectively. 
If OPWFCOR — 0 the correction model is not used. If OPWFCOR = 1, 
the 2 input files, discussed earlier, are read by the subroutines INPTR1 and 
INPTR2 for rotor- 1 and rotor-2, respectively. The following code, contain- 
ing the required input file formats was added to INPTR1 (the same coding 
applies to 1NPTR2): 

IF (OPWFCOR .EQ. 1) THEN 
DO 799 N = 1,12 

799 READ (11,810) DUMMY 

DO 800 J = 1 ,MPSI 
READ (11,810) DUMMY 
DO 800 K = 1 ,MPSI*4 

800 READ (11,820) DC0RR(I , J ,K) , I = 1,3) 

810 FORMAT (A30) 

820 FORMAT (13X ,3(3X ,F12 . 6) ) 

DO 830 N = 1,3 
830 READ (12,810) DUMMY 

READ (12,860) (VCORRH(I) , I = 1,3) 

DO 850 J = 1 ,MPSI 
READ (12,810) DUMMY 
DO 850 IR = 1 ,MRA 
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850 READ (12,860) (VC0RR(I , IR, J) , I = 1,3) 

860 FORMAT ( 19X,3(3XF12 .6) ) 

ELSE 

ENDIF 

where MPSI is the number of azimuth stations (usually 36), DCORR is the 
(input) additional wake distortion, VCORR.H and VCORR are the addi- 
tional velocities at the hub and over the rotor disk, respectively. (VCORR 
is the input additional velocity over the disk.) Note that the unit numbers 
are 11 and 12 for the distortion and velocity files, respectively. For rotor-2, 
the above coding is the same in INPTR2, except the unit numbers 11 and 
12 are replaced by the unit numbers 21 and 22. The above coding provides 
the necessary input file formats. Also note t hat additional wake distortion 
information must be provided for 4 wake spirals (see the “DO 800 ...'* loop). 
The data for the corrections is stored in the common blocks WFCORR 1 and 
WFOORR2 for rotors 1 and 2, respectively. 

In the subroutines WAKEN 1 and WAKEN2, the additional velocity is 
applied by the following additional coding: 

IF (OPWFCDR .EQ. 1) THEN 
DO 777 I =1,3 
DO 777 KR = 1 ,MRA 

777 VIND(I ,KR,L) = VIND(I,KR,L) + VC0RR(I ,KR,L) 

ELSE 

ENDIF 

where VI ND is the wake induced velocity at the rotor disk. 

Since the position vectors of the wake are always calculated as needed, 
the total additional wake distortion vector is also generated as needed. The 
total wake geometry is calculated in three places in the low resolution part 
of CAMRAD.Modl; once in the influence coefficient calculations (WAKECl 
and WAKEC2), once in the geometry print subroutines (GEOMPI and GE- 
OMP2), and once in the CFD interface BVI wake location calculation (GE- 
OMB VI). In WAKECl and WAKEC2, the following was added: 

REAL DSUM(3) 

IF (OPWFCDR .EQ. 1) THEN 
JDCNT = LM - K1 
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IF ( JDCNT .GT. 0) GOTO 668 
JDCNT = JDCNT + MPSI 
GOTO 664 

665 CONTINUE 

DO 670 I = 1,3 
DSUM(I) = 0. 

DO 671 KAGE = 1,K1 

DSUM(I) = DSUM(I) + DCORR( I, JDCNT, KAGE) 

671 CONTINUE 
670 CONTINUE 

DO 666 1=1,3 

RSECTE(I) = RSECTE(I) + DSUM(I) 

666 RTET(I) = RTET(I) + DSUM(I) 

ELSE 

ENDIF 

where DSUM is the total additional wake distortion, RTET is the position 
vector of the tip vortex, RSECTE is the position of the secondary tip vortex. 
In the subroutine GEOMBVI, the same coding is used as above, except the 
variable K1 is renamed K, and the variable RTET is named RTV. 

In GEOMP 1 and GEOMP2, code was added to include the additional 
distortion, and to allow' for printing the wake at a resolution higher than 
10 azimuth steps. In order to accomplish this, the wake geometry and the 
distortion are linearly interpolated to a higher resolution by the following 
code modification (the high resolution modifications will be discussed in a 
later section): 

REAL DCQRHI(3) ,DC0RL0(3) 

IF (0PWFC0R .EQ. 1) THEN 
JDCNTL0 = JREFL0 - K0UNT 

664 IF ( JDCNTL0 .GT. 0) GOTO 665 
JDCNTL0 = JDCNTL0 + MPSI 
GOTO 664 

665 CONTINUE 

JDCNTHI = JREFHI - K 
668 IF (JDCNTHI .GT. 0) GOTO 667 
JDCNTHI = JDCNTHI + MPSI 
GOTO 668 
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667 CONTINUE 
ENDIF 

DO 801 I = 1,3 
IF (OPWFCOR .EQ. 1) THEN 
IF (KOUNT .EQ. 0) THEN 
DCORLO(I) = 0. 

ELSE 

DCORLO(I) = 0. 

DO 811 KAGE = 1, KOUNT 

DCORLO(I) = DCORLO(I) + DCORR(I , JDCNTHI ,KAGE) 

811 CONTINUE 
ENDIF 

DCORHI(I) = 0. 

DO 816 KAGE = 1,K 

DCORHI (I) = DCORHI(I) + DCORR( I , JDCNTHI , KAGE) 

816 CONTINUE 
ELSE 

DCORLO (I) = 0. 

DCORHI (I) = 0. 

ENDIF 

RWT(I) = (RWA(I) +DCORLO(I) ) *FACT + 

(RWD(I) +DCORH I ( I ) ) * FA CT2 
RWSO(I) = RWKB(I) *FACT1 + RWKE(I) *FACT2 
801 RWSI(I) = RWKC(I) *FACT1 + RWKF(I) *FACT2 

(extensive code added here to account for 
vortex rollup and spinup - see source code for 
details . ) 

2.9.4 Sample Case Discussion 

The effect of the tunnel and fuselage corrections on the lift coefficient 
are shown in Figure 2.14. The body used in this case is the sting/fuselage 
depicted in Ref. [9]. The tunnel used in this case is the I)NW tunnel 
also depicted in Ref. [9]. Figure 2.14 (a) shows C\ contours without the 
tunnel/fuselage correction while Figure 2.14 (b) shows contours of the same 
case, except the tunnel/fuselage correction model is used for a typical rotor 
wind tunnel “fuselage” in an open wall wind tunnel. Figure 2.15 is the same 
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information, again plotted radially at several azimuth stations. For this case, 
the average shaft tilt plus an incremental velocity change method was used. 
To account for the average shaft tilt correction, the shaft tilt was changed 
from 5.3 aft for the case without the correction applied to 4.233 aft for t he 
case with the correction applied. Thus the average shaft tilt portion of the 
correction was 1.067 forward tilt. The average change in advance ratio for 
this case was negligible {^correction < -001) and thus was not applied. 

2.9.5 Extensions to High Resolution 

During the wake influence coefficient calculations in the high resolution 
part of the code, the wake posit ion is interpolated from known quantities. To 
avoid unnecessary interpolation, the additional velocity over the rotor disk 
is initially linearly interpolated to the high resolution radial and azimuthal 
locations. This interpolated value is superimposed on the wake induced ve- 
locities in the subroutines VINDCAL1 and VINDCAL2 in the same manner 
as was done for the low resolution wake induced velocity in WAKEN1 and 
WAKEN2. 

In the high resolution wake influence coefficients calculation, the wake 
distortion vector at a higher resolution is interpolated from the low resolu- 
tion information. To avoid unnecessary interpolation, the additional wake 
distortion is added to the low resolution results before interpolation to the 
high resolution. This addition is done in the subroutines WKC1INT and 
WKOINT in a manner similar to that done in the low resolution subrou- 
tines WAKEO and WAKEC2. 

2.10 Low Resolution Loading Output File 

Several subroutines were modified in CAM RAD. Mod] in order to output 
information to a file at the end of the trim iteration. Subroutines LOADR1 
and LOADR‘2 were modified to include calls to the new subroutines PRFILl 
and PRFIL2, respectively. For rotor-1, the output information is written to 
unit number 7 and for rotor-2 the output information is written to unit 
number 8. The output information, is as follows: 

1. Radial stations, ra(r) 

2. Angles of attack, o(/y^) 

3. Mach number, A/(r, 0) 
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Figure 2.15: The local lift coefficient a*s a function of span showing the effect 
of tunnel/fuselage correction. 
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4. Tangential velocity, Ut(r, td) 

5. Perpendicular velocity, U p (t\ ij y ) 

6. Radial velocity, U r (t\ VO 

7. Inflow angle, ^(r, V’) 

8 . L ift coefFi ci en t , Ci ( r , ) 

9. Drag coefficient, Crf(r, c ? ) 

1 0 . M omen t coeffi cie n t , ( f m (r, c' 1 ) 

11. Pitch angle, 6(t\ 0) 

12. Tip flap, /i(V0 

13. Maximum bound circulation, F(<e) 

Except for the ra(r), and P ($) arrays, the following format is used: 

WRITE (XXX, ’(A30)’) YYY 
DO 25 I = 1 ,NRAD 

25 WRITE (XXX, * (12F11.5) ’ ) (ZZZ(I,J) , J = 1,NAZM) 

where XXX is the unit number 7 or 8, YYY is one of the variable names 
listed above, ZZZ is the variable associated with the name YYY, I is the 
radial station index, J is the azimuth station index. The radial station 
output has the following format: 

WRITE (XXX, 669) 

669 FORMAT (IX, 'RADIAL STATIONS J ) 

WRITE (XXX, 670) ( R A ( I ) , I = 1 ,MRA) 

670 FORMAT (1X,10F12.4) 

and the tip flap and maximum bound circulation have the format: 

WRITE (XXX, 888) 

888 FORMAT (IX,’ TIP MOTION - FLAPPING’) 

DO 889 J = 1 ,MPSI 

889 WRITE (XXX, 890) PSI ( J) , ZZZ( J) 

890 FORMAT (2FU.5) 

where XXX and .1 have the same meaning as before, PSI is the azimuth 
location, and ZZZ is either /J or T. 



2.11 Wake Geometry and Blade Position Output 
Files 

2.11.1 Introduction 

A modification to several subroutines was made in order to write to out- 
put files the wake geometry and blade position at an azimuthal resolution 
that is not necessarily the same resolution as used in the CAMRAD.Modl 
trim procedure. The purpose of printing this information into the output 
files is for plotting of the wake geometry and blade position information with 
user supplied plotting programs. For example, if the trim was performed 
at 10 azimuth steps, the wake and blade position may now be written out 
at, say, 5 steps. The modification involves changes to the subroutines GE- 
OMP1 and GEOMP2. The changes are identical in these subroutines except 
for the unit numbers and the common blocks used. If the wake geometry 
print variables in namelist NLLOAD are being used, the wake geometry 
and blade position information are written to files at the resolution defined 
by the variable MPSIWGP in the namelist NLHIRES. This variable is the 
number of azimuth steps per revolution for the wake geometry printout. 
For example, if MPSIWGP = 72, the wake geometry and blade position 
would be written to files at a 5 per azimuth step resolution (n'.e., (360 per 
revolution )/(72 steps per revolution)). 

2.11.2 Code Changes 

The tip vortex wake geometry information at the high azimuthal res- 
olution is derived from the low resolution information by linear azimuthal 
interpolation. The coding added toGEOMPl to output requested informa- 
tion is as follows: 

KMAX = MAXO (KFW , KDW , KNW , KRW) 

IF (MPSIWGP .LE. 0) GOTO 2 

DO 777 J = 1, MPSIWGP 

XPSI = FL0AT(J)*360 . /FLOAT (MPSIWGP) 

IZEFt = 0 
XPHI = 0. 

CALL FINDRAD (...) 

CALL GE0ME1 (...) 

CALL GE0ME1 (...) 

PSIRAD = XPSI*PI/ 180 . 
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CALL WAKEB1 (...) 
DO 800 I = 1,3 
RWT(I) = TDUMB(I) 
RWSO(I) = TDUMB(I) 
RWSI(I) = RDUMB(I) 


(extensive coding added to implement rollup 
(and spinup) . 

(tunnel/fuselage corrections added here.) 

WRITE (XXX, 778) IZER.XPHI , J ,XPSI ,RWT(1) , 

1RWT(2) ,RWT(3) ,RSEC(1) ,RSEC(2) ,RSEC(3) 

DO 777 K = l.KMAX 

XPHI = FLOAT (K) *360 . /FLOAT (MPSI) 

KOUNT = K-l 

IF (ABS(l.-FACTl) .LT. 0.001) KOUNT = K 
CALL GE0ME1 (...) 

CALL GE0ME1 (...) 

WRITE (XXX ,778) IZER.XPHI , J ,XPSI , 

1RWT(1) , RWT(2) ,RWT(3) ,RSEC(1) ,RSEC(2) ,RSEC(3) 

777 CONTINUE 

778 FORMAT (2X,I8,F10.2,I8,F10.2,6F12.5) 

where XXX is the unit number 13 for rotor-1 (GEOMP1) or unit 23 for 
rotor-2 (GEOMP2), and the subroutine F1NDAZ is a new subroutine used 
to find the linear interpolation factors for a given azimuth location. For 
use in GEOMP2, the same coding is used, except the calls to GEOME1 are 
changed to calls to GEOME2, and calls to WAKEB1 are changed to calls 
to WAKEB2. Actual arguments to the subroutines denoted bv “(...)” here, 
may be found in the source code. The coding added to GEOMP1 to output 
the blade position information is as follows: 

DO 886 JJ = 1 ,MPSIWGP 

XXPSI = FLOAT (J J ) *360 . /FLOAT (MPSIWGP) 

XXPSI = XXPSI+PI/180. 

CALL WAKEB1 (...) 

DO 886 I = 1,3 
XR00TB(I) = RDUMB(I) 



XTIPB(I) = TDUMB(I) 

886 CONTINUE 

DO 888 JJ = 1 , MPSIWGP 

XXPSI = FLOAT( J J) *360 . /FLOAT (MPSIWGP) 

XXPSI = XXPSI*PI/180 . 

CALL WAKEB1 (...) 

DO 888 I * 1,3 
DO 887 II = 1 ,NRAD 

887 XMIDB(I , II , JJ) = MDUMB(I,II) 

888 CONTINUE 

DO 889 I * 1,3 
DO 889 JJ = 1, MPSIWGP 

889 WRITE (YYY ,890) XROOTB(I , JJ) , (XMID(I , II , JJ) , 

1 11=1, NRAD) ,XTIPB(I, JJ) 

890 FORMAT (2X,10F12.5) 

where YYY is t he unit number 14 for rotor-1 (GEOMP1) or unit 24 for rotor- 
2 (GEOMP2). Again for GEOMP2, the calls to WAKEB1 are replaced by 
VVAKEB2 and the actual arguments may be found in the source code. 

2.11.3 Known Caveats 

The wake geometry printer plot does not include any of the modifications 
listed in this section - only the printed output wake geometry file includes 
these modifications. 

2.11.4 Extensions to High Resolution 

The high resolution calculations are independent of these wake geometry 
output modifications. These modifications are simply for the purpose of 
printing the wake geometry into an output file for plotting. 


2.12 Tip Core Size Modifications 

2.12.1 Introduction 

In the early development of CAMRAD.Modl, it was determined that, a 
single-sized tip vortex was not adequate to model certain effects such as blade 
vortex interactions. Several of the early modifications are described in this 
section. First, the original single tip vortex model was modified to allow core 
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size changes as a function of wake age. This single vortex model produced a 
relatively smooth velocity profile radially out from the center of the vortex. 
It was seen in several experimental investigations that the vortex velocity 
profile was, in some cases, quite different from the well known ‘"Scully- type" 
vortex profile. As an attempt to study the effect of changing the velocity 
profiles, and to have these profiles related to the loading on the blade, a 
dual core model was implemented. The dual core model consisted of two 
concentric tip vortices (see Figure 2.16) such that each contained a portion of 
the maximum bound circulation. More recently, a model of the vortex rollup 
process was implemented that supersedes this section. However, if the rollup 
model is not being used, the models discussed here may be used. This section 
is mainly for historical value and was written to document the single tip 
vortex core variable size model and t he dual core model in CAMRAD.Modl 
The vortex rollup model, which is described in a subsequent section, is the 
recommended model. 

2.12.2 Single Core Modifications 

The single tip vortex core model in CAMRAD.Modl was modified to 
allow a variable core size depending on the age of the vortex segment. First, 
the single core model is chosen by selecting the variable OPDCORE = 0 in 
namelist NLHIRES. Then the options for the variable tip core size are chosen 
via the variable OPTVCOR, also in namelist NLHIRES. If OPTVCOR = 0, 
then the original single core model is retained. If OPTVCOR = 11, then a 
step function is applied to the core size. The core size, in this case, starts 
at the size defined by the input variable, CORE(l), in namelist NLWAKE. 
The core size increases to the value RCORINC (in namelist NLHIRES) at 
the vortex age PHIINC (also in namelist NLHIRES). If a core size function 
other than a constant or a step function is desired, OPTVCOR = 12 may 
be chosen. This choice applies a tenth order polynomial in wake age to the 
core size. Any agewise core size variation function may be approximated 
by this truncated series. A summary of single core size choices is listed in 
namelist NLHIRES in Chapter 5. 

2.12.3 Sample Case Discussion 

As an example of the modification, the lift coefficient from the 10 case 
from previous sections is shown in Figure 2.17, along with the same case 
including this modification. Figure 2.17 (a) is the same 10 azimuth case 
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with a constant core size as before. Figure 2.17 (b) is the same case, but 
the single core step function is applied. Figure 2.18 is the same information 
plotted radially at several azimuth stations. For this case, the step increase 
in core size is made after the equivalent of one revolution of wake age. That 
is, the step function is at a wake age of 360 (PHIINC = 360 ), and the 
core radius was increased from CORE(l) = .0212R to RCORINC = .091R. 
The loading for the large core case is seen to be much smoother (Figures 
2.17 (a) and 2.17 (b) ) in the rear of the disk than in the small core case, as 
expected. Even though the core size increase in this example is unrealistic, 
the effect of the modification is shown well. The core size may, in reality, 
vary as a more general function of wake age. In such a case, the function 
may be fitted with the tenth order polynomial and the case may be run with 
OPTVCOR = 12. 

2.12.4 Dual Core Modifications 

A modification was added to CAMRAD.Modl to include a “dual core 
tip vortex" model. This model essentially alters t he velocity profile of a 
vortex core by superimposing two vortex cores of different sizes (/.r., inner 
and outer core) and strengths. In CAMRAD.Modl, a “core model factor 1 ' 
is defined which multiplies the influence coefficients for particular vortex 
segments. This factor can be used as a tool to simulate velocity distributions 
in viscous vortex cores. Several choices, listed in namelist NLHIRES (see 
Chapter 5), may be made when using the dual core model. If the variable 
OPDCORE = 1, then the dual core model will be used in the code. If this 
choice is made, several other parameters in namelist NLHIRES must be set 
by the user. For example, the variable OPCORAC specifies the manner in 
which the inner core size is calculated. If OPCORAC 1 = 0, the inner core 
and outer core sizes are user specified constants. If OPCORAC = 1, the 
inner core size is determined from internal calculations based on the angle 
of attack of the last aerodynamic collocation point along with user specified 
constants. 

Similar to the single core modifications, there is an option to vary the 
core size with wake age for the dual core model. The core size variation 
function choices, for the dual core option, does not include a step function 
option, but does include a tenth order polynomial in wake age for both 
the inner and outer vortex cores. The parameters in this core size function 
determine the usage. There is not a variable, as in the single core model, 
used to turn on or off the core size function. Of the parameters listed below, 
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(a) C{ for 10 wake resolution with constant size single core. 



(b) Ci for 10 wake resolution with single core step function. 


Figure 2.18: The local lift coefficient as a function of span showing the effect 
of wake model specified with single core step function. 
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if all are set to zero, the result is that no dual core agewise size changes 
occur. These parameters are set to values other than zero to produce the 
desired core size changes with wake age. The functions are as follows: 


jVl 


n = / 


inner' + (l 

inner S) T ^ ^ inner t n0 


n — 1 


( 2 . 29 ) 


N2 

t o — ^ outer T dottier ,0 + douter,n4 >n (2..‘{0) 

n=l 

where d mner<u is the variable DCCORFAO in NLHIRES, <Z inner<n is the set of 
variables DC'CORFA(n) in NLHIRES, d f , u/(?r ,o is the variable DCCORFBO 
in NLHIRES, douter.n is the set of variables DCCORFB(n) in NLHIRES, N l 
is the variable NDCOOFA in NLHIRES, N2 is the variable NDCC’OFB in 
NLHIRES. t'i nner is the current inner core radius, and v outer is the current 
outer core radius. 

If the dual core model is used, the strength if the vortex must be split 
between the two such t hat the total strength remains the same. To accom- 
plish this, the inner core influence coefficient is multiplied by a factor, /i, 
and the outer core is multiplied by (1 - ji). The value of fi implemented is 
proportional to the blade chord at the last aerodynamic collocation point on 
the blade, to the angle of attack at the last aerodynamic collocation point 
on the blade, and to the velocity at the last aerodynamic collocation point 
on the blade. The constant of proportionality is a user specified constant. 
The calculations of i are done internally at the vortex leading and trailing 
edges. Also, the value of f) is non-dimensionalized by the maximum bound 
circulation associated with the current vortex segment. The value of li is 
thus: 


a = 


( y 7r cl Tv ) 


( 2 . 31 ) 


where ^ is the user specified constant GAMACST in namelist NLHIRES. c is 
the blade chord at the last aerodynamic collocation point, U is the velocity 
at the last aerodynamic collocation point, a is the angle of attack at the last 
aerodynamic collocation point, and V ma x i* the maximum bound circulation 
at the current azimuth location. 
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2.12.5 Sample Case Discussion 

As an example of this modification, Figure 2.19 shows the lift coefficient 
contours of the standard 10 "case along with the same case using the dual 
core model. In this case, a constant inner and outer core size is used (OP- 
CORAC = 0). No core size “age function'’ was used and a value of 0.1 was 
used for GAMACST. Figure 2.20 shows the same information, plotted at 
several azimuth angles. 

2.12.6 Code Modifications 

Code changes for this modification are made in subroutines WAKRC1 
and WAKEC2. At each occurrence of the single core influence coefficient 
subroutine VTXL, there is now an IF statement involving the variable OPD- 
CORE. If OPDCORE = 0, the single core model is used as described earlier. 
If OPDCORE = 1, the dual core model is used as described previously. The 
dual core model uses a new subroutine VTXL2 to calculate the dual core 
influence coefficients. VTXL2 is identical to the VTXL except the core size 
effect on the influence coefficients is calculated in the calling routine instead 
of inside VTXL2. 

2.12.7 Known Caveats 

1. The dual core model and the rollup model should be mutually exclu- 
sive, but the code does not check to see if only one of these is being 
used. 

2. In the original version of CAMRAD, the wake influence coefficient 
calculations assumed a constant vortex core radius. Since the core 
sizes in this section are variable, strictly speaking, the variation of 
core radius should be accounted for in the integrals that are used to 
derive the influence coefficients. However, as an approximation, these 
effects are neglected and are expected to be small. 

2.12.8 Extensions to High Resolution 

The same modifications are included in the HIRES portion of the code. 
For the single core model, the step function or the tenth order polynomial 
may be applied either in the low resolution only, or in the low and high 
resolution portions, depending on the value of OPTVCOR. The core size 
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expansion functions applied to the dual core model are always in effect in 
both low and high resolution. 

2.13 Modification for Tip Vortex Bursting 

2.13.1 Introduction 

A modification was made to the low resolution part of CAMRAD.Modl 
to include vortex bursting for the single tip vortex core model. When a 
vortex is sufficiently close to a blade, bursting may occur and the vortex size 
and/or strength (circulation) altered. In CAMRAD.Modl, for bursting to 
occur, two criteria must be met: (I) a vortex must either cross the blade, or 
have an endpoint sufficiently close to the blade (i.c., within a user specified 
azimuthal tolerance, PSITOL), and (2) the distance between the blade and 
the vortex at the blade-vortex crossing point must be less than or equal 
to a user specified tolerance, ZTOL (Az/R). To locate the blade-vortex 
crossing point it is assumed that the lead-lag displacement of the blade is 
negligible and that the blade has a no sweep. With these assumptions, the 
blade is modeled as a straight line at each azimuth angle. This line model 
of the blade may then be compared to the location of all wake segments to 
determine which vortex segments cross, or are within an azimuthal tolerance 
of, a blade (see Figure 2.21). 

The procedure for the bursting calculation starts with the initialization 
of the burst vortex core size to the input core size. Since the burst vortex 
strength is calculated as a fraction of the maximum bound circulation, the 
fraction is initialized to unity so that the full strength vortex is used at 
the beginning of the calculation. A blade test parameter is initialized to 
zero. This parameter, which monitors which blade has burst which vortex 
segments, is used in subsequent wake-trim iterations to avoid re-bursting by 
the same blade-vortex encounter. Initializing this parameter to zero means 
“no blade has burst any vortex segment”. 

At each azimuth step in the influence coefficients calculation, the burst 
core size, strength fraction, and blade test parameter are initialized to the 
appropriate value determined by the azimuth angle and wake age. For ex- 
ample, for an azimuth angle of 0 = 100 and a wake age of <f) = 100 , the 
burst core size is initialized to the value at t = 90 with a wake age of 
<!) = 90 . This initialization scheme convects the burst parameters with the 
burst vortex as the burst vortex is convected downstream. At the current 
azimuth location, all vortex segments of the reference blade are tested to 
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Figure 2.21: Schematic of a blade- vortex intersection. 


determine if there is a vortex crossing of any blade on the rotor or if the vor- 
tex endpoints are within the azimuthal tolerance PSITOL. Of the following 
criteria (Equations 2.32, 2.33, and 2.34), if the first criterion is met, there 
is a blade-vortex crossing; if the second or third criterion is met, the vortex 
is within the azimuthal tolerance PSITOL of the blade: 

{V'bladf yie) { tblade Vie.) — 0* (2.32) 

1 y blade ty'Ve | < PSITOL (2.33) 

< PSITOL (2.34) 

where t blade., tie, and tie are the azimuth angles of the line model of the 

blade, of the leading edge of the vortex , and of the trailing edge of the 

vortex, respectively (see Figure 2.21). 

If a blade azimuth-vortex crossing is detected, the intersection point is 
located by first calculating the radial location of the intersection along the 
blade. If the radial location of the “crossing” is less than the innermost 
radial station on the blade or if it is outboard of the blade tip, then the 
“crossing” is ignored since there is no blade segment involved. With the 
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radial station of the intersection known, the vertical coordinate of the the 
blade and of the wake is found at the intersection by linear interpolation 
between known blade and wake coordinates. The difference in these vertical 
coordinates is the blade- vortex intersection “miss-distance”. If the miss- 
distance is greater than the user specified tolerance, ZTOL, the intersection 
does not burst the segment. If it is less than, or equal to the vertical tol- 
erance, the vortex is burst. The bursting changes the vortex core size and 
the vortex strength. Also, the blade test parameter is set equal to unity so 
that in a later wake-trim iteration, the same blade-vortex intersection does 
not generate an additional burst of the same vortex. 

If the vortex segment is within the specified azimuth tolerance, the “in- 
tersection” point is calculated as the closest point on the blade to the closer 
of the leading edge or trailing edge of the vortex segment. Again, if the 
“intersection” point is inboard of the innermost radial station, or outboard 
of the tip, the intersection is ignored. Otherwise, the blade height is deter- 
mined by linear interpolation of known blade information, and the vortex 
height is assumed to be the vertical coordinate of the endpoint that is in- 
volved in the “intersection”. As before, the “miss-distance” is the difference 
in the blade z-coordinate and the vortex z-coordinate. 

Currently, a simple bursting model is used. For a bursting event, the 
core size and strength are multiplied bv the factors, CORMULT and CIR- 
MFLT. respectively. Then the burst core size and strength are saved. As 
the calculation proceeds, the same vortex may be burst again by another 
blade. H owever, the same vortex may not be burst again by the same burst 
event during a subsequent wake-trim iteration. 

2.13.2 Sample Case Discussion 

As an example of this modification, the 10 case without bursting is 
compared to the case with bursting. This example is intended to show the 
effect of no bursting compared to the case were all vortices have been burst. 
This is not t he intended usage of the model, but illustrates that bursting may 
have a significant effect on the predictions. Figure 2.22 shows a comparison 
of lift coefficient contours for these two cases. The case in Figure 2.22 (b) 
uses the bursting model (OPBURST = 1) with only the elements “crossing” 
a blade being burst (PSITOL < 0.0). The vertical tolerance in this case 
is set to a relatively large value (ZTOL = 0.1) in order to dramatically 
show the effect of the bursting model. This value effectively bursts all wake 
elements on the rotor disk. For this forward flight case, all wake elements 
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are within a 0.1R vertical distance of the blades. Figure 2.23 illustrates the 
same information, but plotted as a function of radius for several azimuth 
angles. 

2.13.3 Code Modifications 

Changes were made to INPTWl and INPTW2 to include a namelist 
read of the new namelist NLBIJRST. This namelist is read after NLYVAKE 
for each rotor. (The namelist NLBIJRST parameters are listed in Chapter 
5.) Also, included in INPTWl and INPTW2 is the initialization of the ar- 
rays COREBLO and CIRCBLO. These are the burst core size and strength 
fraction arrays, respectively. COREBLO is initialized to the input core 
size. CIRCBLO is initialized to unity. New common blocks, BURST and 
BURST2, for rotor-1 and rotor-2, respectively, were also added. Changes 
were made to subroutines WAKECl and WAKEC2 to include vortex burst- 
ing. The common blocks BURST and BURST2 were added to these subrou- 
tines. Loops were added to calculate the effects of bursting. One of these 
loops calculates and stores the blade position at all radial and azimuthal 
stations. Another loop initializes the arrays COREBLO, CIRCBLO, and 
INBSTLO at each new reference blade azimuth location as described in 
above. INBSTLO is the blade test parameter. It is an array in the common 
blocks BURST and BURST2 used to prevent bursting of a vortex that was 
burst by the same burst event during a previous free wake-trim iteration. A 
loop was added inside the wake age loop of the influence coefficients calcu- 
lation to determine if a bursting event occurs and if so, to set COREBLO, 
CIRCBLO, and INBSTLO to appropriate values. 

2.13.4 Known Caveats 

1. The circulation burst model (as opposed to the core size burst model) 
is not operational in CAMRAD.Modl. 

2. The burst model (circulation burst or core size burst) has not been 
exercised. 

3. The burst model is implemented for the single tip vortex core model 
only. 

4. The burst model is not implemented for the rollup model. 
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(a) Ci for 10 wake resolution wit hout bursting model. 



Figure 2.23: The local lift coefficient as a function of span showing the effect 
of the wake bursting model. 
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5. The influence coefficient calculations in the original version of CAM- 
RAD were derived for a constant vortex core size. Here, the core size 
changes with wake age and azimuth. The influence coefficients used 
here do not include the effect of a variable core size with wake age or 
azimuth angle. 

2.13.5 Extensions to High Resolution 

The same modifications above apply to the high resolution portion of 
the code. The changes needed for the high resolution application are that 
the arrays COREBLO, CIRCBLO, and INBSTLO are larger and are named 
COREBHI, CIRCBHI, and INBSTHI. The tolerances are the same as in 
the low resolution portion. The blade test parameter INBSTHI insures that 
subsequent interactions in the high resolution far wake do not re-burst, the 
same vortex segment by the same blade. 

2.14 Modifications to Use Input Blade Motion 

2.14.1 Introduction 

A technique was developed to allow the user to input elastic blade motions 
directly into CAMRAD.Modl and use the input blade motion values to 
calculate the resulting aerodynamics. These elastic blade motions may be 
obtained from other analyses or measurements. The current modification 
was tested using measured blade motions, but predicted motions could have 
been used just as well. If the input blade motions are used, the internal 
calculation of blade motions via an iterative harmonic analysis of the rotor 
equations of motion, is overridden. Several new subroutines were included 
to allow measured blade motions to be used in CAMRAD.Modl. 

If using a cantilever blade model (HINGE = 1 in namelist NLRTR), 
MOTNIN-FL and MOTNIN-P read the flap/lag elastic motion and the 
pitch motion, respectively, from the files provided by the namelist variables 
FLAPFILE, LAGFILE, and PITCHFILE (from the new namelist NLMEAS 
discussed later) which contain the elastic blade motion in a TECPLOT (see 
Ref. [12]) format. If using a hinged rotor (HINGE = 0 in namelist NL- 
RTR), the same procedure is followed as above except MOTNIN-FLA is 
used instead of MOTNIN-FL. In addition, the variables FLAPO and LAGO 
define the measured mean flap and lag angles for the hinged blade for use 
in MOTNIN-FLA. The hinged blade motion option has not at present time 
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been tested. If the option HINGE = 2 (articulated blade model) is chosen 
from namelist NLRTR, an error message is printed from subroutine RAMF 
and the code stops. This is because the articulated blade option (HINGE 
= 2) contains no elastic blade motion - only rigid flap and lag. 

All three files, FLAPFILE, LAGFILE, and PITCHFILE, contain mea- 
sured elastic motion at each measured radial station, each in a separate 
TEC-PLOT “zone” (see below). The number of radial and azimuth stations 
read are determined by the variables RIN and MPSIIN, respectively, in t he 
namelist NLMEAS. For example, the file FLAPFILE contains the lines: 

zone 

index, flap 
index, flap 
index, flap 


where “index” is a number associated with the azimuth location of the elastic 
flap value, “flap"*. (The value of “index” is not actually used by the analysis 
- the input parameter MPSIIN is used to determine the number of lines in 
the file and thus the azimuth location of the “flap'' value.) 

Each zone is a list of the measured values (flap, lag, or pitch) at the cur- 
rent radial station and contains the same number of lines that corresponds 
to the number of measured azimuth stations. The flap and lag values MUST 
be input in centimeters. The pitch values MUST be input in degrees. For 
example, the file FLAPFILE contains RIN zones and MPSIIN entries in each 
zone. The first zone in FLAPFILE contains the azimuthal time history of 
the elastic flap for the first measured radial station; the second zone contains 
the azimuthal time history at the second measured radial station, etc. This 
zone convention is repeated for all the measured radial stations. This same 
convention is followed in all three input files, i.e the FLAPFILE, LAG- 
FILE, PITCHFILE files. (Note on interpolation vs. extrapolation: To avoid 
extrapolation, the first zone in each file should be a “dummy” zone contain- 
ing all zeros for a radial station inboard of the innermost radial station used 
in the CAMRAD.Modl analysis. Likewise, the last zone could be a zone 
containing information at r/R = 1.) 

The subroutine IUNI is an interpolation routine used to interpolate data 
onto a desired grid, and SFTMOD is a modified version of the WOP WOP 
Fourier Transform subroutine, SFT. 



The modification is implemented as an option controlled by the variable 
IMODEIN in the namelist NLMEAS. If IMODEIN = 0, the normal CAM- 
RAD. Mod I modal analysis blade motion prediction is used. If IMODEIN = 
l, the input elastic motions are used instead. The maximum dimensions for 
input elastic motions are set to 50 radial stations and 2049 azimuth stations. 

2.14.2 Method for Bending Modes 

In order to use measured blade motions in CAMRAD.ModI, the motion 
must be cast into a suitable form. Since CAMRAD.ModI uses a modal 
analysis to describe the mode shapes of the rotor blade, it is required to 
cast the measured deflections in the form of harmonics of modal amplitudes. 
Since the flap/lag mode shapes are decoupled from the torsion mode shapes 
in CAMRAD.ModI, it is possible to split the problem into 2 parts: (1) the 
flap/lag deflection analysis, and (2) the torsion analysis. First, the flap/lag 
analysis will be discussed. 

The modal equation for the flap/lag deflection is as follows: 

nbm 

~{r. t) = m{r)qk{t) (‘ 2 . 35 ) 

Ar=l 

where, 

5(r, 0) = measured deflection (flap,lag) 

fjk( r ) — bending mode shape of the kth mode (flap, lag) 

(fk(v y ) — modal amplitude of the kth mode 
nbm = number of bending modes 

The unknown in Equation 2.35 is To solve for qk(p)< a system of 

nbm equations is generated at each azimuth location, by multiplying (dot 
product) the modal equation above by a mode shape, fj m (r ), and integrating 
radially over the span to eliminate the radial dependence. The result is a 
nbm X nbm linear system at a given azimuth: 

r\ nbm -\ 

/ i]m(r) • z(r, Ip) dr = V qk{p) / i} m (r) ■ f} k (r)dr (2.36) 

Jo k = i Ji > 

where 1 < in < nbm. The integration has been moved inside the summa- 
tion, since both are linear operators and ^(t ) has been moved outside the 
integration, since it is not a function of radius. In matrix form. Equation 
2.3(j becomes: 
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[A m k]{qk($)} = {B m } (2.37) 

wliere each element of A m k — /J rf m (r) • fjk{r)di\ and each element of B m — 
Jo Vm(f') • t)dr. The solution of Equation 2.37 is as follows: 


Ulk{t)} = [A mk ]- l {B m } (2.38) 

After this step, the modal amplitude that produces the measured de- 
flection, at the current azimuth, is known. This procedure is repeated for 
MPSI azimuth stations. Since CAMRAD.Modl uses the complex Fourier 
coefficients of the modal amplitudes instead of the modal amplitudes di- 
rectly, the modal amplitudes just calculated must be decomposed into its 
Courier components. I he modal amplitude for the k-tli bending mode may 
be represented as a series summed over the number of harmonics used: 


mharm 

9*(V’) = (f^nc cos ( n t) + PnJ sin(ntf’)) (2.39) 

n = 0 

and the coefficients of the series fine and /fiV are related to the complex 
coefficients by: 



fi W 

j-'nc 


; l3 W 

1 tJns 


(2.40) 


2.14.3 Method for Torsion Modes 

The elastic torsion motion is analyzed in a fashion similar to the bending 
motion analysis in the previous section. However, the torsion motion analy- 
sis is a scalar operation rather than a vector operation. The modal equation 
for the torsion motion is: 


ntrn 


o(r, t/>) = ^2 Ca- ('■)/>*• (0) 

A=1 


where, 

0(r, VO — measured torsion deflection 
Ck(r) = torsion mode shape of the kth mode 
Pk($) — modal amplitude of the kth mode 
ntrn — number of torsion modes 


(2.41) 



The solution for Pk(i’) follows the solution for qk{t) in the previous section, 
except the modal equation is multiplied by the mode shapes, C m (r), rather 
that taking a dot product. The resulting equations are: 

I ntm ,1 

/ C m(r)6{r'il')dr =Y^Pk(4') / Cm{rKk{r)dr (2.42) 

fc=i Jo 

where 1 < m < ntm. The matrix form of which is: 


[A mk ]{p k W} = {B m } (2.43) 

The solution of which is: 


{Pk(t)} = [A m k]- l {Bm} (2.44) 

The series for the modal amplitudes is: 


mharm 

Pk(il') = Y Wnc ] fos(n^) + 0 { n k J sin (np)) 

n — 0 

The complex coefficients of which are: 



»(*) 




(2.45) 


(2.46) 


2.14.4 Code Modifications 

A new namelist, NLMEAS, was added for rotor-1 only. (That is, the 
input blade motion modification is currently available only for rotor-1 in 
C'AMRAD.Mod 1.) The subroutine INPTN was modified such that, if the 
namelist NLWAKE is read for rotor-1 (OPREAD(2) = 1 in NLTRIM), the 
subroutine INPTM1 is also called. This new subroutine reads the namelist 
NLMEAS. The variables in namelist NLMEAS are found in Chapter 5. 

The new subroutines, MOTNIN_FL (or MOTNIN_FLA) and MOT- 
NIN_P, read the flap, lag, and torsion measured motion, interpolate the 
motion onto the grid used in CAMRAD.Modl, generate a linear system of 
equations at each azimuth, solve the linear system at each azimuth, Fourier 
analyze the modal amplitudes, and calculate the complex coefficients by 
storing the real Fourier coefficients in the appropriate complex arrays. 

The call to subroutine INRTM1 in subroutine RAMF was changed to 
"CALL 1NRTM1 (IMODEIN)" and the call to MOTNR1 was changed to 
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“CALL M0TNR1 (IMODEIN)”. The subroutine INRTMl was changed such 
that if IMODEIN — 1, the motion for bending and torsion is set to the 
previously calculated complex Fourier coefficients. The line: 


BETA ( JN , JRQW ) =B 


was changed to: 


IF (IMODEIN .EQ. 1) THEN 
BETA( JN , JROW) = BETA( JN , JROW) 
ELSE 

BETA(JN, JR0W)=B 
ENDIF 


and the line: 


THETA (JN, JROW) =B 


was changed to: 


IF (IMODEIN .EQ. 1) THEN 

THETA (JN, JROW) = THETA (JN , JROW) 

ELSE 

THETA (JN, JROW) =B 
ENDIF 


The subroutine INRTMl was changed such that if IMODEIN = 1, 
the motion for bending and torsion is not updated by the internal CAM- 
RAD. Modi motion analysis. The line: 

42 BETA(N1 ,1) = BETA(N1 , I) + DEL*KH 

was changed to: 

IF (IMODEIN .EQ. 1) THEN 
BETA(N1 , I) = BETA(N1 , I) 

ELSE 

BETA(N1 , I) = BETA(N1 , I) + DEL*KH 
ENDIF 

42 CONTINUE 
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and the line: 


44 THETA (N1 , I) = THETA (Nl, I) + DEL*KH 

was changed to: 

IF (IMODEIN .EQ. 1) THEN 
THETA (Nl ,1) = THETA (Nl , I) 

ELSE 

THETA (Nl, I) = THETA (Nl , I) + DEL*KH 
ENDIF 

44 CONTINUE 

2.14.5 Extensions to High Resolution 

The HIRES portion of the code interpolates the needed motion from 
known low resolution information. Thus, there are no additional modifi- 
cations required to the method in order for the high resolution (HIRES) 
portion of the code to work properly. 

2.15 Modifications to Use Input Normal Force Co- 
efficients 

2.15.1 Introduction 

In a manner similar to the user-provided blade motion input, a modifica- 
tion has been made to allow the user to input a file containing normal force 
coefficients. C n , into CAMRAD.Modl. These may be from another analy- 
sis or from measurements. Since measured values of C n at a blade section 
are normally calculated from an integration of measured surface pressures 
around the section, and since accurate dynamic drag coefficients and mo- 
ment coefficients are normally not available from the pressure data, in this 
modification it is assumed that the lift coefficient, C’j, is approximated suffi- 
ciently bv C n - It i s assumed that the drag coefficient, C ( i, and the moment 
coefficient, C m , are sufficiently determined by the normal CAMRAD.Modl 
airfoil table “look-up” procedure. With these assumptions, the user-input 
C' n ' s are interpolated to the resolution needed by the low resolution part of 
the CAMRAD.Modl analysis. These values are then used to replace the 
internally calculated values of ( 
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The input file for C n , called CNFILE in namelist NLMEAS, is the same 
format as described for the FLAPFILE, LAGFILE, and PITCHFILE files. 
The use of this modification is controlled by the variable IAEROIN in 
namelist NLMEAS. If IAEROIN = 0 the C n input file, CNFILE, is not 
used. If IAEROIN = 1 the C n input file, CNFILE, is used. The maximum 
dimensions lor input C f n T are set to 50 radial stations and 2049 azimuth 
stations. 

2-15,2 Code Modifications 

A new namelist, NLMEAS, was added for rotor-1 only. (That is, this 
modification is available only for rotor-1 in CAMRAD.Modl.) The variables 
in namelist NLMEAS are found in Chapter 5. 

If IAEROIN = 1, the subroutine INPTM1 calls the new subroutine 
MSAEROl to read and interpolate the CNFILE data. Then the inter- 
polated data is used in either AEROFl or AERBED1 (depending on the 
aerodynamic model chosen by the variable OPBED in namelist NLBE1) for 
rotor-1). The common blocks AEROMS and AEROIN were added to both 
of these subroutines. 

2.15.3 Known Caveats 

This modification has not been exercised nor has it been tested thor- 
oughly. 

2.15.4 Extensions to High Resolution 

This method has not been implemented in the high resolution calculations 
to date. 

2.16 Indicial Aerodynamics in Low Resolution 
CAMRAD.Modl 

2.16.1 Introduction 

Chapter 4 discusses a computer program called the Indicial Post- 
Processor (IPP), which is used as one of two methods to process the high 
resolution information from HIRES (Chapter .3), in determining the rotor 
blade loading. To study the effects of the indicial aerodynamics models on 


81 



the rotor trim solution, a low resolution indicia] aerodynamics model was 
implemented in CAMRAD.Modl. Since a modification must fit into the 
overall scheme of CAMRAD.Modl, it was not possible to directly use the 
Indicia! Post- Processor code for the low resolution implementation. Also, 
the Indicia! Post -Processor was developed and is tuned for high resolution 
solutions. It has several features that would be inconsistent for use with a 
low resolution solution. The low resolution implementation of the indicia] 
aerodynamics is the topic of this section. 

For simplicity, Figure 2.24 schematically shows an abbreviated flowchart 
of the original CAMRAD trim loop. The outer rectangle represents the 
subroutine, TRIM. The three inner rectangles represent the three stages 
involved in a typical trim solution procedure. These three stages are the 
uniform inflow stage, the rigid wake stage, and the free wake stage. At each 
of these stages, a trimmed rotor solution is obtained and the trimmed solu- 
tion is used to initialize the following stage (if one exists). The subroutines, 
TRIMI arid WAKECl, listed in the rectangles represent the major tasks in- 
volved in trimming the rotor. Subroutine WAKECl calculates the influence 
coefficients for rotor- 1. Subroutine TRIMI calculates the trim solution for 
the current wake stage. 

Figure 2.25 expands subroutine TRIMI to further illustrate its function. 
In this figure, the outer rectangle represents a call to the subroutine TRIMI 
from subroutine TRIM. In TRIMI, the subroutine RAMF iterates on the 
rotor blade circulation and motion with fixed controls until the circulation 
and motion root-mean-square (rrns) change from one revolution to the next 
is less than an input tolerance criterion. The function of subroutine TRIMI is 
to call RAMF which first calculates the forces and moments on the rotorcraft 
with the initial guess for the trim controls, then calculate a "derivative 
matrix'* for use by the Newton-Raphson solution procedure, and finally uses 
a Newton-Raphson method to iteratively solve for the required rotor controls 
for the trimmed flight condition. 

Figure 2.26 expands the RAMF subroutine to further illustrate its func- 
tion. Subroutine RAMF first calculates the blade modes, then iterates on a 
circulation loop which in turn iterates on a motion loop. Inside the motion 
loop, an azimuth loop calculates the blade position and motion with subrou- 
tine MOTNB1. Then the blade aerodynamics are computed using AEROF1. 
Inside AEROF1, there is a radial loop that calculates the aerodynamics for 
all radial stations on the blade at the current azimuth. 

If the low resolution indicia! aerodynamics option is chosen, the effect of 
the opt ion is to replace the azimuth loop in Figure 2.26 with the one in Figure 
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Uniform Inflow Stage 


Rigid Wake Stage 


Free Wake Stage 


Figure 2.24: The basic trim loop of C'AMRAD.Modl. 
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Calculate forces and moments 
with initial control settings 

Calculate derivative matrix by 
incrementing each control 

Free Wake Stage 


Figure 2.25: The TRIM1 loop of C'AMRAD.Modl. 
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RAMF 

Calculate Modes 
Circulation Loop 

Calculate Wake Inflow 

Motion Loop 

Calculate Transfer Function 

Azimuth Loop 

I MOTNB1 
AEROF1 


Span Loop 


AEROS1 


AEROT1 




test motion convergence 
test circulation convergence 


test trim convergence 


Figure 2.26: The RAMF loop of CAM RAD. Modi. 
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2.27. This replacement loop has several differences compared to the original 
loop. First, a new subroutine, TWAL, calculates the near wake effects for 
the reference blade. Second, AEROF1 is replaced by AERBED1 to calculate 
the aerodynamics of the blade. Inside the radial loop in AERBED1, the 
subroutine AEROS1 has been replaced bv AEROS1B, and a new subroutine 
SEPRATEl has been added to calculate the effects of leading edge and 
trailing edge separation. More detail about the subroutines and the changes 
are given below. 

2.16.2 Code Modifications 

First, namelist reads for NLBED were placed in the subroutines INPTR1 
and INPTR2 for rotor- 1 and rotor-2, respectively. Contained in NLBED 
are the new parameters OPBED, HCOR, and ICURV. These and other 
parameters are listed in Chapter 5. 

If OPBED 1 or OPBED2 = 1, the parameter OPSTLL in CAM- 
RAD.Modl is automatically set equal to unity internally, for the correspond- 
ing rotor (1 or 2), regardless of the input value for that parameter. This is 
done since the indicial aerodynamics includes a simple dynamic stall model. 
Therefore, for the indicial aerodynamics to use the 2-D airfoil tables cor- 
rectly, the static stall option (OPSTLL = 1) must be chosen. This static 
st all opt ion in CAMRAD.Modi equates to interpolating the 2-D airfoil table 
information at an unmodified (that is, unmodified by dynamic stall param- 
eters) angle of attack. In addition, the subroutines INPTR1 and INPTR2 
print a. message to standard output reminding the user that the parameter 
OPSTLL has been set equal to unity internally. 

Second, changes were made to the subroutines INPTW1 and INPTW2 
so that if OPBED 1 or OPBED2 = 1, the extent of the near wake (KNW) 
is checked for the corresponding rotor (1 or 2). If KNW is not equivalent 
to 90 , a warning message is printed to standard output. The value of 
KNW is not changed, and the program continues to execute. It is left to the 
user to correct the input parameter. In addition, INPTW1 and INPTW2 
check that the original CAMRAD.Modi near wake is turned off («>., WK- 
MODL(2) thru WKMODL(5) are set = 0 so that the near wake is not 
“double counted"). If not, a warning message is printed to standard output. 
The parameters are not changed and the program continues to execute. It 
is left to the user to set the WKMODL input parameters correctly. 

Third, if OPBED = 1, the TRIM subroutine calls the subroutine 
1NITBED1 and INITBED2 for rotor-1 and rotor-2, respectively. These 
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Azimuth Loop 

MOTNB1 

TWA1 

AERBED1 

Span Loop 


AEROS1B 


AEROT1 


SEPRATE1 


Figure 2.27: Replacement azimuth loop for the low resolution indicia! aero- 
dynamics option in CAMRAD.Modl. 
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subroutines initialize the parameters needed in the judicial aerodynamics 
calculations. 

Fourth, the calls to subroutine AEROF1 from MOTNR1 for rotor-1 and 
AEROF2 from MOTNR2 for rotor-2 were replaced by the following: 


IF (0PBED1 .EQ . 1) THEN 
IF (ITWA1 .EQ. 1) CALL TWAl(JPSI) 

CALL AERBED1 (...) 

ELSE 

CALL AER0F1 (...) 

ENDIF 

and by: 

IF (DPBED2 .EQ. 1) THEN 
IF (ITWA2 .EQ. 1) CALL TWA2(JPSI) 

CALL AERBED2 (...) 

ELSE 

CALL AER0F2 (...) 

ENDIF 

respectively. 

The subroutines AERBED1 and AERBED2 are the indicial aerodynam- 
ics equivalent to the subroutines AEROF1 and AEROF2. The variables 
1TWA1 and ITWA2 are flags that are initially set to zero (internally) so 
that the indicial trailed near wake is not calculated for the case of uniform 
inflow. After the uniform inflow stage, the flags are set equal to unity (in- 
ternally) so that the indicial trailed near w r ake will be included by calls to 
subroutines TWA1 for rotor-1 and TWA2 for rotor-2. 

2.16.3 Subroutine Descriptions 

1. INITBED1 and INITBED2 

Both of these subroutines initialize the parameters needed bv the in- 
dicial aerodynamics subroutines (AERBEDl and AERBED2) at the 
beginning of the trim loop. 

2. TWA1 and TWA2 (Trailed Wake Algorithm (TWA)) 
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These subroutines use a modified version of the trailed wake algo- 
rithm of T.S. Beddoes (see Ref. [13]). These calculate the downwash 
at radial stations, RA, along the rotor blade at a specified azimuth 
station, JPSI, due to a trailed near wake system of vortices extending 
90 behind the reference blade. The parameters WKM()DL(2) thru 
VVKMODL(5) need to be set equal to zero, as discussed earlier. The 
downwash due to the A th vortex segment is determined at the 7th ra- 
dial location. The "current vortex segments’" are located at the radial 
segment endpoints, RAE. Placing the vortex segments at these loca- 
tions is a modification to Beddoes’ Trailed Wake Algorithm (TWA). 
The total instantaneous downwash, W new , at a given radial station is 
then the sum of downwashes from all current vortex segments plus a 
contribution from the previous downwash, exponentially decreased by 
an amount equivalent to convecting the vortex location by one time 
step. The convection geometry has two options, straight or curved, as 
discussed earlier. The basic form of the equations (Ref. [14] and [15]) 
are as follows: 


W new ( 1 i VO — ( A w ( / , A ) + \ w ( /, 7\ ) ) (2.47) 

A 

X W (L K) = X w%ol d{I, K)<-* At + 1.359 D w (2.48) 


y;,(/, K) = v WtM {i, K)e~ 4 * At - 0.359 d w (2.49) 


/ ?(A 1CA 

V 4 nh ) (\ + Y£Ly/2 


VAt 

h 


(2.50) 


C c 


h 2 


h 2 + h c 


(2.51) 


where h is the distance from A' th vortex to /th radial station, h c is 
the input HCOR core size, l is the freestream velocity encountered 
by section, (= /*-(-// sin(V0)i At is the time step (1 azimuth step), 
7(A) is the strength of Ath vortex, <j> is the decay factor depending 
on ICURV, h, and V. 


Usage of the core factor C c in Equations 2.50 and 2.51 is also a modi- 
fication of the TWA of Beddoes. This factor is used to avoid the well 
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known singularity at the center of a potential vortex. This factor is a 
“Scully-type” viscous core factor. Once the instantaneous down wash, 
linen., is known, a velocity deficit function is applied to allow the 
downwash to “build up” , since the vortex does not instantaneously 
produce a downwash field. The resulting downwash, IT app , is then 
used by subroutine AERBEDL TWA2 does the same for rotor-2. The 
deficit function is of the form: 

W app (I , r) = Wnew(I, t) - DEF(I) (2 .52) 

DEF(I) = DEFouWe^r 1 + (W new ( 7, 0) - W ntw {I^ old )) (2.53) 

where c is the local blade chord. 

3. AERBED1 and AERBED2 

Once the instantaneous downwash, W app , is known from TWA1 and/or 
TWA2. all the velocities on the blade are known. The aerodynamic 
loads may then be calculated by AERBED1 and AERBED2 for rotor- 1 
and rotor-2, respectively. The indicia] modifications basically involve 
replacing the angle of attack, with an effective angle of attack that 
includes the features of modified classical unsteady thin airfoil theory. 
That is, the angle of attack will be a quasi-steady angle of attack, 
minus a lift deficiency function in the form of an angle of attack defi- 
ciency function that accounts for the shed wake of an airfoil (see Ref. 

[W]).’ 

The integrated downwash approach used to calculate the quasi-steady 
and of attack in Ref. [14], cannot be used correctly in the low reso- 
lution calculations of CAMRAD.Modl. This integration is replaced 
by an angle of attack calculation that assumes the the downwash is 
comprised of a uniform downwash along the chord due to rotorcraft 
motion, wake inflow, and airfoil motion (except pitch rate) plus a linear 
downwash due to pitch rate motion of the airfoil. A picture illustrat- 
ing this assumption is shown in Figure 2.28. The quasi-steady angle of 
attack needed in this method is the angle of attack based on the values 
of the velocity normal the the chord, U n , the velocity in the chord wise 
direction, f/ c , the velocity due to pitch rate, U p , and the velocity due 
to the near wake, W app , at the 3/4 chord point: 
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Velocity 



Blade section 


(a) uniform downwash along chord due to rotorc.raft motion, 
wake inflow and airfoil motion. 



1/4 chord 

Blade section 


Velocity 

(b) linear downwash due to pitch rate motion of the airfoil. 


Figure 2.28: Downwash components on the airfoil. 


. _l (l- n T ff a n?) T ( ' n \ . _ ,, 

o ;{/4 = tan ( J (2.54) 

Once o :i j A is known, the deficiency functions found in Ref. [14, 15, l(j] 
may be applied to calculate the effective angle of attack. Again, since 
the computation is performed in low resolution, the noil-circulatory 
or impulsive terms defined in the mentioned References are not in- 
cluded. These impulsive components account for the apparent mass 
terms seen in classical unsteady aerodynamics and are high resolution 
effects. Thus including these in the low resolution sections would not 
be consistent with the assumption of low resolution analysis. 

After computing the angle of attack, AERBEDl and AERBED2 call 
the subroutines, AEROS1B and AEROS2B, respectively. These sub- 
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routines, once called, calculate the lift, drag, and moment coeffi- 
cients (C'i, Cj, C m ) by calling the subroutines AEROT1 and AEROT2, 
respectively, which interpolate the aerodynamic table information. 
These ( C\i, and C m values have no time delay effects due to trail- 
ing edge or leading edge separation. These effects are accounted for 
in AERBED1 and AERBED'2 which call the subroutines SEPRATE1 
and SEPRATE2. 

4. SE PRATE 1 and SEPRATE2 

These subroutines calculate the lift and drag coefficients including the 
effects of trailing edge separation (TES) and leading edge separation 
(LES). The values from AERBED1 and AERBED2 before calls to 
SEP RATE 1 and SEPRATE2 assume no TES or LES. The effect of 
TES is to lag the lift behind the angle of attack. This lag occurs 
because the chordwise location of the TES point, does not instanta- 
neously follow the changes in surface pressure during motions of the 
airfoil. This effect is explained in detail in the mentioned References. 
The TES theory coded in SEPRATE1 and SEPRATE2 follow these 
references. However, in the original publications, the TES chordwise 
point is assumed to be known a priori ; this is not the case here. Be- 
fore the lagged TES point can be found, the unlagged TES point must 
be determined by solving the TES point equation in the references for 
the variable FNP. The resulting equation is: 



(2.55) 

(2.56) 


FNP is the value of the unlagged TES factor. Ci potential * s fh e P°~ 
tential normal force coefficient. C'/ )C 8i is the normal force coefficient 
from the airfoil table lookup performed in AERBED1 (or AERBED2). 
FNP is then lagged and applied to calculate the TES contribution, 
C'nnf, to the normal force, as described in the mentioned References. 

SEPRATE1 and SEPRATE2 then calculate the LES contribution, 
C nnv , to the normal force coefficient. This calculation is a simplified 
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version of the dynamic stall model developed in the mentioned Refer- 
ences. This simplified version accounts for one vortex release during a 
dynamic stall event. This is a limitation only for deep, sustained stall. 

The value of the normal force and chordwise force are then determined 
by the equations in the mentioned References. Since SEPRATE1 and 
SEPRATE2 need to return C\ and Cj , the chordwise and normal force 
coefficients are converted to these quantities before returning. It is 
assumed that the moment coefficient is appropriately represented by 
the value obtained from the interpolated table information from sub- 
routines AEROT1 and AEROT2. 

2,16.4 New Common Blocks 

The following subroutines contain one or more new subroutines for the indi- 
cia! aerodynamics: 1NITBED1, INITBED2, TRIM, MOTNR1, MOTNR2, 
AERBED1, AERBED2, TWA1, TWA2, INPTR1, INPTR2, SE PRATE 1, 
SEP RAT E2, INPTW1, and INPTW2. The new common blocks are as 
follows: BED1, BED2, BEDOLD1, BEDOLD2, BEDVEL1, BEDVEL2, 
BEDTWA1, BEDTWA2, BEDCON1, BEDCON2, BEDUSE1, BEDHSE2, 
BEDSEP1, BEDSEP2, DUCK1, and DUCK2. The list of variables in these 
new common blocks is extensive. For a listing of these variables, refer to the 
subroutines in the CAMRAD.Modl source code. 


2.16.5 Known Caveats 

1. The indicial aerodynamics model should not be used with the CFD 
interface since the TWA can not truncate the near wake as is required 
bv the CFD interface. 

2. Recently it has been discovered that there is a flaw in the TWA model 
in CAMRAD.Modl. The symptoms are that the loads do not decrease 
as they should near the blade tip, and the loading is not smoothed 
appreciably in a "‘‘near wake” fashion and can cause problems in the 
large core calculations of the rollup model (discussed later). The exact 
cause of the problem is not known at this time. Therefore use of the 
low resolution indicial aerodynamics option is not recommended at 
this time. 
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2.16.6 Extensions to High Resolution 

The indicia] aerodynamics, for high resolution is implemented as a sepa- 
rate post- processor program to account for the near wake given a far wake 
“forcing function”. The indicia! post- processor and its implementation will 
be discussed in a later chapter. 

2.17 Modifications for a Vortex Rollup Model 

2.17.1 Introduction 

In the original version of CAMRAD, the vortex wake consisted of several 
components. First, a near wake lattice model was used for the wake imme- 
diately behind the reference blade. The far wake for all blades (starting at 
the end of the near wake lattice model for the reference blade) consisted of 
a tip vortex and an inboard far wake “panel”, usually represented as one 
shed and one trailed vortex in the center of the panel. Normally, both were 
present with large vortex core radii. The geometry of the inboard wake was 
determined by a prescribed /rigid wake model. The tip vortex distorted ge- 
ometry was computed from either a prescribed /rigid wake model or from a 
free wake model. In calculating the influence coefficients for the tip vortex, 
it was assumed that the tip vortex geometry was determined solely bv the 
blade tip position and the tip vortex geometry model. The tip vortex rigid 
wake geometry model assumes that the distorted wake geometry depends 
only on the convection of the vortex endpoints by the freestream velocity. 
The free wake model assumes that the wake geometry is determined by the 
freestream convection of the vortex endpoints and by the vortex self-induced 
velocities. Also, the tip vortex strengt h was determined as a function of the 
maximum bound circulation on the blade, regardless of the distribution of 
bound circulation on the blade. In order to incorporate more physical prin- 
cipals into the wake structure/geometry analysis, a new vortex model was 
introduced — the “rollup model”. The rollup model modifies the wake ge- 
ometry calculations which are used in determining the influence coefficients. 
The rollup model does not directly affect the free wake calculations, but it 
modifies the results. 

Several parts of the rollup model will be discussed in the following sec- 
tions. But first, the method of implementation will be discussed. Then the 
rolled-up vortex positions will be discussed, followed by a discussion of the 
vortex position "phase-in” models. Following that will be a discussion of 
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the multi-core vortex model used for both tip and secondary vortices. And 
finally, the vortex “spin” model will be discussed. 

2.17.2 Method of Implementation 

The rollup model (Ref. [3]) was introduced as an additional loop in 
the Trim loop of GAMRAD.Modl. Originally, the Trim loop consisted of 
several wake stages: uniform inflow, rigid wake, and free wake. The rollup 
procedure was added as an additional step; the trim sequence now proceeds 
as follows: uniform inflow, rigid wake, free wake, free wake with rollup model 
(see Figure 2.29). At the end of the original three wake stages, the “free wake 
with rollup model” stage is entered. The first step in this stage is to calculate 
the rolled-up vortex positions for the tip and secondary vortices (this is part 
of the “ROLLUP Calcs”). These locations will be used to position the tip 
and secondary vortices relative to the original wake location determined by 
the free wake geometry calculations. Figure 2.30 shows the original positions 
of the wake endpoints and those that are shifted inward (in the rotor tip 
path plane) by an amount determined by the rollup calculations. T here 
is also a vertical (perpendicular to the tip path plane) shift that will be 
discussed later. The second step is to calculate the strengths for each vortex 
core in the multi-core vortex model. This step is also part of the “ROLLUP 
Calcs, in Figure 2.29. With the vortex multi-core strengths and vortex 
locations known, the wake geometry used in the computation of the influence 
coefficients is then modified by shifting the w r ake endpoints by a calculated 
amount based on the rollup model and by using the multi-core model for 
the vortex core structure (labeled “Rollup Stage” in Figure 2.29). If needed, 
the first several wake endpoints may be phased-in to their final location 
over a prescribed wake age interval. Once the wake influence coefficients are 
known, the controls-motion-circulation interations proceed as before. Next, 
new airloads are calculated in the trim loop using the new rollup vortex 
locations, and the process is repeated for a user specified number of rollup- 
trim iterations. 

2.17.3 Rolled-up Vortex Positions, Part 1 

The rolled-up vortex model for the inward shift of the vortices (“Part 1”) 
incorporates elements of the Betz inviscid rollup model (Ref. [17]) which was 
originally derived for trailed vorticity of fixed wing aircraft. Also included 
are adaptations from the works of Donaldson and Bilanin (Ref. [18]) and 
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Figure 2.29: The TRIM loop with the new Rollup Calculations included. 
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Bliss (Ref. [19]). In this model, the bound circulation distribution along 
the blade span is used to define and locate an axisymmetric tip vortex that 
has a circulation (strength) distribution which varies in the vortex radial 
direction (outward from the center of the vortex). The u rolled-up” position 
of the vortex (i.e., the final location of the vortex) is calculated in the new 
model to be the fully rolled-up vortex position. This is similar to the Betz 
method in that the final span wise location of the vortex far downstream 
is at the wing’s spanwise centroid of vorticity. However, the rollup model 
also defines a “secondary” vortex, inboard of the tip vortex, when certain 
criteria imply that there should be more than one vortex trailed into the 
wake. The present analysis is limited to a maximum of two vortices (a 
tip and a secondary). The inboard wake treatment is the same as original 
CAMRAD. The model addresses only the final rolled-up location of the 
vortex endpoint in planes nominally parallel to the rotor disk. The final 
vertical rolled-up location will be discussed later. From the Betz rollup 
model, it is seen that the tip vortex far downstream of a fixed wing should 
be placed at the spanwise centroid of vorticity; whereas, in CAMRAD. Modi, 
the tip vortex location, excluding the rollup effects, is determined by a free 
wake analysis. Due to the complication of the particular free wake geometry 
analysis in CAMRAD.Modl, a direct modification of the free wake geometry 
analysis was not attempted. Instead, a superposition approach was taken. 
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In this approach, the free wake geometry is determined without the rollup 
model, then the geometry is adjusted to account for the rollup effects. This 
adjustment takes the form of an inward radial shift in the wake endpoint 
locations according to the rollup location calculations. 

In the actual determination of the vortex rolled-up spanwise location, it 
is assumed that the circulation distribution is smooth radially. This leads 
to a smooth vorticity distribution. However, in typical helicopter BVI con- 
ditions, there can be significant spanwise circulation variations caused bv 
perpendicular BVIs. It is assumed in this modeling that these spanwise vari- 
ations caused by perpendicular BVIs do not contribute to sustained trailed 
vortices. This is, these spanwise loading variations contribute to local trailed 
vorticity, but do not create “long-lived” vortex filaments such as seen in a 
tip vortex. Since the modeling presented here is limited to two vortices (a 
tip and a secondary vortex), it is desirable to use a smooth bound circu- 
lation distribution to determine the final rolled-up positions of the tip and 
secondary vortices. Determining a smooth bound circulation distribution 
(/. ( eliminating effects of the perpendicular BVIs) is the role of the “large 
core calculation”. This large core calculation involves a “side’ calculation 
whereby the wake influence coefficients are found using a “large” vortex core 
radius (typically about 0.3 rotor radii). Once the large core wake influence 
coefficients are known, the resultant loads and thus circulation may be com- 
puted. The circulation from this “side” calculation is dubbed the “large core 
circulation” (also known as the “fat core circulation” in [3]). Figure 9 from 
reference [3] illustrates a comparison between a. full spanwise circulation 
distribution and the resultant large core spanwise circulation distribution 
(called the “fat core distribution” in Figure 9 of reference [3]). Here it is 
seen that the large core circulation has essentially removed the effects of 
perpendicular BVIs. 

In addition to the large core calculation, the effects of blade rotation on 
vortex stretching must be taken into account. Figure 2.31 shows a rotating- 
blade tip vortex with straight line segments being emitted at a radial location 
y t . To illustrate a concept, vortex line segments, signifying vorticity shed 
from locations y, are shown outboard and inboard of y t . Each segment 
length depends on the spanwise location for a given azimuth step size and 
on the rotational and inflow velocity at location y. Each filament segment 
has an induced field velocity given by 

srsi , 1P „ 

6v= (2.50 
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Figure 2.31: Rotor blade with line segment modeling of trailed vorticitv due 
to bound circulation distribution. 
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where r is the perpendicular distance to tlie vortex. In the rollup concept 
considered here, the vorticity filaments are drawn to and are entrained into 
the vortex at location y t , with the vortex length being defined at y t . When 
this is done, stretching (or compression) of the vorticity must occur in pro- 
portion to the ratio of its original segment length to its new segment length 
at y t . This serves, in principle to maintain the same fiv contribution at an 
observer located at a distance r in the wake. This concept is implemented 
by multiplying t he large core spanwise circulation distribution by a weight- 
ing factor. F. The weighting factor for the large core spanwise circulation 
distribution for the tip vortex is as follows: 


F = (1 — y) + /tsin xp ( 2 . 5 K) 

The weighted large (fat) core spanwise circulation distribution is also shown 
in Figure 9 of reference [3]. This weighted distribution is subsequently used 
to determine the locations y t and y s , the spanwise final rolled-up locations of 
the tip vortex and secondary vortex, respectively. r I he weighted large core 
spanwise circulation distribution, typically a smooth distribution, is used 
to locate the final rolled-up locations of the tip and secondary vortices. A 
spanwise vorticity centroid function can be found bv the following equation: 


mv ) = — ‘ — r ^„d„ 

yU> 7(y) -7(0)70 drt 
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(2.59) 

(2.60) 


where y is the blade spanwise coordinate pointing inboard starting at the 
tip, y(y) is the spanwise vorticity centroid as a function of //. and y(y) is the 
weighted large core circulation as a (unction of y. In practice, the derivative 
of 7 is determined by a forward difference scheme starting at the tip of the 
blade, and the integral is calculated using the trapezoidal rule. 

The current rollup model implements the above integral in three typical 
cases. "Case 1" (see Figure 2.32) assumes that the circulation distribution 
increases monotonically from a value of zero at the tip to a value of r +mal . 
at a radial station, y+ max (also labeled as point "A”). In this case, the 
minimum value of vorticity between the tip and y+ mar is actually at y+man 
and the value of the vorticity there is zero. For this case, only a tip vortex is 
needed. In actual implementation and for coding convenience, the secondary 
vortex exists along with the tip vortex for all cases. However, in Case 1, 
the secondary vortex location is identical to the tip vortex location and its 
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strength is zero. The fully rolled-up position of the tip vortex is defined 
as the centroid of vorticity between the tip and point A from the integral 
below: 


Vt = 


i f y=A fry(*/) 

y(A) Jo Or) 


qdi] 


(2.61) 


where y t is the vorticity centroid located between the tip and the point A. 
and 7 is the weighted large core circulation, as discussed above. 

If the function y{y) does not increase monotonically, it is assumed that 
two vortices are required; “Case 2 ” and “Case 3” deal with the two such 
possibilities. In Case 2 (see Figure 2.33), a non-zero minimum occurs in 
the vorticity distribution between the tip and y+ mar . For this case, y+ mar 
is labeled "B and the location “A is determined in the actual discretized 
radial solution as the radial station just prior to the station at which the 
function y(y) fails to increase monotonically. The final tip vortex rolled- 
up location here is determined, from Equation 2.61, to be the centroid of 
vorticity between the tip and point A. The secondary vortex location is found 
to be the centroid of vorticity between point A and point B as follows: 
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(2.62) 
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If the value of y s is ill-defined by this equation, the value ofy s is taken as a 
point half-way between points A and B. 

For "Case 3’ (see Figure 2.34), there exists a negative minimum in the 
circulation distribution, V- max , between the tip and the positive maximum 
bound circulation, r+ ma j-, which leads to an inflection point in the vorticity 
distribution. The maximum in circulation occurs at a radial station y +ma<r , 
labeled “B’\ And as before, point “A” (in the actual radially discretized 
solution) is located at the radial station just prior to the station at which 
the function y(y) fails to increase monotonically. Again, the tip vortex final 
rolled-up location is determined by the equation for y i above between the 
blade tip and point “A"; the secondary vortex final rolled-up location is 
determined by the equation above for y s between the points “A" and “B’\ 
In this case, the tip and secondary vortices will have the opposite circulation 
sense (/.e., the tip vortex will have a negative strength and the secondary 
will have a positive strength). 
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Figure 2.32: Final rolled-up location - Case 1. 
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2.17.4 Rolled-up Vortex Positions, Part 2 

The previous subsection discussed the final rolled-up location of the tip 
and secondary vortices in a plane nominally perpendicular to the rotor shaft. 
The final vertical rolled-up location in now considered. Previously in CAM- 
RAD. Modi, the tip vortex was assumed to originate at the tip of the blade. 
The original location of the tip vortex was defined by locating the blade tip 
in space, t hen setting the coordinates of the zeroth wake age to that coor- 
dinate. At subsequent wake ages, the wake endpoint location in space was 
determined by a combination of the blade tip position when the endpoint 
was “deposited” in the wake, the convection of the endpoint by the free 
stream, displacement of the endpoint bv a free wake distortion vector, and 
displacement of the wake endpoint by an amount prescribed by the rollup 
model. The previous subsection discussed the additional endpoint displace- 
ment term from the rollup in the horizontal plane. This subsection discusses 
the additional term for the vertical coordinate. For the rollup model, this 
vertical term is needed to account for the fact that certain vortices do not 
leave from the tip of the blade. 

For a. vortex (tip or secondary) that is inboard of the blade tip, the final 
vertical rolled-up position term is equal to the vertical location at the span- 
wise final rolled-up location on the blade. For example, of there were no 
coning or bending of the blade, the vortex final rolled-up location would be 
constructed from a spanwise location, and a vertical final rolled-up location 
equal to zero. This is because the tip and the spanwise rolled-up location 
are at the same vertical coordinate for this scenario. As another example, 
if a rigid, articulated blade were coned upward, the spanwise rolled-up lo- 
cation would be below the tip of the blade. Therefore, the wake endpoint 
would be shifted downward by the difference in these two coordinates (see 
Figure 2.35). With this additional term, the vortex appears toeminate from 
the blade at the location marked with an “X” in the figure, instead of the 
location marked with an “O”. 

2.17.5 Rolled-up Vortex Position Phase-in 

In the previous subsections, the final rolled-up locations of the tip and 
secondary vortices are discussed. Provisions are made to allow the vortex 
to migrate to or “phase-in” to the final rolled-up position from some initial 
spanwise location on the blade. Several models for this have also been 
developed and implemented in C'AMRAD.Modl. First, a simple functional 
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Figure 2.35: Z-coordinate for rolled-up vortices in rollup model. 


form was implemented to specify a phase-in scheme. Second, a model based 
on fixed-wing work of Spreiter and Sachs was attempted. This second model 
was abandoned after initial testing. In CAMRAD.Modl, the coding for this 
model still exists, but the variable controlling its usage has been hard coded 
such that the option cannot be used. Thus, the Spreiter/Sachs option will 
not be discussed in this documentation. 

I he prescribed phase-in model in ( AMR AD. Modi has several forms 
which can be applied to the tip and/or the secondary vortices. First, the 
variable ISECPII in namelist NLROLL determines if the secondary vortex 
will be involved in the phase-in process. Normally, the secondary vortex does 
not participate in the phase-in process as it is assumed that the secondary 
vortex will “form” and remain at the location y s . Therefore, the secondary 
vortex is normally considered to always be at its final rolled-up position. The 
tip vortex, on the other hand, is assumed to be “created” at some location 
on the blade outboard of the secondary vortex. The tip vortex is assumed 
to then phase-in to (/.e., migrate to) its final rolled-up position. Figure 
2.36 illustrates a wing with a “tip” vortex trailed. In the figure, note that 
at some downstream wake age, the vortex is located at its final rolled-up 
position. Conceptually, the phase-in function is smooth and continuous. In 
practice, the phase-in function is applied at given wake segment endpoints 
of fixed wake age; each vortex segment between these endpoints is a straight 
line. Also, note that the vortex originates from a user-specified fraction of 
the final rolled-up location relative to the tip of the blade. Between these 
two locations, a tenth-order polynomial of wake age, with user-specified 
coefficients, is used to phase-in the tip vortex. Normally, the initial location 
is assumed to be at the tip of the blade and the phase-in occurs linearly in 
age (by using only the first term in the phase-in polynomial). In a typical 
case, the phase-in is assumed to be complete (*>., the tip vortex is at its 
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Figure 2.36: Rolled-up location “phase-in” model. 


final rolled-up location) after one rotor revolution. 


2.17.6 Multi-Core Vortex Model (Tip Vortex) 

Another feature of the rollup model is the “multi-core vortex model”. 
From classical fluid mechanics, the velocity field induced by an infinite ideal 
vortex (a solution to Laplace’s equation) is as follows: 


r 

2 nr 


(2.63) 


where T is the circulation associated with the vortex, and r is the radial 
perpendicular distance to the point where the velocity is desired. This vortex 
is irrotational (except at the center), and has a velocity distribution that 
approaches infinity as the center of the vortex is approached. In “real” 
fluids, viscous and turbulent effects become dominant near the center of the 
vortex and thus prevent infinite velocities at the vortex center. Historically, 
many models have been introduced to account for the viscous effects near 
the center of the vortex. These models typically assume a particular velocity 
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distribution based on a vortex core size, r c . For example, a Rankine vortex 
model for the velocity is similar to that of an ideal vortex for 7 * > 7 ' C1 but for 
/■ < /'('■, a solid body rotation is assumed. This assumption leads to a linear 
velocity increase from zero at the center of the vortex to the value v = ~— 

Z7T7Y 

at the edge of the vortex core. Another example is a model dubbed here 
as the “Scully” vortex model and has an induced velocity distribution as 
follows: 


v 


2irr y r 2 -f r c 2 J 


(2.64) 


where T, ?*, and r c are defined as discussed above. The term in the paran- 
theses is known as the core factor and it serves to smoothly transition the 
induced velocity profile v from the ideal vortex to zero at the center of the 
vortex, with a peak induced velocity at the distance r c from the center of 
tiie vortex. The functional form of the core factor is arbitrary to a certain 
degree, but does provide a smooth transition in the velocity profile as a 
function of distance from the vortex center, unlike the Rankine vortex that 
has an abrupt change from a linear function to a 1 fr function at the edge 
of the vortex core (that is, at /* = r c ). 

Previously in CAMRAD.Modl, the tip vortex core was modeled (typi- 
cally) using the Scully vortex model above. The tip vortex core size was a 
single, constant input quantity. The present vortex core modeling for the 
tip and secondary vortices implements a multi-core model which ties the 
strength of the tip and secondary vortices to the large core circulation dis- 
tribution on the blade at each azimuth. Figure 2.37 illustrates the model 
of the structure of the far wake vortices which depend on the large core 
circulation distribution at each blade azimuth station. The vortices shown 
are taken as fully developed; the intermediate rolling-up process and aging 
process is not modeled here. The structure of the tip and secondary vortices 
are represented in Figure 2.37 as “sets” of concentric vortices of varying 
radius. 

In reality, the tip vortex strength would vary in a continuous manner 
radially outward from the center of the vortex. However, a discretized so- 
lution is desired in CAMRAD.Modl. Thus, the strength of the tip vortex 
will be calculated at a set number of user input core sizes. Up to ten in- 
put core sizes ((/y) p ) are possible {p is the index of the input core sizes; 
p = 1,2, 3, ..., 10). The only requirement on these core sizes is that they are 
spaced in such a manner that they adequately resolve the viscous core region. 
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Figure 2.37: Multi-core model (tip and secondary) shown at fully rolled-up 
locations. 

Using the discretized version of the spanwise vorticity centroid function as 
discussed earlier in this section, the spanwise vorticity centroid function be- 
comes, y(y;), where y, are the discretized radial locations. The discretized, 
weighted large core spanwise circulation, 7 (y s ), is assumed to be related to 
y(y t ) as shown in Figure 2.38. The circulations are linearly interpolated to 
the actual user input core sizes, (r c ) , where the subscript, p is the index 
of t he user input core sizes. This assumed relation holds for the tip vortex 
until the point labeled “A" is reached for any of the three possible “Cases” 
discussed previously. Any input core size that is outside of this interpolation 
range is assigned a strength of zero. The interpolated circulation values are 
represented by the symbol 7 ,,, where, again, p is the index of the user input 
core sizes. The st rengt h, 7 ' for the vortex of core radius r p is determined as 
the difference in circulation values between vortices of core size r p and r p _ { . 
For the tip vortex, this difference is as follows: 

h't)p = (7 t) P - (7t) P -i (-2-65) 

Note that if p = 1 the term (7f) P -i = 0. 

The most inner core radius (? 7 ) p for the tip ((r a ) p for the secondary) 
represents an approximate viscous core radius. This inner core size may be 
determined in two ways. First, it may be a constant, user specified core size 
that would be the smallest core size expected from the particular rotor being 
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Figure 2.38: Large core, weighted circulation as a function of y(yi) and 
input core sizes (r c ) p . 


examined. Or second, it may be determined by an empirical model. (The 
empirical model is discussed later as an option to the multi-core model.) 
The above definition of the multi-core model concept allows one to match 
nearly arbitrary large core circulation distributions to trailed vortex st rength 
distributions. It is seen that although some choices can be made in the code 
to define the inner core radii of the tip and secondary vortices (as options), 
the importance of the core radii as tuning parameters is significantly reduced 
when compared to the single core model. 

In order to maintain the same Sv contribution at an observer at r in the 
wake, the strengths of 7 ' are “unweighted” by the following factor: 

F = 71 =rr^ — (2.66) 

(1 -*/,) + // sin t 

The strengths of the all of the vortices in the multi-core representation of 
the tip vortex are now known, and the velocity field of this vortex can be 
calculated. The following formula is used to calculate the velocity field at 
any point due to influence of the multi-core tip vortex: 



(2.67) 


where p is the summation index, P is the number of vortex cores in the multi- 
core model, r is the perpendicular radial distance to the vortex center, (r c ) p 
is the p th user input core size, n is a user input integer used to vary the 
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vortex core model and the term in parantheses is the core factor. Note the 
similarity in the functional form of the velocity profile to the form of the 
Scully type vortex model discussed previously. Also, note that for different 
values of n many of the common core models are recovered. For example, if 
n = f . the Scully model is retained. If n — >■ oo, the Rankine vortex model is 
recovered. Typically a value of n = 2 is used. 


2.17.7 Multi-Core Vortex Model (Secondary Vortex) 

As discussed earlier, if a “Case 2” or “Case 3” situation is detected, a 
secondary vortex is assumed to exist at a final rolled-up spanwise location y s . 
Like the tip vortex, the strengths of the multiple concentric vortex cores that 
compose the secondary vortex are needed to calculate the velocity field due 
to the presence of the secondary vortex. Whereas the tip vortex multi-core 
strength distribution was related to the function y{y), the secondary multi- 
core distribution is related directly to differences in spanwise circulation on 
the blade. First, a spanwise origin, labeled point “C” in Figures 2.33 and 
2.34, is needed that lies between points A and B. The amount of vorticity 
outboard of point C and the amount of vorticity inboard of point C are 
assumed to be additive in the multi-core model for the secondary vortex. 
That is, the vorticity outboard of point. C will provide a vortex strength 
distribution in the multi-core model of the secondary, as will the vorticity 
inboard of point (’. These effects are the superimposed. Point C is located by 
finding the maximum rate of change of the weighted large core circulation, 
\(h, /dy\. In practice, the derivative is calculated using a forward difference 
scheme: 


— « 2i — liri (2.08) 

dy yi - //,- i 

where the subscript i starts at the location just inboard of point A and 
continues until point B. Also, in practice, the location of the maximum 
vorticity is determined by successively testing the calculated value of \dj/ dy\, 
as above, with 1.1 times the previously calculated value. If the radial station 
associated with the maximum slope is found in this manner, it is saved as 
point C'. If a maximum is ill-defined bv this process, point C is determined 
bv the radial station midway between the radial stations assigned to point 
A and to point B. Figure 2.39 illustrates points A, B, and C for a generic 
“Case 3” . The vertical axis is the relative circulation value at the given radial 
stations t The second (lower) horizontal axis represents the interpolation 


110 



of the relative circulation values onto distances equivalent to the user input 
secondary vortex core sizes (r c ) p outboard (between A and C) and inboard 
(between B and C) of the point labeled C. First, the circulation at each of 
these core sizes (/>) p in the outboard region is determined by the relative 
circulation distribution at that radial station. In the case shown (i.e.. Case 
d), these values will be negative since the circulation outboard of point C 
is less than that at point G. Next, the circulation inboard of point C is 
determined in a manner similar to that used in the outboard region. These 
circulations at equivalent vortex core radii (r c ) p are then subtracted from 
one another as follows: 

*)v li'c) Aboard ^ L ') p, outboard (2.69) 

where now, j p is the circulation value associated with the input user core 
size (r c ) p . The strength (vorticity) of each of these core sizes is now taken 
as the difference in the circulation between the cores as follows: 

( f s)p il$)p — (Ts)p — 1 (2. ( 0) 

where p is, as before, the user input vortex core size index. As before, if 
V = li the term 7 P _ i — 0 . Again, in order to “unweighC the strengths of the 
multi-core model, the (7')^ terms are multiplied by the following function 
F: 


0 “ Vs) +/isin ; 

Note that except for the term this F is just the inverse of the F used 
to weight the large core spanwise circulation distribution. Now the induced 
velocity field due to the secondary multi-core vortex is calculated in the 
same manner as the velocity field induced bv the tip multi-core vortex: 


P =P 
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(2.72) 


2.17.8 Multi-Core Vortex Model Options 

There are actually three multi-core model options in CAMRAD.Modl. 
The first, described previously, is the default option in CAMRAD.Modl 
when using the rollup model. This default option uses an array of constant 
core sizes for the multi-core tip vortex and the multi-core secondary vortex. 
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The second option, controlled by the variable 1C0RYCB in namelist 
NLROLL, is to use a single, empirically determined core size for the multi- 
core model. This option is referred to as the “single variable core, multi-core 
model”. This option uses the following empirically derived core sizes: 


(r c ), = 0.015 + 0.075y ( (2.73) 

(r c ) s = 0.015 + 0.075(y s - y A ) (2.74) 

where (r c ) t is the inner most tip vortex core size (p = 1), (r c ) s is the inner 
most secondary vortex core size (p = 1), y ( is the fully rolled-up location of 
the tip vortex as calculated above, y s is the fully rolled-up location of the 
secondary vortex as calculated above, y A is the radial location of point A, y 
is the distance from the tip of the blade to the inboard point y on the blade. 
With these core size definitions, the core sizes for the tip and secondary 
vortices vary azimuthally. The strengths of the core for each the tip vortex 
and secondary vortex are determined as before with the index p = 1. 

The third option, refered to as the “variable multi-core, multi-core 
model”, is a hybrid of the two previous models. First, the user defines a set 
of constant radius vortex core sizes as done in the default option. Then the 
“variable multi-core, multi-core model” will first, calculate the distribution 
of strength in the constant core sizes as in the default option. With the 
strength distribution known, the empirical core sizes discussed previously 
are assumed to be the minimum viscous core radius. With this assumption, 
the “closest” core size in the constant core size array is reset to the mini- 
mum core size. Any strength inside of the minimum core size is summed 
and assigned to the newly defined minimum core size. Since the minimum 
core size will be between two core sizes in the constant core size array, an 
interpolated amount of strength is also removed from the core size that is 
just larger that the minimum core size. The amount that is added to the 
minimum core size is then removed from the larger core. 4 he purpose of 
this last step is facilitate a smooth transition between two core sizes as the 
minimum core size increases and decreases azimuthally. For example, if the 
minimum core size is only slightly smaller that a particular core size in the 
constant core size array, the minimum core will receive all strength inside 
its radius and most of the strength from the core size to which it is closest. 
Subsequently, the core size that is closest to the minimum core size will have 
almost zero strength. This is illustrated in Figure 2.40. The user input con- 
stant core sizes in this Figure are rl, r2,and, r3. The empirically determined 
core size is (»>),. In practice, the core size r2 is set to the value (/>),. This 
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Figure 2. .19: Spauwise locations of A, B, and ( in multi-core model for 
secondary vortex. 

core (r c ) f will have the summed strength of (1) rl and (2) r‘2 and (3) an 
interpolated strength between cores r'2 and r3. The strengths of cores rl 
and r'2 are subsequently set to zero for the remainder of the calculations; 
the strength of core r3 is reduced by the same amount that was gained by 
(r c ) t in the interpolation step above (step (3)). 

2.17.9 Multi-Core Model Caveats 

The wake influence coefficient calculations in the original version of CAM- 
RAI) were developed for a constant vortex core size. In the multi-core model 
options where a core size varies azimuthally, these same influence coefficient 
calculation precedures are still used. That, is, the effect of a. linearly varying 
vortex core size on the value of influence coefficient for each wake segment 
endpoint is not taken into accounted. 

2.17.10 Vortex Pair Spin Model 

The present use of the Scully free wake model in C'AMRAD.Modl puts 
constraints on the way that the wake geometry can be modified to account 
for effects of the rollup model. At present, any particular wake endpoint of 
a tip and secondary vortex is associated with a wake endpoint from the free 
wake geometry calculation. This free wake geometry endpoint is adjusted 
after the free wake geomtry calculation to include the rollup model vortex 
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Figure 2.40: Variable multi-core, multi-core model interpolation schematic. 

position calculations (discussed previously) in the influence coefficients cal- 
culations. Note that the free wake geometry calculation does not account for 
the fact that in some parts of the rotor wake, there exist, both a tip and sec- 
ondary vortex. Figure 2.41 depicts a portion of a wake filament w hich splits 
into a tip and secondary vortex over part of its length. The two vortices in 
this case are shown with opposite (but not necessarily equal) strengths, I s 
and I\. A model was implemented to account for the mutual influence of 
the tip and secondary vortices on each other. This model, called the “spin 
model, is used to rotate the wake endpoints to a new position based on a 
two-dimensional model of the influence of the tw'o isolated vortices on each 
other. For example, in Figure 2.41, both vortices wi 1 1 tend to rotate upwaid 
due to the others' influence. If the secondary vortex were stronger than the 
tip vortex, the tip vortex would be moved higher that the secondary. Figure 
2.42 shows a two dimensional scenario used to calculate effects of the spin 
model. The following geometric relations can be derived using the figure: 

P t = Pt(r t ,z t ) (2-75) 

P s = P s (r s ,z s ) (2.76) 
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(2.77) 



Q = 


n 

cos 0 Z 



(2.78) 

(2.79) 


AC =10 -Cal (2.80) 

where P \ and P s . are the vectors to the tip and secondary vortex endpoints at 
a given azimuth location. At a given azimuth location, the velocity influence 
of each vortex endpoint on the other is calculated using a two dimensional 
infinite multi-core vortex model: 


P =P 

= E 


(7 .•?)/> AC 


1 (AC) 2 + ((^)p)seco7irfarj/ 


(2.81) 


p=p 

v- = E 

v - 1 


+ {{ r c)p)tip 


where \ t is the velocity influence at the tip vortex due to the secondary 
multi-core vortex, and where \' s is the velocity influence at the secondary 
vortex due to the tip multi-core vortex. This model does not include the 
effect of the exponent n that was discussed in a previous section (/.e., a 
"Scully type vortex core model is still assumed here). Using the above 
quantities, new quantities can be calculated: 


t sCt — l/C 
v a - Vt 
VjSl 
C t ~~ C 0 

Vt - y s 

Pt-Ps 

J t *r.s 

p, - p,| 

n| 

=0 = (G - 0) sin 0 o 




(2.83) 

(2.84) 

(2.85) 

(2.86) 

(2.87) 

( 2 . 88 ) 
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r 0 = (C> - C<) cos ^o 


(2-«9) 

where is the rotor rotational rate, x, y,z with subscripts t,s are the co- 
ordinates of the wake endpoints being examined for the tip and secondary 
vortices, respectively. It is necessary to identify three “Cases’" in order to 
correctly apply the spin model to a particular situation. (Note that these 
cases are not associated with the three Cases discussed in earlier subsec- 
tions.) Once the spin Case is identified, several parameters can be evaluated. 
If first statement is true, then the following two quantities are set: 


Cast 1 : C < Cs < 0 

(2.90) 

0 t ,„ = 0. 

(2.91) 

0 t = uJT c T 

(2.92) 

0s,o = 0. 

(2.93) 

0 S = u)T c T 

(2.94) 

Cast 2 : £ s < Q < Q, 

(2.95) 

01.0 = X 

(2.96) 

0 t = U )T C T + 7T 

(2.97) 

0s.o — 

(2.98) 

0 g = UT C T + 7T 

(2.99) 

Cast ;{ : C < Co < Ct 

(2.100) 

0,.o = 0. 

(2.101) 

0, = u>T c T 

(2.102) 

0„,o = * 

(2.103) 

01 = u !T C T + n 

(2.104) 


where t c is a user defined multiplier for the rate of spin. Another user input 
constant, r„ defines the wake age at which the spin calculations begin. For 
example, of r 0 = 0.0. the spin calculations begin immediately. Once the spin 
Case has been identified, the following relations are applied to modify the 
tip and secondary vortex endpoint locations to account for the spin “effects” 
discussed above: 
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Figure 2.41: Vortex "spin” model schematic. 


/ cos (00 + 0 t ) - cos( 0 „ -I- 6 Uo ) cos V' \ 

P' = Pf + |0 ~ Col co8(0,j + 0? ) - cos(6>„ + 0, „) sin iji I (2.105) 

\ «i»(0o + 0/) - sin (0 O + *,,„) ) 


( cos(0 o + 0 S ) - cos(0 o + 0 V ,) cos \ 

^ = Pa + 10 - Col cos(0 o + 0 S ) - cos (0 O + 0 S:O ) sill V (2.106) 

\ si"(0, + 0) - sin(0 6 + 0,,o) / 

where P' t and P' s are the new coordinates of the vortex endpoints included in 
the spin model, P t and P s are the new coordinates of the vortex endpoints 
before the spin model. 

2.18 Namelist Reading Subroutine Changes 

2.18.1 Introduction 

Many new variables and capabilities have been added to CAM- 
RAD. Modi. The original version of CAMRAD included a BLOCKDATA 
capability for input of many code parameters that were not normally al- 
tered in the course of a prediction task. Variables that are changed rou- 
tinely were modified using namelist inputs. The newly added variables in 
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Figure 2.42: Vortex locations for “spin” model from previous figure. 

C'AMRAD.Mod 1 are normally given values using the namelist inputs using 
new namelists. In some cases, existing namelist reading subroutines were 
modified to input the new namelists; in other cases, new namelist reading 
subroutines were added. This section describes each of the changes and ad- 
ditions. Some of the changes are discussed in other sections, but this section 
is intended to provided a concise listing of the changes. 


2.18.2 Subroutine INPTN 

This subroutine controls the reading of all namelists and reads namelist 
NLTRIM. Two new variables, FAC’TM and OPMXFWG, were added to 
NLTRIM; default values for these two variables were also added. Two 
new namelist reading subroutine calls were added to this subroutine. The 
new subroutine INPTCFD was added immediately following the call to 
INPTVVl. INPTCFD reads the namelist NLCFD which contains parameters 
controlling the CFD interface. Immediately following the call to subroutine 
INPTCFD, a new call to subroutine INPTM1 was added. INPTM1 reads 
the namelist NLMEAS which controls the usage of measured blade motion 
and measured C„ information in C'AMRAD.Mod 1. 
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2.18.3 Subroutines INPTRl and INPTR2 

INPTRl is one of the most, highly modified of the namelist reading sub- 
routines. Originally this subroutine read only namelist NLRTR. In addi- 
tion. it now reads the following namelists (1) NLHHC , (2) NLHH02, (3) 
NLHIRES, (4) NLBED, and (5) NLSWP. NLHHC and NLHHC2 rontol 
the usage of Higher Harmonic Control; NLHIRES controls the usage of the 
HIRES section of CAMRAD.Modl; NLBED controls the usage of the In- 
dicial Aerodynamics in the low resolution section of CAMRAD.Modl; and 
NLSWP controls the usage of the aerodynamic sweep correction in the high 
and low resolution sections of CAMRAD.Modl. All of these namelists must 
exist in the script file after NLRTR, even if they are empty. Many of the 
variables in these namelists are given default values in INPTRl. All of the 
above changes apply to subroutine INPTR2. 

2.18.4 Subroutine INPTWl and INPTW2 

INPTW1 originally read only the wake namelist NLWAKE. Now, it also 
reads the namelists NLBURST and NLROLL. NLBURST controls the pa- 
rameters in the vortex bursting model. NLROLL controls the usage of the 
vortex rollup model. INPTWl now also performs many checks to be sure 
that certain incompatible variable combinations are not attempted by the 
user. Some non-functioning options are also flagged in this subroutine with 
error messages printed to the standard output file. All of the above changes 
apply to subroutine INPTW2. 

2.19 Fuselage Aerodynamic Tables 
2.19.1 Introduction 

CAMRAD.Modl was modified to include changes to the fuselage aerody- 
namics made by Sikorsky under contract to NASA Langley. These changes, 
originally developed by Sikorsky and subsequently implemented in CAM- 
RAD.Modl, modify the fuselage aerodynamics in CAMRAD.Modl by re- 
placing the empirical analysis by “look-up” tables. With this model, the 
empirical aerodynamic formulae for the wing, the body, the horizontal tail, 
and the vertical tail, are replaced by a table “look-up’ of aerodynamic char- 
acteristics for a new "wing-body-tail . These aerodynamic characteristics 
for the wing- body-tail are applied at a wing-body-tail location instead of 
at separate locations for the wing, the body, the horizontal tail, and the 
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vertical tail. See reference [20] for details. Several changes, mainly common 
block size corrections, were made to the original Sikorsky implementation. 

2.19.2 Code Changes 

1. The new subroutine BODYAT from Sikorsky was added. 

2. The variable WBTTAB was added to the common block BADATA 
in subroutines BODYA, BODYAT, BODYF, FILEI, 1NITB, 1NPTN, 
PERF, PRNTB. ROTNET, and TRIMP. 

3. The variable WBTTAB was initialized to zero in subroutine INPTB. 

4. A subroutine call to BODYAT was added to subroutine BODYF if 
WBTTAB is set to one. 

2.19.3 Extensions to High Resolution 

This modification has no bearing on the high resolution calculations. 


2.20 Machine Dependencies 

2.20.1 Time 

Original modifications to CAMRAD.Modl were made on a VAX with a 
VMS operating system. The code was ported to a DEC ALPHA workstation 
operating with a UNIX (OSF/1) operating system. Changes were made to 
the code to allow proper operation of the time and date stamps used in the 
code output. In the executive program, CAMRAD, in the input preparation 
program. BLOCKFILE, and in the airfoil preparation program, AIRFOIL, 
all calls to the system subroutine (TIME were changed to the following: 
CALL 1TIME (IHOUR, IMIN, ISEC) 

and the information is encoded into a variable, such as JDTIME, using the 
following: 

INTEGER JDTIME (2) 

ENCODE (8,901 .JDTIME) IHOUR, IMIN, ISEC 
901 FORMAT (12, 1H:, 12, 1H: , 12) 

In the subroutine TIMER, all calls to ITIME were changed to the following: 
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INTEGER I ARRAY (3) 

CALL ITIME (I ARRAY) 

and the requested program time in seconds is calculated using: 

T = IARRAY (3) *3600 . + IARRAY (2) *60 . + IARRAY(l) 

2.20.2 Date 

In the executive program, CAMRAD, the input preparation program. 
BLOCKFILE, and the airfoil preparation program, AIRFOIL, all calls to 
the system subroutine IDATE were changed to the following: 

CALL IDATE (IMONTH, IDAY, I YEAR) 

and the information is encoded into a. variable, such a,s JDDATE, using the 
following: 

INTEGER JDDATE (2) 

ENCODE (8, 902 .JDDATE) IMONTH, IDAY, I YEAR 
902 FORMAT (12, 1H/, 12, 1H/, 12) 

2.20.3 Dimension Statements 

Some FORTRAN compilers do not allow a dimension statement in a 
subroutine to be similar to the following examples: 

DIMENSION A(l) 

DIMENSION AA(NM, 1) 

Therefore, in order to make the code more portable, the dimension state- 
ments of the form of the prior examples were changed to the following: 

DIMENSION A(*) 

DIMENSION AA(NM,NM) 

when the intent of the dimension “1” is to pass an arbitrary length array to 
the subroutine. This replacement with the character tells the program 
to pass the array with the same dimensions as it has in the calling routine. 
In the second example, the actual dimension, NM, is used instead of the 
“ 1 ”. 
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2.20.4 Debug and Input Data prints 

Originally the “INPUT DATA" section printed to standard output de- 
pended on the variable LEVEL. For example, if LEVEL = 0, the wake input 
quantities were printed for uniform inflow only. To get the input data printed 
for the free wake portion, one was required to execute the program through 
the free wake stage. This dependency was eliminated by removing several 
IF statements in the subroutines PRNTW1 and PRNTW2. (Actually, they 
were “commented out", not deleted.) Now, all the input parameters for the 
wake are printed, regardless of the value of LEVEL. 

The COMPLEX variable KEPSI was removed from the debugging 
namelist DBINC in subroutine INITC. This was done because the DEC AL- 
PHA did not allow a COMPLEX data type to be printed using a namelist 
with other data types present . This removal impacts only the printing of the 
variable KEPSI in the initialization stage of the program if the debugging 
flag DEBUG (3) = 2. 

2.20.5 File Handling 

The subroutine FILEV was altered to operate in a UNIX environment. 
The original OPEN statements that have VAX-specific keywords, such as 
READONLY, were altered to remove non-standard keywords. In addi- 
tion. the original logical file names, determined from VMS system calls, 
are now obtained from the UNIX operating system by the system call 
“getenv (environmentvariable, actualfilename)" . The environment variables 
are as follows: INPUTFILE, AFTABLE1, AFTABLE2, RESTARTFILE, 
E1CENFILE, LPOUTDB, LPOUTPP, LPOUTPUT, LPOUTLIN, NL1N- 
PUT. AFTABLE, AFDEC’Kl, .... AFDECK10. 

Other files are opened without machine dependent elements in the OPEN 
statements. Filenames in these OPENs are hard coded, such as “AL- 
PHA P.l) AT". Other unit numbers are used in various places and each of 
these unit numbers must be linked to a file via the “In -s command. All 
unit numbers used are listed in the User s Manual. 

2.20.6 Logicals and DATA statements 

In subroutines that create printer-plots, DATA statements were originally 
used to set values of variables that are declared as a LOGICAL*! type. This 
is not allowed by all compilers. To make the code more portable, these dec- 
larations were modified. INTEGER variables were created corresponding to 
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the affected LOGICAL variables. These INTEGER variables have the same 
name as the LOGICAL variables with an “IP prepended. Then the INTE- 
GER variables are initialized using the DATA statement and the LOGICAL 
variables are EQUIVALENCEd to these INTEGER variables. For example 
the declarations of the type: 

LOGICAL* 1 THING 
DATA THING /1H-/ 

were changed to: 

INTEGER IITHING 
DATA IITHING /1H-/ 

EQUIVALENCE (IITHING, THING) 

This type of change was made in the following subroutines: BODEPP. 
FLUTM, GEOMP1, GEOMP2, HISTPP, NOISR1, NOISR2, POLRPP, 
S I ABM. TRCKPP, TRIMP, and in the airfoil preparation program. 

2.21 Miscellaneous Changes and Bug Corrections 

1. In the subroutines FILEV, an error was corrected. The variable BKL- 
DAT was corrected to BLKDAT. This error effected only the reading 
of a BLOCK DATA program as an ASCII input file. Normally, the 
input data is read as a binary file. 

2. In the subroutines FLFTRl and FLUTR2, a bug was corrected. The 
variable OP I ION was corrected to OPFLOW. This bug had an effect 
only in the flutter analysis. 

3. 4 wo subroutines, PRNTHRl and PRNTHR2, were added to print 
information contained in namelists NLHHC, NLHHC2, NLHIRES, 
NLBED, NLSVVP, NLBURST, NLROLL, and NLCFD to standard 
output if the parameter NPRNTT = 1 in namelist NLTRIM. 

4. The common block WORK was defined originally in two different 
places in the code. The WORK common blocks in GEOMP1 and 
GEOMP2 were renamed WORKGl and WORKG2 to eliminate the 
bug. These common blocks were not used outside of GEOMP1 and or 
GEOMP2. 
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5. Where practical, many of the exponent operators, ** , from the most 
inner loops in CAMRAD.Modl have been changed to multiplications 
since in some instances, multiplication is many times faster than ex- 
ponentiation. 

b. Two subroutines. INITHR1 and 1NITHR2. were added to subroutine 
1NIT to initialize radial HIRES parameters. Addition of these subrou- 
tines reduces the number of input parameters required to run the high 
resolution part of CAMRAD.Modl. 


124 



Chapter 3 

HIRES 

3.1 Introduction and Solution Procedure 

3.1.1 Introduction 

For many rotorcraft analysis tasks, a low temporal (azimuthal) and spa- 
tial (spanwise) resolution analysis of the aerodynamics and dynamics of the 
rotorcraft is sufficient. For example, in a performance analysis of an entire 
vehicle, 15 degree azimuthal resolution and 15 radial stations on the blade 
span may be sufficient to determine the average forces generated by the 
rotor. These average forces may be used to determine a steady state config- 
uration for the rotorcraft. There are several comprehensive rotorcraft tools 
available to analyze these configurations and flight conditions. One such 
example of a comprehensive rotorcraft code is the original version of CAM- 
RAD. Written using a low resolution azimuthal and spanwise discretization, 
limits were originally imposed such that azimuth step sizes of 15 and larger 
were required to be used (10 and larger if not using the free wake analysis). 
Also, the maximum number of spanwise locations on the reference blade 
was thirty. But, a high radial resolution and high azimuthal resolution are 
needed to calculate high resolution airloads for use in prediction of loading 
noise and Blade- Vortex Interaction (BVI) noise. 

Experience has shown (Ref. [8]) that an azimuthal resolution of 1 , 
and a radial discretization of 75 radial stations along the span, are often 
adequate to resolve the BVI unsteady aerodynamic loading required for 
the acoustic analysis. One possible (conceptual) method to increase the 
resolution of a low resolution code is to simply redimension all relevant 
arrays in the comprehensive code and execute the code at that increased 
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resolution. One factor limiting the success of such an approach is the large 
increase in computation time (CPU time) that would be required to trim 
the rotor. Some of the increase in computation expense would be the fact 
that, internal to the code, many revolutions of the rotor(s) are required to 
converge the internal motion and circulation loops of the code. Increased 
resolution in these loops would slow the convergence of the comprehensive 
code. In addition, free vortex wake models can become unstable as the 
vortex length decreases, due to a higher resolution trim solution. However, 
even if this approach were practical, the general complexity of the free wake 
algorithm implemented in CAMRAD would make such global free wake code 
changes prohibitive. 

Another approach is to allow the rotorcraft. t.o achieve a trimmed state at 
a low resolution, effectively not changing the comprehensive code, then ap- 
ply a post-trim analysis to “reconstruct” a higher resolution solution from 
the low resolution solution. A true post- processor type of analysis would 
apply a “stand-alone” code to the output information of the low resolu- 
tion performance code. For example, a separate post-processor code might 
interpolate all output information up to a higher temporal and spacial reso- 
lution, or might execute a totally different analysis, such as a CFD analysis, 
to operate on the low resolution information. In this chapter, a modified 
post -processor type of analysis is used for the high resolution portion of the 
code. Since the high resolution portion of the code, as implemented, is not 
a completely stand-alone code, this approach may be thought of as a hybrid 
post- processor an alysis. 

With a hybrid post-processor analysis, there are a number of paths avail- 
able to determine the high resolution loading solution. One such path, using 
an external CFD code, was discussed in the section on the CFD interface 
in Chapter 2. This current chapter will focus on another method that has 
become known as HIRES. This method includes all of the coding invoked 
at the end of the trim loop in CAMRAD. Modi (excluding the CFD inter- 
face, the ROTO.NET/ WOP WOP interfaces, the Flutter analysis, and the 
Transient analysis). 

3.1.2 Solution Procedure 

As discussed previously, a calculation of high resolution wake induced 
velocities, and thus loading, would not be feasible in the trim loop of CAM- 
RAD. Modi. Thus, a modified post-processor type of analysis (known as 
HIRES) was added to the end of the trim loop to process the information 
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into a high resolution solution. This solution procedure effectively uses the 
low resolution wake solution and kinematically reconstructs the wake at az- 
imuthal locations between known low resolution locations. HIRES also uses 
the modal analysis results from CAMRAD.Modl to define the blade position 
at a given high resolution azimuth location. With this information known, 
interpolations are use to calculate the high resolution information. The 
interpolation technique conceptually involves (1) linearly interpolating the 
wake geometry in wake age and azimuth, such that the tip vortices kinemat- 
ically translate from one known low resolution location to the next known 
low resolution location, (2) linearly interpolating the blade shape radially 
between known low resolution collocation points on the blade, (3) linearly 
interpolating the blade bound circulation to a higher azimuthal and radial 
resolution. With these quantities known, the wake influence coefficients 
may be determined at a high azimuthal and radial resolution. Then, with 
the influence coefficients known, the wake induced velocities, and thus the 
loading, may be calculated at a high azimuthal and radial resolution. The 
actual details of the HIRES implementation can differ from the above con- 
ceptual implementation. These differences are due to practical issues such 
as code execution time constraints, coding convenience, etc. The details of 
the solution procedure are presented here. 

One task involved in obtaining high resolution airloads is to calculate 
the wake induced velocities at a high resolution. This calculation requires 
knowledge of the blade position, the bound circulation, and the wake ge- 
ometry at a high azimuthal (temporal) and radial (spatial) resolution. The 
azimuthal blade position may be determined at any azimuthal location with- 
out modifications to the code. This is true because CAMRAD.Modl uses 
a modal analysis and stores the harmonics of the generalized coordinates of 
the blade mode shapes at user specified low resolution (spanwise) collocation 
points on the reference blade. Storing these quantities in this manner allows 
one to compute the coordinates of any low resolution (spanwise) collocation 
point at any azimuth station. A high spanwise resolution is achieved by 
linearly interpolating coordinates amongst the low resolution radial colloca- 
tion points. In order for this interpolation technique to closely approximate 
the true blade shape, it is required that as many low resolution spanwise 
collocations points be used as possible. 

The high resolution bound circulation is obtained initially by linearly 
interpolating the low resolution bound circulation to a high azimuthal and 
radial resolution. Since the bound circulation is used in determining the 
wake induced velocity, the interpolated values will be used to initiate the 



high resolution solution. 

In determining the wake geometry at a high resolution, the locations 
of the tip (and possibly secondary) vortices and the locations of the in- 
board wake elements are needed. For the tip (and secondary) vortices, the 
high resolution vortex end points are located by kinematically translating 
from one low resolution location to another. Linear interpolation is used 
to locate these intermediate endpoints (see later section entitled “Vortex 
Segment Location”). In addition, the inboard wake elements are located by- 
found by linear interpolation between low resolution endpoint locations, (see 
later section entitled “Vortex Segment Location”). With the high resolution 
blade coordinates known, the high resolution bound circulation known, and 
the high resolution wake geometry known, the wake induced velocities and 
aerodynamic loading may be computed at a high resolution. 

3.1.3 Implementation of Solution Procedure 

Historically, the first high resolution modification was the development 
of a far wake model. This model included only the effect of the tip vortex 
starting at the blade tip and existing for the same number of wake spirals 
as in the low resolution portion of the code. As such, the entire inboard 
vortex wake and lattice near wake were not modeled. This model was used 
in Reference [8]. At the next HIRES development stage, the inboard far 
wake elements were modeled. Several new models were also implemented in 
the far wake solution including a vortex segmentation model to effectively 
smooth the straight 10 vortex segments automatically as needed. 

In the subsequent HIRES development stage, two methods were devel- 
oped simultaneously to account for the near wake behind the reference blade. 
One method was implemented as an option inside the high resolution por- 
tion of CAMRAD.Modl and the other was developed as an independent 
post- processing code. The near wake model inside the high resolution por- 
tion of the code is a vortex lattice model roughly similar in nature to the 
vortex lattice model in the low resolution portion of the code. This method 
will be discussed in a later section of this chapter. The other method, which 
is an independent post-processor, called the Indicia! Post-Processor (IPP), 
implements a time domain indicia! aerodynamics formulation derived from 
the methods of T. S. Beddoes and Gordon Leishman (Ref. [13, 14, 15, 16]). 
This method will be discussed in Chapter 4. 

As discussed earlier, Figure 1.1 shows an outline of the CAMRAD.Modl 
code with HIRES represented as a rectangle. The arrow into the HIRES 
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Figure 3.1: CAMRAD. Mod 1/HIRES Schematic 


rectangle conceptually shows that HIRES is executed after the Trim loop of 
CAMRAD. Modi. Figure 3.1 expands the HIRES rectangle of Figure 1.1, 
to illustrate the procedures followed in HIRES. First, the rotor is trimmed 
using the low resolution portion of the code (CAMRAD.Modl), then the low 
resolution circulation is linearly interpolated azimuthally and radially, to a 
high resolution. The wake influence coefficients are then determined using 
only the “far wake". The far wake consists of 10 agewise vortex segments 
and an inboard (rigid) far wake model (to be discussed later). With the 
far wake influence coefficients determined, the wake induced velocities and 
airloads (associated with the far wake only) may be calculated. If determined 
to be necessary (by the user), this “Far Wake Loop” (see Figure 3.1) may 
be repeated. Normally, only once through the Far Wake Loop is required. 
At the end of this iteration loop, the user has the option to use the high 
resolution vortex lattice near wake or the IPP as discussed earlier. 

If the user chooses the vortex lattice near wake model (the option is 
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actually chosen at the beginning of the code execution by setting appropriate 
parameters in namelist NLHIRES), the execution of HIRES automatically 
will continue from the Far Wake Loop into the "HIRES Lattice Near Wake 
box shown in Figure 3.1. Otherwise, the code will stop execution at the end 
of the Far W'ake Loop. This is the path that is chosen (that is, stopping 
after the Far Wake Loop) if the 1PP is to be used for the high resolution near 
wake calculations. Note that if a second rotor is also included, the entire 
high resolution procedure may be repeated for rotor-2. This is accomplished 
by replacing the “1” in each box of Figure 3.1 by a “2”. 


3.2 Initialization 

There are a number of variables used in HIRES that are linearly interpolated 
from the low resolution solution of CAMRAD.Modl to a high resolution 
(azimuth ally and/or radially) and do not change during the HIRES phase 
of C AMRAD.Modl. These interpolated variables are calculated before the 
Far Wake Loop to avoid unnecessary interpolations during the computation- 
ally intense portions of HIRES. This calculation stage before the Far W r a,ke 
Loop is called “Initialization”. (Note that some initialization of HIRES pa- 
rameters takes place at the same time the low resolution parameters are 
being initialized. These initialized parameters are ones that depend on the 
values of user input low resolution parameters such as blade chord, twist, 
( tc. These high resolution parameters are found by radially interpolating 
(linearly) the low resolution parameters to the high resolution radial stations 
provided by the user in namelist NLHIRES (variable RAEINT). This ini- 
tialization of the radial parameters is accomplished in subroutines INITHR1 
and INITHR2 which are called just after the call to subroutine INITR2 in 
subroutine IN IT.) 

As seen in Figure 3.1, there are four subroutines before the far wake 
loop is begun. The four subroutines handle the initialization of the high 
resolution parameters ( e.g bound circulation, vortex burst parameters, 
etc.) that depend on the low resolution solution. The functions of the 
initialization subroutines are described below. 

Subroutines INITRAD1 and INITRAD2 (for rotor-1 and rotor-2, respec- 
tively) linearly interpolate, azimuthally and radially, the additional velocity 
distribution over the rotor disk from the tunnel/fuselage correction method 
to a high resolution for use in the high resolution induced velocity calcula- 
tion. 



Subroutines INITCRCl and INITCRC2 (for rotor- 1 and rotor-2, res pec- 
lively) initialize the high resolution bound circulation and maximum bound 
circulation, azimuthally and radially, by linear interpolation. For the far 
wake calculations, only the maximum bound circulation is needed at a high 
resolution. The high resolution bound circulation distribution will only bo 
used if the internal high resolution vortex lattice near wake model is used. 

Subroutines INITBST1 and INITBST2 (for rotor-1 and rotor-2, respec- 
tively) initialize the high resolution tip vortex burst parameters from the 
low resolution solution. 

Subroutines ROLLUPH1 and R0LLUPH2 (for rotor- 1 and rotor-2, re- 
spectively) calculate parameters needed for the high resolution tip vortex 
roll u p calculations. These parameters are linearly interpolated from the low 
resolution rollup calculations to high resolution values. 


3.3 High Resolution Far Wake 

3.3.1 Introduction 

After the initialization of the HIRES parameters, the far wake loop is 
executed. A subroutine, CIRCCAL1 (or CIRCCAL2 for rotor-2), is called 
to lag the circulation in the far wake loop (if the far wake loop is to be 
executed several times, that is). This lag factor is used similar to all of 
the other iteration lag factors in CAMRAD.Modl. However, experience has 
shown t hat normally one Far Wake Loop iteration is sufficient and the lag 
factor is simply set to unity. 

After the subroutine CIRGGAL1 is called, the far wake influence co- 
efficients are calculated at a high resolution by the subroutine WKG1INT 
(WKG2INT for rotor-2). This subroutine comprises the bulk of the far wake 
high resolution calculations. Once the influence coefficients are calculated, 
the wake induced velocities are calculated using the high resolution circu- 
lation distribution. These velocities are then combined with velocities due 
to blade rotation, blade motion, etc. at each spanwise section. Once these 
total far wake velocities are known, the angle of attack and Mach number 
may be computed. These angles of attack and Mach numbers are then used 
in an airfoil table to determine the airloads due to far wake effects. 
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3*3.2 Far Wake Influence Coefficients 


The high resolution far wake influence coefficients are used to calculate 
the normalized induced velocity due to each vortex element in the wake 
system at a particular point on the rotor blade. For the high resolution far 
wake influence coefficients calculation, the far wake consists ol the bound 
vortex of all blades except the reference blade, the inboard far wake of each 
blade, and a far wake tip vortex from each blade. The bound vortex for each 
blade, excluding the reference blade, consists of a straight line vortex from 
root to tip with a vortex core size of 25% of the mean blade chord. This 
model follows the low resolution model of the bound vortices. The inboard 
far wake consists of a two vortex segments, one in the trailed direction and 
one in the shed direction, on each vortex “paneF. The panel extends from 
the blade root to the blade tip and is convected agewise using a rigid wake 
model. The far inboard trailed vortex is placed in the radial center of the 
panel and the far inboard shed vortex is placed on the azimuthal center of 
the panel. The modeling of each of these vortices is similar to the modeling 
of the low resolution inboard far wake. 

The tip vortex for all but the reference blade extends agewise from the 
blade tip to the number of wake spirals specified (see namelist NLWAKE) 
in the low resolution portion of the code. For the reference blade, the tip 
vortex and the inboard far wake begin at an agewise location determined by 
the parameter KNWINT specified in namelist NLHIRES. 


3.3.3 Vortex Segment Location 

The locations of the tip and inboard vortices must be determined at 
azimuthal stations in between the known low resolution results. First, the 
far inboard wake is linearly interpolated to the current azimuthal location, 
t/\ between two known low resolution stations, V’/o an d t/Vn* For example, 
Figure 3.2 shows two known low resolution far wake u p ane l s ’' (which are 
normally represented in CAMRAD.Modl as two vortex elements with large 
core sizes) at the azimuth locations, Wu and The solid lines represent 
two consecutive known low resolution azimuth stations. Note that for each 
low resolution location, the wake age (<z>) starts at zero. The dotted lines 
represent the current high resolution azimuthal location where information 
is unknown. At the current azimuth, </% the vortex location is determined 
by linearly interpolating between two points of equal wake age, <j>. For 
example, the new vortex endpoint coordinates (labeled A) are determined 



by interpolating between the endpoint coordinates at (4 J ki* 0 = 0) and 
(0 /lm 0 = 0 ). All wake endpoint locations of the inboard far wake are 
similarly determined. 

The tip vortex in the far wake is determined by a slightly different in- 
terpolation scheme than that used for the inboard rigid far wake. In order 
for the tip vortex to be translated between one low resolution location and 
another, a scheme shown in f igure 3.3 was implemented. Instead of interpo- 
lating being two wake endpoints with the same wake age, the interpolation 
is done along convection lines of the tip vortex. It can be seen that in the 
inboard far wake interpolation scheme, the trailed and shed vortex lines are 
effectively “pulled* around the azimuth with the blade instead of being “de- 
posited in the wake. I his simplistic model is acceptable for the far wake 
because the influence coefficients are much less sensitive to changes in in- 
board far wake models than to tip vortex models. However, to include more 
physically realistic modeling of the tip vortex problem, the tip vortex is as- 
sumed to be “deposited ” by the blade into the wake. To accomplish this, the 
tip vortex endpoints are interpolated along straight convection lines between 
the low resolution tip vortex endpoints. For example, in Figure 3.3, as in 
Figure 3.2, t he solid lines represent known low resolution azimuth locations 
and the dotted lines represent the current high resolution azimuth location 
desired. Here, the coordinates of the tip vortex are obtained by interpo- 
lating between the known tip vortex endpoint coordinates at the locations 
(Phi* 0=10 ) and (t pi a , 0=0 ). This scheme effectively translates, or 

convects, the vortex from one low resolution endpoint location to the next. 
Note that with this scheme, an addit ional tip vortex segment (from the blade 
tip to point A) must be included to connect the interpolated endpoint at 
location A to the tip of the blade at the location 0. This additional vortex 
segment is included in the analysis. 

This modeling of the far wake and tip vortex is carried out for all blades. 
The far wake of the reference blade begins at the wake age specified by the 
parameter KNWINT in namelist NLHIRES. As a historical note, KNW1NT 
could be set equal to zero to recover an earlier version of the high resolution 
calculation procedure which included only the tip vortex and the inboard 
far wake. Since the tip vortex and inboard far wake in this case begin at 
the reference blade, no near wake model would be required. However, this 
is a crude approximation to the near w r ake. 



Sheet Outside 
Edge 



Edge 


Figure 3.2: Far wake inboard geometry interpolation illustration. 
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3.3.4 Tunnel/Fuselage Corrections 

The tunnel/fuselage correction model is included in determining the tip 
vortex location at a high resolution. The usage of this model is similar to 
the low resolution application that is described in Chapter 2 and no new 
theory is presented here. In the high resolution scheme, the tunnel/fuselage 
additional wake distortions are added to each low resolution wake endpoint 
coordinate before interpolation is made to the current location. Since the 
addition of the tunnel/fuselage components of the additional wake distortion 
takes place before the interpolation, the w r ake geometry components of the 
tunnel/fuselage correction model do not need to be intepolated separately 
to a high resolution. However, the additional velocities over the rotor disk 
are linearly interpolated immediately after input for use in the high reso- 
lution portion of CAMRAD.Modl. All of these processes are internal to 
CAMRAD.Modl and no additional user intervention is required. 

3.3.5 Rollup Model 

The rollup model modifications to the low resolution portion of CAM- 
RAD.Modl were in discussed in Chapter 2. The high resolution implemen- 
tation of the rollup model is only an extension of the low^ resolution imple- 
mentation. No new theory is necessary. As implemented, the rolled-up tip 
and secondary vortex locations are added to the low resolution wake end- 
point locations (as was described for the tunnel/fuselage corrections above) 
before interpolating to the current azimuth angle and wake age in HIRES. 
If the default option for the multi-core model is being used for the tip and 
secondary vortices (/.<., the array of constant size cores), then the low 
resolution multi-core vortex core sizes are used at the low resolution wake 
endpoints as is done for the vortex locations described above. However, 
if the “single, variable core, multi-core model 1 ’ or if the “variable multi- 
core, multi-core model” is being used, the low resolution results for these 
internally calculated core sizes are interpolated to a high resolution in the 
subroutines R0LLURH1 and ROLLUPH2 (for rotor-1 and rotor-2, respec- 
tively). These new interpolated arrays are then used in the calculation of the 
high resolution wake influence coefficients. Again, all of these interpolations 
and connections are made internally in CAMRAD.Modl and no other user 
intervention is required. 
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3.3.6 Vortex Segmentation 

In determining the wake influence coefficients of the high resolution tip 
vortices, it was discovered that under certain circumstances, the tip vor- 
tex endpoints, being connected by straight line segments, aligned such that 
the junctions where the vortex endpoints are connected produced an artifi- 
cial unsteady effect in the aerodynamic loading. These artificial effects were 
reduced by decomposing the offending vortex segments into several new seg- 
ments (see Figure 3.4). To accomplish the decomposition, a criterion was 
set such that if a tip vortex segment is within a 10 azimut h angle of the ref- 
erence blade, then that vortex segment along with the two vortex segments 
on either side of it are subdivided into five segments. The resulting six vor- 
tex endpoint coordinates and strengths are linearly interpolated from the 
original four endpoint coordinates and strengths bv the following formulae: 


q(\) = Q(l) (3.1) 

9(2) = |<?(1)+|q(2) (3.2) 

f /(d) = “Q(2) + -Q(3) (3.3) 

9(4) = ~Q{ 2) + ~Q(3) (3.4) 

9(5) = ^Q(3) + -Q( 4) (3.5) 

9(6) = Q(4) (3.6) 


where q is the new divided segment position or strength and Q is the original 
position or strength. This vortex segmentation is controlled by the param- 
eter OPSEGD in namelist NLHIRES. If OPSEGD — 0, no segmentation is 
performed. If OPSEGD = 1, segmentation is performed. 

3.3.7 Aerodynamic Collocation Point Shifts 

The aerodynamic effects of a swept planform in the high resolution so- 
lution procedure is addressed. With a swept planform (f.c., sweep of the 
quarter-chord fine of the blade), the aerodynamic collocation points are spa- 
tially displaced. Thus, a position change is made to the collocation point 
coordinates in the influence coefficients calculation (see Figure 3.5). 
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Figure 3.4: Illustration of vortex segmentation. 
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Figure 3.5: Swept planform collocation point shifts. 
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No modifications were made to alter the location of the tip vortex due 
to blade tip position changes nor were modifications made to change the 
structural model of the blade. Since the tip vortex location is unmodified 
bv the collocation point shift, care must be used to make sure that the far 
wake of the reference blade begins at a sufficient wake agewise distance from 
the reference blade. The definition of sufficient, depends on the rotor being 
analyzed. For example, if the far wake were to start at the reference blade, an 
aft shifted collocation point would “see” tip vortex elements ahead of the aft 
shifted blade tip since the vortex begins at the original (unshifted) reference 
blade tip. Normally, the far wake does not extend up to the reference blade 
tip OPNECV = 1, or KNWINT > 0). 

If the internal high resolution vortex lattice near wake is used with the 
shifted collocation points, however, the shifted location problem just dis- 
cussed is not an issue during the lattice near wake calculations. This is 
because the near wake lattice model geometry is determined by the blade 
geometry. The shifted collocation point method is not included in the low 
resolution portion of the code. 

The shift in collocation points is a user specified shift in the lead-lag 
direction (parallel to the hub plane) and does not include a. vertical shift. 
The parameters DRPROOT. DRPT1P, and DRP. These parameters are 
found in namelist NLHIRES. The meaning of these variables are that (1) 
DRPROOT is the shift of the blade root (positive aft), (2) DRPTIP is the 
shift of the blade tip (positive aft), (3) DRP(IOO) is the shift of the high 
resolution collocation points (positive aft). 

The equations used to determine the new collocation point coordinates are 
as follows (see Figure 3.5): 


x = xo + Ac sin ip (3.7) 

y = y 0 - Accos t/' (3.X) 

s = -o (3.9) 

where Ac is DRPROOT, DRPTIP, or DRP depending on the current loca- 
tion on the blade. 

3.3.8 Far Wake Loading 

Once the wake influence coefficients are known, they are assumed to be 
invariant throughout the subsequent loading calculations. This is similar 
to that found in the low resolution portion of CAMRAD.Modl. As was 
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seen in Figure 3.1, after the influence coefficients calculation in subroutine 
WKC1INT. subroutine MTNRINT1 is called to calculate the aerodynamic 
loading. This subroutine is a simplified version of the low resolution sub- 
routine MOTNR1. The simplifications are possible since MTNRINT1 is 
not required to be involved in a trim loop. Its only functions are to call the 
blade position subroutine MTNBINTl (the high resolution equivalent to the 
low resolution subroutine M0TNB1), call the airloads calculation subrou- 
tine AEF1INT (the high resolution version of the low resolution subroutine 
AEROF1), and to print out the results to a file. The file write parallels the 
low resolution output file and is described in the following subsection. 

If deemed necessary by the user, the far wake loading information calcu- 
lated may be used to calculate a new circulation distribution which in turn 
could be used to recalculate the induced far wake velocity. This loop may 
be seen in Figure 3.1 as the return arrow from below MTNR.INT1 to above 
CIRCGAL1. Normally, this return path is not used. 

3.3.9 Output Aerodynamic Information 

A high resolution aerodynamic information file is output for each rotor in 
the same format as described for the low resolution aerodynamic information 
file described in Chapter 2 with the following exceptions: 

1. The number of radial stations is MRAINT, 

2. The number of azimuth stations is MPSIINT, 

3. There are three comment lines at the top instead of two, 

4. The moment coefficient is excluded, 

5. The blade flapping deflection is excluded, 

6. The interpolated maximum bound circulation is the value interpolated 
from the low resolution solution, 

7. The value of maximum bound circulation after the high resolution 
calculation is included. 

The unit numbers for the far wake aerodynamic files for rotor- 1 and rotor-2 
are 17 and 57, respectively. 
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3.3.10 Output Induced Velocity Information 

lor each rotor, the far wake induced velocity is written to a file (in 
the subroutine VINDCAL1) at the end of the far wake phase of the high 
resolution calculation. The file contains all three components of the induced 
velocity for all blade sections at all azimuth stations in the following format: 

IF(IWR.EQ . ITERINT) WRITE (18 ,830) 

830 F0RMAT(1X, ’INDUCED VELOCITY (ROTOR l):’) 

DO 850 1=1,3 
DO 850 JJJ=JFIRST, JLAST 
850 IF (IWR.EQ. ITERINT) WRITE(18 ,860) 

1 (VINDINT(I.JR.JJJ) , JR=1 ,MRAINT) 

860 F0RMAT(10F12.6) 

These velocity files are written to unit numbers 18 and 58 for rotor-1 
and rotor-2, respectively. 


3.4 High Resolution Lattice Near Wake Model 

3.4.1 Introduction 

As discussed in the previous section, the far wake consists of the bound 
vortices of all blades except the reference blade, the tip vortices of all blades, 
and an inboard rigid wake of all blades. The far wake does not start, im- 
mediately from the blade but is offset in wake age bv the value KNWINT 
(see namelist NLH1RES). From the blade to the first far wake elements, 
a near wake model is employed. The near wake model can be viewed as 
a transfer function as shown in Figure 3.6 and the far wake aerodynamic 
loading can be viewed as a forcing function for the near wake model. The 
near wake transfer function has been formulated by two separate means: 
(1) a vortex lattice model, which uses, among other things, the far wake 
aerodynamic loading (circulation) information, and (2) an external “lndi- 
cial Post -Processor” (IPP) code, which uses the total far wake velocities, the 
wake induced velocity, and the blade pitch at all blade radial and azimuthal 
stations. The vortex lattice near wake model is internal to the high reso- 
lution part of CAMRAD.Modl and is discussed in this section. The IPP 
code is external to the CAMRAD.Modl code (but is still part of the system 
called “HIRES”), and is discussed in a later chapter. Either of the two near 
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Figure 3.6: Schematic of near wake transfer function. 

wake models may chosen, but not both; however, note that the lattice near 
wake model has not been exercised or validated. 

In this section, the near wake lattice model is presented. The model 
uses as input, some portion of the far wake information and as an output a 
near wake induced velocity which when combined with the far wake induced 
velocity may be used to calculate the aerodynamic loading on the reference 
blade. The latt ice near wake loop involves calculating influence coefficients of 
all wake elements in the near wake lattice, calculating the velocity induced at 
blade collocation points by these vortex elements, and calculating the total 
aerodynamic loading (and circulation) at. the blade collocation points with 
the newly calculated near wake induced velocities and previously calculated 
far wake induced velocities. Each aspect of this model is discussed in detail. 

3.4.2 Lattice Geometry 

The high resolution near wake lattice is comprised of vortex elements 
placed on wake “panels behind the reference blade. These panels extend 
in a wake agewise direction which is perpendicular to the reference line 
(the unswept, straight quarter-chord) of the reference blade, as shown in 
Figure 3.7. The “side edges" of the panels are defined by the location of 
the collocation points on the reference blade. The edges of the panels in the 
radial direction are determined by the radial shape of the blade. Therefore 
the panel shape is always determined by the blade shape (radially, vertically, 
and azimuthally). So, if the collocation points are shifted due to sweep, as 
discussed earlier, the near wake is adjusted so that it is always attached 
to the blade. The equations used to determine the near wake panel edge 
locations, x nw , y nw , and z nw , associated with each radial station, i, are as 
follows: 

*nw(i) -- r nw (i) cos( t’ - <j>) 

Vnw(i) = i' nw (i) sin(V’ - <t>) 

~mr ( *) — ~b ( C. li' <f>) 
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(3.10) 

(3.11) 

(3.12) 






(3.13) 


l'nw{>) = \f x l + y'h + ~fc 

where y and £5 are the coordinates of the blade at the i-th radial 
station, V' is the reference blade azimuth location, and (p is the wake age. 
This calculation is performed for all wake elements of the ages 4> — 0 to 
KNWlNT*At\ 

With the panel edge locations and therefore the panel corner locations 
known, the vortex elements representing the panel may be located. Each 
panel’s trailed vortex is located in the spanwise center of the panel. That, is, 
the trailed vortex on each panel is located by connecting the two midpoints 
of the leading edge of the panel and the trailing edge of the panel (see Figure 
4.8). The shed vortex on each panel is found by determining the fore-aft 
midpoint on each panel side edge. The shed vortex on each panel is located 
by connecting the two midpoints just determined, then shifting the shed 
vortex aft by a distance of one quarter chord (c/4). At the root arid tip, the 
trailed elements are placed in the middle panel, as for all the others, with 
the outside edge being defined by the end of the blade (whether root or tip). 
The equations used in determining the vortex locations are as follows: 


(^2 + i'4) 
VA = 2 

(3.14) 

(n + r 3 ) 
’■» = 2 

(3.15) 

(fij + r-O 

2 | r 4 - r 3 | 

(3.16) 

(^2 + fi) (^2 - n) * <k+\ 
2 \ri - n| 

(3-17) 


where r are position vectors to the locations denoted bv the subscript, and 
d t is the local chord divided by four. 

3.4.3 Vortex Strengths 

The strength of each vortex element (trailed and shed) must be deter- 
mined. The trailed vortex strength is calculated from the radial derivative 
of circulation and the shed vortex st rength is calculated from the azimuthal 
derivative in circulation. The change in circulation is found by differencing 
the bound circulation across a particular panel. For example, the bound 
circulation at each point A, B, C, and D are found by the following: 
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Figure 3.8: Close-up of near wake panel geometry. 


} .-A 7 2 74 

(3.18) 

7 B = 71 - 73 

(3.19) 

7 c = 7 4 - 73 

(3.20) 

7 D = 72 - 7i 

(3.21) 


where '71,72,73, and 74 are the circulations at the corners of the panels. 
These circulations are equal to the blade bound circulation at the azimuth 
for which the vortex endpoint was “released” from the blade. The circulation 
distribution between points A and B and between points C and I) may 
optionally be a stepped distribution with a. step in circulation at the panel 
midpoint or may be instead a linear circulation distribution. The vortex core 
model may also be either a distributed core model or a concentrated core 
model. These options are controlled by the namelist NLHIRES variables 
MDLSNW, MDLTNYV, OPCSNW, and OPCTNW. The core size of the 
t railed and shed vortex elements may be an input constant for each direction, 
or may be calculated internally as controlled by the NLHIRES variables 
CORETNW and CORESNW. 


3.4.4 Total Loading 

With the trailed and shed wake geometry known, the vortex core size 
known, and the type of core (none, concentrated, or distributed) known, the 
influence coefficients of the trailed and shed elements on a particular blade 
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collocation point may be found. The influence coefficients are then mul- 
tiplied by the appropriate circulation, as determined above, to obtain the 
velocity contribution from that vortex element acting at the blade colloca- 
tion point. As the procedure continues, all velocity contributions of all near 
wake vortex elements for a blade at a fixed azimuth are determined. These 
velocities are added to the appropriate far wake velocities to determine a 
total velocity at the blade collocation point due to all wake elements. Once 
the total velocity is known for a blade at the current azimuth, the aerody- 
namics at that azimuth are determined. The azimuth is incremented and 
the process is repeated. Normally, it is expected to take more than one rev- 
olution (ITERNW > 1) to converge the near wake solution since during the 
first revolution (iteration), the starting transients for the near wake solution 
must be allowed to decay to an acceptable level. 


3.4.5 Option to Reduce Number of Panels 

NLHIRES contains several near wake option parameters. There is an 
option to reduce the number of panels in the near wake as a predetermined 
function of the wake age, 0. In principal, this technique should reduce t he 
computation time in the lattice near wake model. By default, the parameter, 
OPNWMIN = 0, which suppresses the reduction in the number of panels. 
If OPNWMIN —1, the number of panels is reduced by a factor of two every 
5 of wake age. (Actually, the number of panels is exactly one-half the 
previous number if the previous number of panels is even and the number 
is one half of the previous number of panels + 1 if the number of panels is 
odd.) For example, if there are SO panels at the wake age 0=0 degrees, 
then starting at 0 = 5 there will be 40 panels, and at 0=10 there will be 
20 panels, etc. The reduction in panel number will continue until either t he 
minimum of two panels is reached or until the end of the near wake region 
is reached, whichever comes first. 

3.4.6 Circulation Update Option 

There exists an option to update the circulation being used with the 
wake influence coefficients to calculate near wake induced velocities at blade 
collocation points. The NLHIRES parameter OPNWCRC is used to control 
this updating. The updating may be performed at the end of each near wake 
iteration (OPNWCRC = 1) or may be updated at each azimuth step as it 
is calculated (OPNWCRC = 0). The default is to update at each azim utli 
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step as the circulation is calculated. This default should reduce the number 
of near wake iterations required. 

3.4.7 Output Information 

Output from the lattice near wake section of the high resolution calcula- 
tion is similar to the outputs from the low resolution calculation and from 
the far wake calculation. Again, there is one file output for each rotor. The 
unit numbers for rotor-1 and rotor-2 are 19 and 59, respectively. The for- 
mat of each file is the same as the far wake file described earlier with the 
following exceptions: 

1. The vertical component (z coordinate) of total induced velocity is in- 
cluded as an additional entry after the MAXIMUM BOUND CIRCU- 
LATION - CIRCN values, 

2. The vertical component (z component) of the high resolution near 
wake induced velocity is included as an additional entry at the end of 
the file. 

The formats of the above items are as follows: 

WRITE(19,895) 

895 FORMAT (IX, ’INDUCED VELOCITY - Z COMPONENT’) 

DO 896 1=1 ,MRAINT 

896 WRITE(19,700) (VINDT0T(3,I , J) , J= JFIRST , JLAST ) 

WRITE(19,897) 

897 FORMAT (IX, ’NEAR WAKE INDUCED VELOCITY - Z COMPONENT’) 

DO 898 1=1 .MRAINT 

898 WRITE(19 ,700) (VNW(3 ,1 , J) , J= JFIRST , JLAST) 

3.4.8 Known Caveats 

It should be noted that the high resolution lattice near wake model discussed 
in this section has not been exercised or validated. 



Chapter 4 

Indicial Post-Processor 


4.1 Indicial Post-Processor 

4.1.1 Introduction 

As discussed in previous chapters, conceptually, the user has an option 
as to which high resolution near wake model to use: (1) the lattice near 
wake model internal to the HIRES code, or (2) the Indicial Post- Processor. 
Chapter 3 discusses the high resolution lattice near wake model that is in- 
ternal to the HIRES portion of CAMRAD.Modl. This chapter discusses 
the Indicial Post- Processor (IPP). The IPP is a stand-alone code that uses 
output information from HIRES (from the “Far Wake Loop” only). This 
far wake information from HIRES is combined with an indicial aerodynam- 
ics method for the near wake to calc ulate the total aerodynamics on the 
reference blade at a high azimuthal and radial resolution. The IPP uses a 
combination of techniques presented in several publications by T.S. Beddoes 
and G. Leishman (Ref. [13, 14, 15, 16]). 

4.1.2 Solution Procedure 

Conceptually, the IPP calculates the unsteady aerodynamic loading at 
each high resolution collocation point on 1 he reference blade using the ve- 
locities at each collocation point (as determined in the Far Wake Loop of 
HIRES) as a gust field. That is, in the Far Wake Loop of HIRES, velocities 
due to blade rotation, velocities due to blade motion, and velocities induced 
by the "Tar-wake vortex system are computed. The IPP then calculates the 
unsteady aerodynamic loading (response) on the reference blade as if there 
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were a, gust field (given by the velocities from the Far Wake Loop of HIRES) 
traveling past the stationary reference blade. 

In the indicia] aerodynamics method, the high resolution near wake is 
divided into two parts: (1) the shed near wake, and (2) the trailed near wake. 
(Note that dividing the near wake in this manner is a common practice (Ref. 
[11]).) Each of these components of the near wake (shed and trailed) are 
addressed in the IPP. That is, the IPP combines the far wake velocities from 
HIRES with the velocities from the near wake trailed vortex system model, 
and uses those velocities to calculate a quasi-steady angle of attack and a 
Mach number. The indicial formulations of Beddoes and Leishman are then 
used to relate the quasi-steady angle of attack and Mach number at a blade 
section to a blade sectional loading. The indicial formulations essentially are 
empirical “curve-fits’" representing the theoretical and experimental relations 
between quasi-steady angles of attack and Mach numbers. This chapter 
discusses in detail how the HIRES far wake information is “processed’" by 
the IPP to produce high resolution loading. 

4.1.3 Conceptual Program Outline 

Figures 1.1 and 1.3 illustrate the computational flow of the CAM- 
RAD. Modi code. It is seen (as described in earlier chapters) that the low 
resolution (azimuthal and radial) CAMRAD.Modl is executed first, followed 
by the HIRES portion of the code. Once these two steps are complete, the 
user may choose to execute the IPP to obtain high resolution loading at the 
high resolution collocation points on the reference blade. Note that the IPP 
is a true post- processor; that is, it is a stand-alone code. Even so, it is still 
considered part of the CAMRAD.Modl /HIRES code system. 

Figure 4.1 expands the IPP box of Figures 1.1 and 1.2 so that the IPP 
may be outlined in a conceptual manner. It is seen that there are three 
inputs to the IPP. The first input is a namelist called INLST. This namelist 
is used to set several parameters used in the IPP. This namelist will be 
discussed later. Second, the far wake information from HIRES is input. This 
information act ually consists of two data files: (1) the far wake data file, and 
(2) the far wake induced velocity file. These two files will be discussed later. 
Third, the binary airfoil table file from CAMRAD.Modl is input. This 
input, too, will be discussed later. 

The IPP consists of two “loops”. The outer loop is a loop over all of 
the high resolution azimuth stations. At each azimuth inside the outer loop, 
there is a loop over all of the collocation points on the reference blade. As 
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Figure 4.1: Schematic of IPP computational loops and inputs. 

discussed earlier, the IPP calculates the loading on the reference blade as 
if it were traversing a gust field given by the far wake information from 
HIRES. I herefore, the reference blade is the only blade considered in the 
IPP. The loop over the collocation points on the blade is where the sectional 
blade loading is calculated. Once the loading is known on the entire span 
at the current azimuth station, a Trailed Wake Algorithm (TWA) is used 
to calculate the influence of the trailed near wake system. It should be 
noted that the influence of the trailed near wake system needs to be known 
before the blade loading can be calculated. But, since the TWA calculates 
the trailed near wake influence based on the blade loading, an interation is 
required. The iteration technique implemented is to lag the TWA calculation 
behind the blade loading calculation by one azimuth station. Thus, once 
the blade loading has been calculated in the loop over all blade collocation 
points, the TWA calculates the influence of the trailed near wake system, 
then uses that information at the following azimuth station. 

4,1.4 Actual Program Outline 

Figure 4.2 shows the outline of subroutines in INDICIAL, which is the 
executive program of the IPP. The first four subroutines are basically ini- 
tialization subroutines, that will be discussed later. The last subroutine 
in Figure 4.2 (CLCALC) is where the “loop over all of the high resolution 













IPP 



Initialization 


Figure 4.2: Major subroutines in the IPP. 

azimuth stations 1 ' (see Figure 4.1) is located. This subroutine will also be 
discussed later. The first stage in the execution of the IPP is initializa- 
tion of parameters. This will now be examined by examining each of the 
initialization subroutines individually. 

4.1.5 Subroutine INPTRD 

The first subroutine in the IPP, called INPTRD, initializes several 
namelist input parameters to default values and reads the input namelist, 
INLST. These input parameters are user defined quantities and are case de- 
pendent. Each of the namelist input variables are discussed later (along with 
their default values). The subroutine INPTRD also checks several input pa- 
rameters for valid values. For example, the number of azimuth stations, 
NAZM, is compared to the maximum number of stations (max. = 720); the 
number of radial stations, NRAD, is compared to the maximum number of 
radial stations (max = 100); and the maximum number of revolutions in the 
solution (/.e., the number of '‘outer azimuth station loops 11 ) are compared 
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to the maximum number of steps (max = 1440). If errors are detected dur- 
ing these comparisons, the code ends execution and prints an error message. 
Since one of the input variables is the file name of the CAMRAD.Modl bi- 
nary airfoil table to be used for this case, the IPP next calls the subroutine 
INPTA1. 


4.1.6 Subroutines INPTA1 

One of the input parameters in the namelist INLST is the file name of the 
CAMRAD.Modl binary airfoil table to be used in this execution of the IPP. 
To read this unformatted binary airfoil table, the IPP calls the subroutine 
INPTA1. This subroutine is identical the the subroutine INPTA1 from 
CAMRAD.Modl and stores the necessary information in arrays in the IPP 
code. 

4.1.7 Subroutine RDFARW 

Two other input parameters are (1) the file name of the far wake informa- 
tion file from HIRES and (2) the file name of the far wake induced velocity 
file from HIRES. Both of these files are read using the subroutine RDFARW. 
As discussed in the “Output Aerodynamic Information” subsection of Chap- 
ter 3, there are many parameters output in the far wake information file. 
The information of interest in this file is (1) blade section velocities parallel 
to the hub plane, (2) the blade section velocities perpendicular to the hub 
plane, (3) the blade section radial velocity, (4) the pitch angle of the blade 
section, and (5) the blade radial collocation point locations. RDFARW reads 
the far wake information file, extracting the above five quantities while skip- 
ping over the information in the file that is not required by the IPP. After 
required information from the far wake information file has been extracted, 
the file is closed. 

Next, RDFARW opens the far wake induced velocity file output from 
HIRES. This file contains the three components of wake induced velocity 
from the far wake vortex system from the Far Wake Loop of HIRES. For the 
IPP, only the vertical (perpendicular to the hub plane) component of this 
velocity is required. Thus, RDFARW extracts only the vertical component 
of far wake induced velocity from this file. Once that is done, the file is 
closed. 

Finally, RDFARW converts all of the newly extracted non-dimensional 
quantites into dimensional quantities for use in the IPP. 



4.1.8 Subroutine UNUC 


The input information read from the far wake information file from 
HIRES is in the hub coordinate system used by CAMRAD.Modl. For exam- 
ple, velocities at each blade section are given in components that are parallel 
and perpendicular to the hub plane. However, the IPP uses velocities that 
are parallel and perpendicular to the blade chord. Subroutine UNUC' uses 


the far wake aerodynamic velocity and pitch information to calculate the 
velocities parallel and perpendicular to the chord of the blade section. Fig- 
ure 4.3 shows the relations among the velocities Up, Up, f '/v, and Uc • The 
following equations are used to calculate Ujp and Uc from Up, Up, and V z \ 


Ufot — \J Uj + Up 

(4.1) 


Utotl = \JU't + (Cp - V;) 2 

(4:2) 


Uc = Utoi cos(0 - <t>) 

( 4.3) 


M = 

SOUND 

(4.4) 


1 

'cc 

o 

II 

(4.5) 


Uj is the velocity parallel to the hub plane at the blade section pUp is 
the velocity perpendicular to the hub plane at the blade section, Utot is 
the total velocity at the blade section, Vtot 2 the total velocity excluding 
the vertical component of induced velocity, V s is the vertical component of 
induced velocity, Uc is the component of the total velocity parallel to the 
chord, l J i\ is the component of the total velocity (excluding the vertical 
component of induced velocity) perpendicular to the chord, M is the Mach 
number calculated from Uc, SOUND is the input speed of sound, 0 is the 
blade pitch angle, and 0 is the inflow angle calculated from Uj and Up at 
the blade section. 

Since the quantities Uj and Up were obtained from a lifting line aerody- 
namics code, it is consistent to use such approximations as are introduced 
in the previous equations. With these approximations, Figure 4.4 illustrates 
the manner in which the IPP views the aerodynamic far wake environment 
in which the quarter chord advances in space (Lc>, 0). Using this figure, 
the aerodynamic environment may be seen as a reference point (here, the 
quarter chord) traveling in a local free stream velocity, Uc, with a perpen- 
dicular velocity determined by Un and V~. This is equivalent to the blade 
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Figure 4.3: Schematic of coordinate system used in subroutine UNUC. 
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V. 



Figure 4.4: Velocities seen by blade section as it travels through gust field. 
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remaining stationary and allowing the far wake velocities to convert past 
the blade as a gust field. 

An additional quantity calculated in UNUC is the pitch rate at all blade 
sections and azimuth locations. This is calculated using a central difference 
formula: 


HiJ) 


0(iJ+])-6(i,j-]) 

2A t 


( 4 . 6 ) 


where 0 is the pitch rate, 0 is the pitch angle, index / is the radial location, 
index j is the azimuthal location, and A/ is the time step between azimuth 
locations. 


4.1.9 Subroutine CLCALC - Introduction 

The aerodynamic loads (forces and moments) are determined in 
CLCALC-. These forces and moments are determined using the indicia! aero- 
dynamic functions in the mentioned References. These indicia! aerodynamic 
functions use an effective angle of attack at a particular instant and relate 
it to the loads at that instant. The loading at any particular instant is split 
into several parts: (1) circulatory loading, (2) impulsive loading, (3) loading 
due to trailing edge separation, and (4) loading due to leading edge sepa- 
ration. The subroutine CLCALC calculates the circulatory loading effects 
and calls other subroutines to account for the other effects (these will be 
discussed later). 

Conceptually. CLCALC 1 uses the information calculated in UNUC and 
applies indicia] aerodynamic functions to the information at each blade sec- 
tion and azimuth to determine the aerodynamic forces and moments (circu- 
latory loads) there. (If the optional impulsive loads are included, CLCALC 
calls IMPS to calculate the impulsive loading effects.) In practice, both the 
circulatory and impulsive (non-circulatory) effects are two dimensional ef- 
fects; the three-dimensional effects of the trailed near wake, a Trailed Wake 
Algorithm (TWA) similar to that of T.S. Beddoes (Ref. [13]) is used to 
model the effects of a trailed near wake. These two dimensional effects are 
split into circulatory and non-circulatory effects. 

To calculate the loads (circulatory and non-circulatory), effective angles 
of attack are required (see mentioned References). These effective angles 
of attack are calculated in the subroutine INTGRL (discussed later). This 
subroutine ret urns (to CLCALC) the effect ive angles of attack for the circu- 
latory lift and moment and for the non-circulatory lift and moment. How- 



ever, these angles of attack calculated in INTGRL are quasi-steady angles 
of attack. The two dimensional indicia] aerodynamic functions modify these 
quasi-steady angles of attack to account for unsteady effects. CLCALC 
accounts for the circulatory loading effects. Then, the subroutine IMPS ac- 
counts for t he noil-circulatory effects. Once these are known, the subroutine 
SEPRATE accounts for the leading edge and trailing edge separation effects. 

4.1.10 Subroutine CLCALC - Coding 

The outer loop in the subroutine is over the number of azimuth steps 
defined by the namelist input variables NAZM and NREVS. The number 
of steps used in the calculation is NAZM*NREVS. The maximum number 
of steps is 1440. At each step of the outer (p) loop (see Figure 4.5), all 
quantities are calculated at each radial station on the blade. 

As stated above, the quasi-steady angles of attack are determined in 
INTGRL. Then, CLCALC- modifies these quasi-steady angles of attack to 
account for unsteady circulatory effects. From the mentioned References, the 
following formulae are used to account for the unsteady circulatory effects 
in an indicia! form: 


C'eff = >]L ~ A -Y - Z 

(4.7) 

x = + AAm. - VL M )f-°^ l2Tl 

(4.8) 

Y = W' TA,/Ta + Mm, - 

(4.9) 

z = ZoUt~ (,AtlT ’ + A 3 ( Vl - i,L old ) e-** l2T » 

(4.10) 

v = <rw' TA(/r ”' + .4 - vm old )(~ aAt/2Tm 

(4.11) 

a = -|L_(1 - A/ 2 ) 

chord 

(4.12) 

.4, = .165 

(4.13) 

A-2 = .335 

(4.14) 

4 3 = .500 

(4-15) 

A m = 1.00 

(4.16) 

T\ = 20. 

(4.17) 

T 2 = 4.5 

(4-18) 

T 3 = 1.25 M 

(4-19) 
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Done 


Figure 4.5: Flowchart for subroutine CLCALC. 
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T m = M/2 (4.20) 

where i)l and t] m are the quasi-steady circulatory angles of attack calculated 
in the subroutine INTGRL. Once the angle of attack a f}J is known, it is 
then modified in a manner similar to that used in CAM RAD. Modi to model 
the yawed flow and blade sweep effects on the aerodynamics. Next, these 
modified (e.g., modified by the sweep correction model) angles of attack 
and modified {(.y., modified by the sweep correction model) Mach number 
are used to calculate the lift and drag coefficients C'i and Cj bv airfoil table 
interpolation. The (7/ and Cj values are used to calculate a normal force 
coefficient, C n and the chordwise force coefficient, C c . For a reference value 
to be used later, a C n , pot( ntial is calculated next by the following formula: 


t - 71 potential 


2 net 


?!S 


y/1 - M 2 


+ Oo 


(4.21) 


The moment coefficient, C m , is then calculated from the indicial aerodynam- 
ics (circulatory effects). The subroutine IMPS is called next to evaluate the 
impulsive (non-circulatory effects) normal force, chordwise force, and mo- 
ment coefficients, (),.€'</,, and C m ,. The impulsive moment is then added 
to the C m value. The subroutine SEPRATE is called next to calculate the 
effects of trailing edge and leading edge separation on the lift, drag, and 
moment coefficients. (The impulsive lift and drag terms are added to the 
previous lift and drag coefficients inside SEPRATE.) Upon return from sub- 
routine SEPRATE, subroutine TWA is called to calculate the effects of a 
trailed near wake due to the current aerodynamic loading calculations. 

At the end of each azimut h step in CLCALC, the following quantites may- 
be out put (depending on the input parameter ILOAD): (1) normal force and 
chordwise force, or (2) C' n M 2 . The forces are expressed in Newtons and the 
C n and M are local values. The local C„ is determined from the following 
definition: 


N = p(Ma) 2 cC n (4.22) 

.V is t he local normal force per unit span, M is the local Mach number, 
r is the local chord, a is the speed of sound, C n is the local normal force 
coefficient. 
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4.1.11 Subroutine INTGRL 


As seen in Figure 4.5, the first task in CLCALC after initialization is 
to invoke INTGRL to calculate the effective quasi-steady angles of attack 
(using a weighted integral over the chord). These are later to be used in the 
circularory lift, circulatory moment, impulsive lift, and impulsive moment 
calculations. The following integrals are used to determine the effective 
angles of attack: 




MO) 

U c 


(1 — cos 0)d0 


Xl ~ 



MO) 

Uc 


sin 2 0d0 


cos 0=1- 



(4.23) 


(4.24) 


(4.25) 


where i]i is the effective angle of attack at the 3/4 chord point for the 
circulatory terms, Xi is the effective angle of attack for the impulsive terms, 
0 is the integration coordinate, c is the chord, x is the chordwise location, 
w(0) is the downwash at the 0 coordinate, and Uc is the chordwise velocity 
at the 0 coordinate. 

To compute these integrals, the downwash velocity w(0) and Uc as a 
function of 0 must be found. CAMRAD.Modl is not capable of providing 
these intantaneous velocity distributions over the chord. Thus the velocities 
are approximated velocities over the chord. It is assumed that the chord at 
any radial station “spans” several points where the velocities are known. For 
example, Figure 4.6 illustrates a chord at a particular radial and azimuth 
station. Note that the azimuthal station is defined to be the quarter-chord 
location of the section. In this figure, the chord is seen to ^span” several 
other azimuth stations which are a distance Hr At apart (where Q is the rotor 
rotation rate, r is the radial station location, At is the time required for the 
rotor to travel from one high resolution azimuth station to the next). The 
velocities from these other azimuth stations are used in the above integrals. 
However, since the integration is over the entire chord, two additional points 
must be introduced to complete the calculations. In general, the velocities 
at the leading and trailing edges of the blade chord are not known. To 
evaluate the velocities at these points, a linear interpolation between the 
two surrounding points is used. Thus, the velocities at points on the chord 
(including the leading and trailing edges of the chord) are known. 
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The integrals of Equations 4.23 and 4.24 are evaluated using a trape- 
zoidal rule which has been modified to include subintervals in the integra- 
tion. A coordinate system in 0 is used as shown in Figure 4.6. It has been 
determined through numerical experimentation that at least three points on 
the chord are required for the integral evaluations to be sufficiently accurate. 
This is sometimes a problem, particularly near the tip of the blade when 
the azimuthal resolution is inadequate. If there are less than three points 
defining the chord, the IPP issues a warning message and then stops exe- 
cuting. (The corrective action is to re-execute CAMRAD.Modl and HIRES 
with an increased azimuthal resolution.) 

An additional n on-circulatory term may optionally be employed. This 
term is not included in the references mention. To simulate effects of blade 
vortex interations (BVI), it is assumed that the impulsive character of the 
BVI is strong at the leading edge of the blade section, but decays toward the 
trailing edge. An approximation to this assumption is to allow an additional 
“step” function over the blade section. The magnitude of this function 
is scaled by the user. The extent of this function over the chord is also 
controlled by the user through input variables (IWT, WEIGHT, XOCOFF). 
The default values set in the code should be sufficient for many cases. These 
input variables are discussed in a later section. 

4*1.12 Subroutine IMPS 

The equations coded in subroutine CLCALC are the circulatory lift and 
moment terms. These are terms related to the circulation present on the 
airfoil and in the wake of the blade section. They serve to modify the 
quasi-steadv angles of attack to account for shed wake effects. To include 
non-circulatory effects, another term is required; these terms are calculated 
in the subroutine IMPS. The equations coded in IMPS for the impulsive 
loading are as follows: 

(4.26) 

(4.27) 

(4.28) 

(4.29) 


= M Hn 
_ -4 

Cm ' = Jf Hm 


Cc 


At 


H r 


Hn = tf« oid e~ Af/T,A ' + (A / - \i n old)< - At,2T ' Kl 
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Figure 4.6: Coordinates used in the integration subroutine INTGRL. 


Urn 

= H,n oht i A ' /7,Al + (A m - \ m „old)( A '/ 2T ' /v '< 

(4.30) 

H c 

= H Cold e~^T,K t + Kc(Uc _ r ^ rf)f .-A ( /2T,/w 

(4.31) 


. . T7T 

l\ = 

AU e 

(4.32) 


T , _ c 1 + 3 M 

* l — * 

a 4 

(4.33) 


I\'i = 1 - e 4y H 

(4.34) 


where y is the radial distance from the blade tip to the current collocation 
point, t is the thickness ratio of the airfoil, a is the sound speed, and c is 
the chord of the current blade section. These additional terms are added to 
the circulatory terms discussed above. 

4.1.13 Subroutine SEPRATE 

Once the effective angles of attack from subroutine 1NTEGRL and the 
impulsive lift components from subroutine IMPS are known, the effects of 
trailing edge and leading edge separation are taken into account. This is 
done in the subroutine called SEPRATE. Many of the parameters involved 



SERRATE are, like the bulk of the indicia! work, governed by empirical 
factors. This subroutine is an interpretation of the mentioned References. 
Since most of the empirical factors in this subroutine are functions of Mach 
number, the first task in SEPRATE is to linearly interpolate these parame- 
ters to the current local Mach number. Once these parameters are known, 
the t railing edge separation effect may be calculated. The trailing edge sepa- 
ration effect is discussed in detail in the mentioned References. In short, the 
trailing edge separation effect is a lag in the surface pressure (and thus lift) 
with respect to the angle of attack of the blade section due to the motion 
of the separation point near the trailing edge. In the references, the t railing 
edge separation effect is calculated from empirical formulae known a priori 
. In this application, the location of the trailing edge separation point is 
calculated by the following equation: 


f = 

J n 

where /' is the effective separation point location, C n $\ is the value of 
normal force coefficient calculated from an airfoil table lookup previously 
made in subroutine CLCALC (this term also includes the impulsive normal 
force coefficient), C n%po ttnti<il is the value of the potential lift coefficient, also 
calculated previously in subroutine CLCALC. With the effective trailing 
edge separation point known, the motion of the separation point is lagged 
using equations 21 and 22 from the mentioned References. Once the lagged 
location of the trailing edge separation point, f ff is known, it may be re- 
applied to the following equation: 


. 


c n . 


,81 


C ii'potential 


(4.35) 
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n,potfntial 


i + v/T 77 ' 


"f C- n i 


(4-36) 


where the normal force coefficient., C' n j. includes the efTect. of trailing edge 
separation, and the effect of the impulsive normal force, C, M . 

The effects of leading edge separation includes the additional lift on an 
airfoil due to leading edge vortex separation. A simplified version of the 
model presented in the mentioned References is used. The simplification is 
that only one is allowed to be shed in the stall region. For blade vortex 
interaction (BVI) calculations, this is not a severe limitation as B Vis do 
not often appear in the stall region of the rotor. The effect of leading 
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edge separation on the normal force coefficient is added to the normal force 
contribution from the trailing edge separation. 

Now, with the total normal force known, other quantities such as chord- 
wise force coefficient and moment coefficient are calculated using the equa- 
tions of the mentioned References. 

4.1*14 Subroutines TWA 

The subroutine TWA uses the same equations and serves the same pur- 
pose as the subroutine TWA1 in CAMRAD.Modl (see Chapter 2 for details). 

4.1.15 Subroutine FINDNN 

The subroutine FINDNN is a subroutine that determines the current 
azimuth index that is between one and the number of azimut h steps, NAZM, 
given an arbitrary azimuth index between one and NAZM *N REVS. 

4.1.16 Subroutine AEROT1 

Subroutine AEROT1 interpolates the input airfoil table information. It 
is the same as the subroutine AEROT1 from CAMRAD.Modl. 


4.2 Indicial Post-Processor Namelist Input Vari- 
ables 

4.2.1 Namelist INLST 

The parameters in INLST are given in Table 4.1. 
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Table 4.1: Description of INLST Parameters 


FWFILE 


name of far wake file (character* 80) 

VINDFILE 


name of far wake induced velocity file (charac.ter*80) 

NR AD 


number of radial stations used in calculation (max 100) 

NAZM 


number of azimuth stations used in calculation (max 720) 

SOUND 


speed of sound (m/s or ft/s) 

RPM 


rotor RPM 

RADIUS 


rotor radius (m or ft) 

CHORD(NRAD) 


blade chord (m or ft) 

NREVS 


number of rotor revolutions to calculate (default = 2) 

RAE(NRAD) 


edges of radial segments (usually same as high resolution 
variable RAEINT) 

icurv 


integer parameter: 


= 0 

use straight trailed near wake 


= 1 

use circular arc trailed near wake (default) 

EP1C0R 


trailed near wake core size: 


> o. 

use this constant core size for all trailed lines (re. R) 


< 0. 

use .5* (panel width) 
(default — —1.) 

THICK 


blade thickness used in impulsive chordwise loading term 
(percent of chord, in decimal form) 

(default = .12) 

SWP(NRAD) 


sweep angle at radial segment locations (deg) 
(default = 100*0.) 

DENSE 


density ( kg/m 3 or slug/ ft 3 ) 

IWT 

= 1 

use constant “step 1 ' function controlled by WEIGHT and 
XOCOFF (default) 


= 2 

use new weighting function (not recommended) 

WEIGHT 


magnitude of impulsive BVI term (default = 1.0) 

XOCOFF 


x/c cut toff location for BVI term (default = 0.25) 

I LOAD 

= 1 

output normal force (default) 


= 2 

output C n M 2 

IVZ 

= 0 

zero-out V z in impulsive loading term 


= 1 

do not zero-out l ~ in impulsive loading term (default) 
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Chapter 5 

Users’ Manual: Variables 
and Namelists 

5.1 Introduction 

I his users manual is intended to contain new input variables and 
namelists for CAMRAD.Modl. The intent is to document the new namelists 
and variables with respect to the original version of CAM RAD (e.g., where 
the new namelists are located relative to the older namelists, etc). For a de- 
scription of all older namelists and variables, refer to [2]. In addition, since 
the current version of CAMRAD.Modl is intended to run on a workstation 
with a UNIX operating (more specifically, it was set up to execute on a 
DEC- Alpha workstation), this portion of the documentation also describes 
methods to execute the code in a UNIX environment (specifically, on a DEC 
Alpha workstation). 

Even though many changes have been made internally to create CAM- 
RAD.Modl from the “base 7 code, CAMRAD, the basic trim loop struct ure 
of the code is most ly unchanged. The changes to the loop structure that have 
occurred are the addition of extra loops. For example, a new “rollup-trim” 
loop was added to execute at the end of the original three stage wake-trim 
procedure. Also, both a CFD interface and a high resolution procedure 
(HIRES) were added to execute at the end of the entire trim procedure. For 
more information on these new procedures see previous chapters. 

Since the goal of the CAMRAD.Modl code is to obtain high resolution 
loading information that can be used in conjunction with at helicopter noise 
prediction code, the remaining original procedures in CAMRAD.Modl (Hut- 
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ter, flight dynamics, and transient) have not been updated, maintained, or 
tested. In some cases, they are not compatible with the new procedures 
introduced in CAMRAD.Modl. However, if one were to execute CAM- 
RAD. Modi in a manner similar to that of the original version of CAMRAD, 
the flutter, flight dynamics, and transient tasks should still be available - 
however, this has not been verified. 


5.2 Summary of Job Preparation 

A typical job preparation is basically unchanged from the original CAM- 
RAD procedure. The first step is to create a BLOCKDATA file containing 
input information for CAMRAD.Modl. Normally, this information is then 
compiled into a binary input file in a manner similar to the original version 
of CAMRAD. An example script for generating a binary input file is pro- 
vided in the "Binary Input File Preparation” section of this manual. Next, 
airfoil data, must be provided. This is available either by converting C81 
tables into a binary airfoil table using the airfoil preparation program or by 
converting generic airfoil information into a binary airfoil table generated 
with namelist inputs to the airfoil preparation program. An example of an 
airfoil preparation script is given in the “Airfoil Table Preparation” section 
of this manual. Once these two binary files are created, the code is exe- 
cuted using a script file. A script template is given in the “Script Template” 
section. 


5.3 Airfoil Table Preparation 

Airfoil table preparation is normally accomplished by creating a “script' 
file to run the airfoil preparation program. An example script for generating 
the binary airfoil table to be used in CAMRAD.Modl from an input C81 
table is as follows: 

! /bin/csh 

# 

# For Bol05 model rotor 

# NACA 23012 airfoil, standard C81 table 

# 

/bin/rm -rf naca23012 . out 
/bin/rm -rf naca23012.tab 
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unsetenv AFDECK1 
unsetenv AFTABLE 
setenv AFDECK1 naca23012 . c81 
setenv AFTABLE naca23012 . tab 

set campath=/ usr2/boyd/Cam/Sources/Camrad_mod2_rcs 
$campath/ airfoil > naca23012 .out «eoj 
NACA 23012 AIRFOIL (STANDARD C81 TABLE) 

&NLTABL 

NA = 1,16,28,88,100,115, 

NM = 1,7,21, 

M = 0. , .6, .95, 

OPREAD = 2, 

&END 

eoj 

exit 

5.4 Binary Input File Preparation 

#! /bin/csh 

# 

# set INPUTFILE environment variable, 

# compile blockdata file, 

# link files together, 

# run input preparation program ( input _prep) 

# 

# 

/bin/rm -rf bol05.out 
/bin/rm -rf bol05.bin 

# 

unset campath 
unset srcpath 

# 

unsetenv INPUTFILE 

setenv INPUTFILE bol05.bin 

set campath= l 7usr2/boyd/Cam/Sources/Camrad_mod2_rcs n 

# 

f 77 -c bol05_blockf ile . f 
f 77 $campath/blockf ile.o bo!05_blockf ile .o 
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$campath/FILEI . o -o input _prep 

# 

input _prep > bo 105. out 
exit 

5.5 Script Template 

This section describes a method to create a binary input file for CAM- 
RAD. Modi. First the user must create a file containing the necessary input 
BLOOKDATA code segment. Then, the binary input data file is created 
using the following example script: 

#!/bin/csh -v 

/bin/rm -f fort.* >& /dev/null 
unset case 
set case=runcase 
unset campath 

set campath=/usr2/boyd/Cam/Sources/Camrad_mod2_rcs 
unsetenv INPUTFILE 
unsetenv AFTABLE1 
setenv INPUTFILE bol05.bin 

setenv AFTABLE1 naca23012 . tab 

In -s ${runcase} . dat fort. 7 

In -s int_${runcase> .dat fort. 17 

In -s vind_${runcase} . dat fort. 18 

In -s wake_${runcase} . dat fort. 13 

In -s blade_${runcase} . dat fort. 14 

${ campath} /camrad.modl . 1 >& ${runcase> .out «eoj 
&NLCASE 

NFRS=-1 ,NFEIG=-1 ,NCASES=1 , 

&END 

&NLTRIM 

VEL = 0.15060, CTTRIM = 0.05610, 

RPM = 1041.000, APITCH = 4.250, 

&END 

&NLRTR &END 

&NLHHC &END 

&NLHHC2 &END 

&NLHIRES 
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OPINT = 1, OPWFCOR = 0, OPSEGD = 1, ITERINT=1 , 
MPSIINT=720 , JFIRST=1 , JLAST=720,MPSIWGP=72, 

COREINT=. 01815, KNWINT=180, OPNEGV =1, 

WMDLINT=2 ,0,0, 0,0, 2,2, 2, 2, 2, 3, 3, 3, 

WKMDL1 =2 ,0,0, 0,0, 2, 2, 2, 2, 2, 3,3, 3, 

MRAINT=75 , 

RAEINT= . 1500 , . 1783, . 1892, .2000, .2109, .2218, .2326,-2435, 
.2544, .2652, .2761 , . 2870 ,. 2979 ,. 3087 ,. 3196 , .3305, 
.3413, .3522, . 3631 ,. 3739 ,. 3848 , .3957, .4065, .4174, 
.4283,.4391, . 4500 ,. 4609 , .4717 , .4826, .4935, .5044, 
.5152, .5261, .5370, .5478, .5587, .5696, .5804, .5913, 
.6022, .6130, .6239, .6348, .6456, .6565, .6674, .6783, 
.6891 , .7000..7109, .7217, .7326, .7435, .7543,-7652, 
.7761, .7869, .7978, .8087, .8195, .8304, .8413, .8521, 
.8630, .8739, . 8848 ,. 8956 ,. 9065 ,. 9174 , .9282, .9391, 
. 9500 , . 9608 ,.9717,1. 0000 , 



0. ,0. 

o 

o 

o 

o 

,0. 

,0. 

,0. ,0. , 


0. ,0. 

o 

o 

o 

o 

,0. 

,0. 

,0. ,0. , 


0. ,0. 

o 

o 

o 




&END 

&NLBED 

OPBED = 

0, &END 




&NLSWP 


&END 




&NLWAKE 


&END 





&NLBURST OPBURST = 0, &END 
&NLROLL 


0PR0LLU=3 , OPLOWR= 1 , OPHIWR= 1 , 

CORELG= . 3 , ITERRUP=2 , ITERFRU=2 , 

ITERLGC=20 ,FLGC0RG=0 . 1 , 

NTIPFCT=1 ,TIPFC0=0 . ,TIPFC=0 .15915,0. ,0. ,0. ,0. ,0. , 

0 . , 0 . , 0 . , 0 . , 

NTC0R=9 ,TIPC0RE=0 .01 , 0.0166, 0.0233, 0.03, 

0.04, 0.05, 0.07, 0.1, 0.2, 0. , 
NSC0R=9 , SECC0RE=0 .01, 0.0166, 0.0233, 0.03, 

0.04, 0.05, 0.07, 0.1, 0.2, 0. , 
ISPIN=1 ,TAUC0=0 . ,TAUC1=1 . , 

IRUZC0R=0 , 0PR0LSS=0 , 

ICORYCB = 2, IFWLGC = 0, COREXP = 4, 

&END 

&NLMEAS IMODEIN = 0, I AERO IN = 0, &END 
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&NLLOAD &END 

eoj 
exit 

5.6 Input Parameters 


Following are the input parameters for CAMRAD.Modl. Parameters that 
are in italic print are new parameters and/or namelists not included in the 
original version of CAMRAD, but are now in CAMRAD.Modl. 
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4.5 NLCASE 


Namelist NLC'ASE 


JOB 

= 0 
#0 

integer parameter defining job type 
for new job (default) 
for old job or restart job 

RSWRT 

= 0 

integer parameter controlling restart file write 
to suppress write 

New Jobs Onlv: 
NC'ASES 


number of rases (default — 1) 

BLKDAT 

II A V 

integer parameter defining input source: 
read input file (default) 
use loaded blockdata and write input file 
use loaded blockdata 

RDFILE 

= 0 
# 0 

integer parameter controlling input file read: 
reau file for first case only 
read file for every case (default) 


OKI Jobs Onlv: 
START 

= 1 

integer parameter defining task: 
for trim restart (default) 


= o 

for flutter restart 


= 3 

for flight dynamics restart 


= 4 

for transient restart 


Note that the trim restart ran be followed by any or all of the other tasks (as defined by ANT\PE): for 
lint ter. flight dynamics, or trasient restart, only that task ran be done. 
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Additions 


NOISFL 


to Namelist NLCASE 


= 0 don't perform ROT ONE T/W OP WOP interface calculations 
(default) 

> 0 don't perform ROT ONE T/W OP WOP interface calculations 



4.6 NLTRIM 


Namelist NLTRIM 
OPREAD(IO) 


NPRNTI 

ANTYPE(3) 

TITLE (20) 
CODE 

OPUNIT 

NROTOR 

WBTTAB 


integer vector detining namelist read structure 
0 to suppress read of the namelist: 

components (new job only): 

(1) NLRTR, rotor 1 

(2) NLWAKE. rotor 1 

(3) NLRTR, rotor 2 

(4) NLWAKE. rotor 2 

(5) NLBODY 


tasks: 

(6) NLLOAD, rotor 1 

(7) NLLOAD. rotor 2 

(8) NLFLUT 

(9) NLSTAB 

(10) NLTR AN 

integer parameter controlling input data print 
0 use "short form'" print (no input data print) 

integer vector defining tasks for new job or 
trim restart 

0 to suppress: 

(1) flutter 

(2) flight dynamics 

(3) transient 

title for job and case (80 characters) 

alphanumeric code for job and case identification 
(4 characters) 

1 for English unit system (ft-slug-sec) 

2 for metric units system (m-kg-sec) 

number of rotors (1 or 2) 


0 default (use empirical models) 

1 use wing-body- tail t aides 
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DEBUG(25) 


integer vector controlling debug print: 
values: 

= 0 no debug print 
= 1 trace print 
= 2 low level print 
= 3 high level print 


elements: 

(1) time (sec) at which debug print enabled 

(2) input, 2-3 (INPTx) 

(3) initialization, 2 (INITC, INITRx, INITB, INITE) 

(4) trim iteration, 1-2 (TRIMI) 

(5) loads. 2 (LOADIx) 

(6) flutter matrices. 2-3 (FLUTM) 

(7) flutter coefficients, 2-3 (FLUTIx, FLUTAx) 

(8) flight dynamics, 2-3 (STABM, STABE) 

(9) transient. 2 (TRANI) 

(10) rotor/airframe motion and forces, 2-3 (RAMF) 

(11) blade modes. 2 (MODE, MODEx) 

(12) inertia coefficients, 2 (INRTCx) 

(13) airframe constants and matrices, 2 (BODYC, 

ENGNC, MOTNCx, BODYMx, ENGNMx) 

(14) induced velocity. 2 (WAKEUx, WAKENx) 

(15) rotor matrices. 2-3 (INTRMx) 

(16) hub /airframe motion and generalized forces. 2 
(MOTNHx. BODYVx, ENGNYx, MOTNFx, MOTNS) 

(17) rotor motion. 2-3 (MOTNRx) 

(18) rotor aerodynamics, 2-3 (AEROFx, AERBEDx ) 

(19) blade section aerodynamics, 3 (AEROSx, AEROSxB ) 

(20) body forces and aerodynamics. 2 (BODYF) 

(21) wake influence coefficients. 2 (WAKECx) 

(22) vortex line and sheet, 3 (VTXL, VTXS. VTXL2 ) 

(23) prescribed wake geometry, 2-3 (GEOMRx) 

(24) free wake geometry, 1-3 (GEOMFx) 

(25) timer, 1 (TIMER)* 


here the *‘x can take the values 1 or 2 depending on which rotor is being examined. 



YKTS 


airspeed (knots), V 
advance ratio, \/(HR) 


Y EL 

Note: 

Input either YEL or YKTS by namelist input. If neither parameter is defined 
in the input namelist. Y = 0. is used. 

YTIP rotor 1 tip speed, HR (ft/s er or 111 /sec) 

RPM rotor 1 rotational speed (rpm) 

Note: 

Input either YTIP or RPM by input namelist. If neither parameter is defined, 
the normal tip speed, YTIPN, is used from namelist NLRTR. Rotor 2 speed is 
calculated from the gear ratio, TRATIO. 

OPDENS integer parameter defining specification of aerodynamic 

environment: 

= 1 given altitude and standard day 
= 2 given altitude and temperature 
= 3 given density and temperature 

ALTMSL altitude above mean sea level (ft or m) 

(for OPDENS = 1 or 2) 

TEMP air temperature (deg F or deg C) 

(for OPDENS = 2 or 3) 

DENSE air density (slug/ ft A or kg //»**) 

(for OPDENS = 3) 

OPGRND integer parameter controlling ground effect analysis 

= 0 out of ground effect 
9 ^ 0 in of ground effect 

HAGL altitude helicopter center of gravity above ground 

for ground effect analysis (ft or m) 

OPENGN integer parameter specifying engine state: 

= 0 for normal operation 

= 1 for autorotation (engine inertia, engine damping, 
and throttle control torque zero: no engine speed 
degree of freedom) 

= 2 for engine out (engine damping and throttle control 
torque zero) 

AFLAP wing flap angle, t> F (deg) 

RTURN for free flight, trim turn rate, t/>- (deg/sec) 

(positive to right) 
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COLL 


initial values for controls (trimmed as appropriate 
collective stick displacement, or A0 yovv (deg), 
positive nose up 


LATCYC lateral cyclic stick displacement, b r (deg), 

positive left 

LNGCYC longitudinal cyclic stick displacement, (deg), 
positive aft 

PEDAL pedal displacement, (deg), 

positive right 

APITCH for free flight case: aircraft pitch angle. Oft (deg), 
jjositive nose uj) 

for wind tunnel case: rotor shaft angle, Op (deg), 
positive nose up 

AROLL for free flight case: aircraft roll angle, 0ft (deg), 
positive to right 

FT and o ft define orientation of body axes relative to 
earth axes) 

ACLIMB for free flight case: aircraft climb angle, #p/> (deg), 
positive up 

A\ AW for free flight case: aircraft yaw angle, i'fp (<leg), 

positive to right 

for wind tunnel case: test module yaw angle, r-p (deg), 
positive to right 

(Off and Iff define the orientation of the velocity axes relative to the earth axes: 
^ riimh = L sin (Off) and Y s ; f / f = V sin(rVr) cos(flpp)) 

MPSI number of steps per revolution in motion and loads 

analysis (max 36). For nonuniform inflow, must be a 
multiple of the number of blades 

MPSIR number of azimuth steps between update of airframe 

vibration and rotor matrices in harmonic motion solution 

MREA number of revolutions between tests for motion convergence 

in harmonic motion solution 

ITERM maximum number of motion iterations 

EPMOTN tolerance for motion convergence (deg) 

ITERC maximum number of circulation iterations 

EPCIRC tolerance for circulation convergence (ACV/rr) 
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DOF (54) 


rotor 1 
rotor 1 
rotor 1 
rotor 2 
rotor 2 
rotor 2 
airframe 
airframe 
drive train 
drive train 


DOFT(8) 


rotor 1 
rotor 2 


MHARM(2) 


MHARF(2) 


LEVEL(2) 

= 0 
= 1 


integer vector defining degrees of freedom used in 

vibratory motion solution. 

values: 

= 0 to omit 
= 1 to use 
order: 

bending: (max 10) q \ . <fr q i» 

torsion: (max 5 ) pu*p\ p 

gimbal : (max 1 ) 

bending: (max 10) q \ , </ 2 . ...,7m 

torsion: (max 5 ) pw.p i p: t 

gimbal : (max 1 ) 

rigid : yr . Z F 

elas t ic : ( m ax 10) q S7 q s , 

rotor/engine speed: r s < 
governer: Ai9, , A(9 </orri , A 0y OVV2 

integer vector defining l)lade bending degrees of freedom 

used for mean deflection (subset of DOF) 

values: 

= 0 to omit 
= 1 to use 
order: 

bending: (max 4) q \ , , q * , <74 

bending: (max 4) q\ . ry 2 . q:\ * q\ 

number of harmonics in rotor motion analysis (max 20) 

( = () for mean only) 

(1) rotor 1 

(2) rotor 2 

number of harmonics in airframe vibration analysis (max 10) 
(harmonics of N/rev) 

( — () for static elastic only) 

(suggest < MHARM/NBLADE, and the same value for both 
rotors if coupled hub vibration used (see OPHVIB) 

(1) rotor 1 

(2) rotor 2 

integer parameter specifying rotor wake analysis level 
(must be consistant with INFLOW) 
for unform inflow 

for nonuniform inflow with prescribed wake geometry 
for nonuniform inflow with free wake geometry 

(1) rotor 1 

(2) rotor 2 
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ITERU 

= 0 

number of uniform inflow wake- trim iterations 
at uniform inflow level, 
to skip 

ITERR 

= 0 

number of nonuniform inflow with prescribed wake- trim iterations 
at nonuniform inflow/prescribed wake geometry level, 
to skip 

ITERF 


number of nonuniform inflow with free wake-trim iterations 
at nonuniform inflow/free wake geometry level. 

NPRNTT 

< 0 
> 0 

integer parameter controlling trim /performance/load print, 
suppress print (last iteration always printed) 
print every NPRNTT-tli iteration 

NPRNTP 

< o 

integer parameter controlling performance print, 
suppress print 

NPRNTL 

< o 

integer parameter controlling loads print, 
suppress print 

MTRIM 


maximum number if iterations on controls to achieve trim 

MTRIMD 


number of iterations between update of trim derivative matrix 

DELTA 


control step in trim derivative matrix calculation 
(stick displacement, deg) 

FACTOR 


factor reducing control increment in order to improve 
trim convergence (typically .5) 

EPTRIM 


tolerance for trim convergence 

OPGOVT 

O * — i CO CO 

II II II II 

integer parameter specifying governor trim 

trim collective stick, 

trim rotor 1 governor 

trim rotor 2 governor 

trim both rotor governors 

trim targets for wind tunnel rases: 

CXTRIM 


Cs/cr 

XTRIM 


X/q (ft 2 or in 2 ) 

CTTRIM 


C't/v or C r / (t 

CPTRIM 


Cr/cr 

CYTRIM 


Cy/cr 

BCTRIM 


i- («log) 

BSTRIM 


■i» ( d(, g) 
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OPTRIM 


integer parameter specifying trim option 

free flight cases: 

0 no trim 

1 trim forces and moments with: 

, <%.. fij,, 0 p 7’. Of'T 

2 trim forces and moments with: 

Orr< i'fp 

3 trim forces, moments, and power with: 

<\s * M FT' Opr FT 

4 trim forces, moments, and power with: 

<V bf,'0FT' t'rr, $ FP 

5 trim symmetric forces and moments with: 

K. ft FT 

G trim symmetric forces, moments, and power with: 
<V 'ft FT' ft FP 
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wind tunnel rases: 


= 10 

no trim 


= 11 

trim Ct/v 

with (\) 

= 12 

trim C't/g 

with Or 

= 13 

trim Cpja 

with fyi 

= 14 

trim fi r < 

with 

= 15 

trim Cr/(T, ji v . 

with b i} .b r .b N 

= 1G 

trim Cr/iT.Cx/tT.Cy/tT 

with 

= 17 

trim CjJaX'x /(T'Cy/tr 

with b{\*b v Xr 

= 18 

trim CiJ(T, C\ /o\ 4 f ., 

with b {) ,b ( .,b M Xr 

= 19 

trim CiJ(T*X/q.Cy/(T 

with b\) % b r *b H 

= 20 

trim CjJa.X/qX'y/fJ 

with bi),b r Xr 

= 21 

trim CjJrr , X/q, 4 r , 

with , b r% <\ s « Oy 

— 22 

trim ;i r 

with b s 

= 23 

trim Ct/< t, /^ f . 

with bi).b s 

= 24 

trim CjJaX'x /a 

with <V<\ S 

= 25 

trim CjJcfX’x ja 

with b{).Hr 

= 26 

trim CfJaX'x /o rj1 r 

with b{),b s X r 

= 27 

trim Cfjcr, X/q 

with b^X* 

= 28 

trim CjJcf, X/q 

with b()Xr 

= 29 

trim Cj jo, X/q, 4 r 

with bi)*b s Xr 
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Additions to Namelist NLTRIM 


FACTM 0. < FACTM <1. lag factor to help convergence 

(default =1.) 


OPMXFWG =0 


= 1 
— 9 

= 3 
= 4 


use original in free wake 

geometry calculations (default) 

(during rollup calculations, this may he 

superceded by ICORYCB variable) 

use positiv e r, m(J . 

use negative r mrl , 

use outboard r mriJ . 

use inboard r mflJ . 


If NLCFD is to be used, the followmg must, be true: 
OPREAD(2)=2 
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4.7 NLTRTR 


Namelist NLTRTR 


TTTLE(20) 

title for rotor and wake data (80 characters) 

TYPE 

rotor identification (4 characters) 

(suggest MAIN, FRNT, or RGHT for rotor 1: 
TAIL. REAR, or LEFT for rotor 2) 

VTIPN 

normal tip speed, UBa (ft/sec or m/s ec) 

RADIUS 

blade radius, R (ft or m) 

SIGMA 

solidity, a = N * c,„ f -nn/(nR) 

GAMMA 

Lock number, - Ml = p {) ar )tlf . (Ul R 4 /I h 

(l)ased on standard density, a = 5.7, and mean chord) 

('. and a are only 

used to calculate the normalization parameters r inf . nn and I/ } ) 

NBLADE 

number of blades 

TDAMPO 

control system collective damping 
(ft -lb/ rad /sec or m-N/rad/s ec) 

TDAMPC 

control system cyclic damping 
(ft-lb/rad/sec or m-N/rad/sec) 

TDAMPR 

control system rotating damping 
(ft-lb/rad/sec or m-N/rad/sec) 

NUGC 

longitudinal gimbal natural frequency, vac* or 
teeter natural frequency, vr (per rev at normal 
tip speed VTIPN) 

NUGS 

lateral gimbal natural frequency, Vc;$ 
(per rev at normal tip speed VTIPN) 

GDAMPC 

longitudinal gimbal damping. Car or teeter 
damping Cf (ft-lb/rad/sec or m-N/rad/sec) 

GDAMPS 

lateral gimbal damping, CVv.s 
(ft-lb/rad/sec or m-N/rad/sec) 

LDAMPC 

linear lag damper coefficient, C^: 
estimated damping if a nonlinear damper is used 
(LDAMPM > 0.): the lag mode has structural damping 
also (GSB) 

(ft-lb/rad/sec or m-N/rad/sec) 

LDAMPM 

maximum moment of nonlinear lag damper, Mf.p: 
(ft-lb or m-N) 

= 0 to use linear lag damper 
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LDAMPR 


lag velocity, Qj.o* where maximum moment of lag 
damper occurs (rad/sec): hydraulic damping below 
( i p and friction damping above. 

GSB(NBM) bending mode structural damping, (j s 

GST(NTM) torsion mode structural damping, g s 

where NBM is the number of bending modes and NTM is the number of torsion modes 


as used in 

DOF. 


ROTATE 

= 1 

integer parameter specifying rotor rotation direction 

as viewed from above. 

counter-clockwise. 

clockwise. 

OPHVIB 


integer parameter specifying hub vibration contributions: 
gravity and static velocity terms always retained. 

(=0 to suppress) 

(1) vibration due to this rotor 

(2) vibration due to other rotor (must = 0 if ih/M\ ^ 1) 

(3) static elastic motion 

BTIP 


tip loss factor, B 

OPTIP 

= 1 
— o 

integer parameter specifying tip loss type: 
tip loss factor 
Brand tl function 

LINTW 

= 0 
#0 

integer parameter specifying twist type: 
nonlinear twist 
linear twist 

TWISTL 


linear twist rate, (deg). 

used to calculate TWISTA and TWISTI if LINTW = 1 

OPITSLD 

= 0 
=: 1 
— 2 

integer parameter controlling use of unsteady lift, 
moment, and circulation terms: 
suppress terms 
include terms 

zero for stall (15 deg < |a| < 165 deg) 

(OPUSLD is set to zero internally if OPBED = 1) 

OPCOMP 

= 0 

integer parameter controlling aerodynamic model: 
incompressible loads 
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INFL0W(6) 


RROOT 

RGMAX 

MRA 

RAE(MRA+1) 

CHORD(MRA) 

XA(MRA) 

THETZL(MRA) 

TWISTA(MRA) 
X AC (MRA) 


Inflow Model 

integer parameter defining induced velocity calculation: 

(must he consist ant with LEVEL) 

(1) at this rotor: 0 for uniform, 1 for nonuniform 

(2) at other rotor: 0 for zero, 1 for empirical, 

2 for average at huh, 3 for nonuniform (only if il-ifil \ — 1) 

(3) at wing-bodv: 0 for zero, 1 for empirical, 

2 for nonuniform 

(4) at horizontal tail: 0 for zero, 1 for empirical, 

2 for nonuniform 

(5) at vertical tail: 0 for zero, 1 for empirical, 

2 for nonuuiform 

(6) at point off rotor disk: 0 for zero. 1 for nonuniform 


root vortex position for wake model, r rooi /R 


rr /R (induced velocity calculation using maximum 

hound circulation outboard of /T,„„ r /R ). 


Blade section aerodynamic characteristics 
number of radial segments (max 30) 

radial stations r/R at edges of aerodynamic segments: 
sequential from root to tip. 

The following quantities are specified at the 
midpoint of the aerodynamic segments: 


blade chord, c/R 


offset of aerodynamic center aft of the elastic axis, .r \j R. 
.r i is the point about which the moment data in the 
airfoil tables is given. 


incremental pitch of zero lift line. Oyj (deg): 
can he included in TWISTA; i 9 Zf is the pitch 
of the axis corresponding to the zero angle of attack 
in the airfoil tables, relative to the twist angle (TWISTA) 


blade twist relative ,75R, 0 tlv (deg) 


offset of the aerodynamic center (for unsteady aero- 
dynamics) aft of tlie elastic axis, .r U -/R 
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MCORRL(MRA) 

MC’ORRD(MRA) 

MCORRM(MRA) 

Mach number correction factor f\j = for lift 

Mach number correction factor f\j = for drag 

Mac'll number correction factor f\i = for moment 

Blade section inertial and structural characteristics 

MRI 

number of radial stations where characteristics defined 
(max 51 ) 

RI(MRI) 

radial stations r/R: sequential from root to tip. 
( R I ( 1 ) = 0. and RI(MRI) = 1.) 

MASS(MRI) 

section mass, m ( slug/ft or kg/m) 

EIXX(MRI) 

chordwise bending stiffness (lb -ft 2 or N -n? 2 ) 

EIZZ(MRI) 

fiapwise bending stiffness (lb -ft 2 or N -in 2 ) 

XI(MRI) 

offset of center of gravity aft of elastic axis, .r j/R 

XC(MRI) 

offset of tension center aft of elastic* axis, .rr/R 
(at the tip, XC should be set nearly equal to XI) 

KP2(MRI) 

polar radius of gyration about elastic axis, k p 2 / R 1 

ITHETA(MRI) 

section moment of inertia about elastic axis, h 
(slug-ft or kg-m) 

G.l(MRI) 

torsional stiffness, GJ (l!>-/t 2 or N-/n 2 ) 

TWISTI(MRI) 

blade twist relative .75R, (deg) 
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OPSTLL 

Stall Model 

integer parameter defining stall model: 

= 0 no stall 
= 1 static stall 
= 2 McCroskev stall 

= 3 McCroskev stall with dynamic stall vortex loads 
= 4 Boeing stall 

= 5 Boeing stall with dynamic stall vortex loads 
(OPSTLL set = 0 internally if OPBED = 1) 

(the stall delay 

r can 

be suppressed by setting TAU = 0.) 

OPYAW 

= 0 
= 1 
— 2 
= 3 

integer parameter defining vaw flow corrections: 

both yawed flow and radial drag included 

no yawed flow (cos(A) — 0.) 

no radial drag (F r = 0.) 

neither yawed flow not radial drag included 

TAU(3) 


stall delay time constants for lift, drag, and moment: 
Tr.*Tp*TM (calculated of < 0) 

ADELAY 


maximum angle of attack increment due to stall delay 
<' ,n(u<if i<uj (deg) 

AMAXNS 


angle of attack in linear range for no stall model 
(deg) 

PSID(3) 


dynamic stall vortex load rise and fall time (azimuth 
increment) for lift, drag, and moment: (deg) 

ALFDS(3) 


dynamic stall angle of attac k for lift, drag, and moment: 
<'<ts (deg) 

ALFRE(3) 


stall recovery angle of attack for lift, drag, and moment: 
(deg) 

CLDSP 


maximum peak dynamic stall vortex induced lift coefficient: 
AC/ t/ , 

CDDSP 


maximum peak dynamic stall vortex induced drag coefficient: 

CMDSP 


maximum peak dynamic stall vortex induced moment coefficient 

ni th . 

KHLMDA 


factor, h f , , for hover induced velocity (typically 1.1) 

KFLMDA 


factor, Ky, for forward flight induced velocity 
(typically 1.2) 

FXLMDA 


factor, fj . , for linear inflow variation in forward flight 
(typically 1.5) 

FYLMDA 


factor, f !r for linear inflow variation in forward flight 
(typically 1.) 
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FMLMDA 


factor, for linear inflow variation due to hub moment 
(typically 1.) 


FACTWU factor introducing lag in Cj<C\i r , and C\i y used to 

calculate induced velocity (typically .5) 

KINTH factor for hover interference velocity at other rotor 

(k -21 or K| 2 ) 

KINTF factor for forward flight interference velocity at 

other rotor (k^i or />* 1 * 2 ) 

(linear variation between KINTH at // — .05 and KINF at // ■ = .1 is used) 

KINTWB factor for rotor-induced interference velocity at 

wing-body. A ir 

KINTHT factor for rotor-induced interference velocity at 

horizontal tail. Kh 

KINTYT factor for rotor-induced interference velocity at 

vertical, K\ 

(A'vr, A'ff, and I\\ equal fraction of fully-developed wake times maximum 
fraction surface in wake) 

HINGE integer parameter specifying blade mode type 

= () hinged 

= 1 cantilever 

= 2 art inflated (flap and lag modes only) 

NCOLB number of collocation functions for bending mode calculations 

(total number of flap and lag, alternating): (max 20) 

NCOLT number of collocation functions for torsion mode calculations 

(max 10) 

NONROT integer parameter. 

/ 0 to calculate nonrotating bending frequencies 

EPMODE criterion on change of collective pitch to update 

blade bending modes, Atfy.-, (deg) 

MASST tip mass (slug or kg); the tip mass can also be 

directly in the section mass distrubution 

XIT offset of tip mass center of gravity aft of elastic 

axis, .r//R 

MBLADE l>lade mass (slug or kg). 

if < 0 integral of section mass used (with mass 
included at r=0. to account for the hub mass) 
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EFLAP 

ELAG 

KFLAP 

KFLAG 

RCPLS 

TSPRNG 


flap hinge offset, f’y/R. (extent of rigid hub for cantilever blade) 
lag hinge offset, eif R. (extent of rigid hub for cantilever blade) 

flap hinge spring (ft -lb/rad or m-N/rad) 

lag hinge spring (ft-lb/rad or m-N/rad) 

hinge spring parameter. 'R., 

hinge spring parameter, ff M „ 


(hinge spring pitch angle is = 0 H „ + ) 


RCTL 


structural coupling parameter, R, 

(effective pitch angle, R6 1 , used to calculate 
blade bending modes (normally R = 1.)) 

NOPB 


integer parameter specifying twist inboard of rp \ 
(= 1 for no pitch bearing) 

WTIN 

— ' <M 

II II 

integer parameter defining control system stiffness input: 
for ho 
for uSf) 

FTO 


control system frequency, ^ 

(per rev, at liotmal tip speed, YTIPN): 
collective 

FTC 1 


cyclic 

FTR 


react ion less 

KTO 


control system stiffness. A# 
(ft-lb/rad or m-N/rad): 
collective 

KTC 


cyclic 

KTR 


react ion less 

KPIN 

= 1 
— 2 

integer parameter defining pitch /bending coupling input: 

for input 

for calculated 

(negative to suppress cosine terms factors in I\p t and I\ r 



root geometry to calculate nitch /bending colliding 



(KPIN = 2 or -2) 1 b 

PHIPH 


pitch horn cant angle, (deg) 

PHIPL 


pitch link cant angle, 0 p / (deg) 
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RPB 

pitch bearing radial location, rppj R 

RPH 

]>itch horn radial location, rppj R 

XPH 

pitch horn length, xpp j R 

ATANKP ( NBM ) 

pi tch /bending coupling tan" 1 (Kp t ) (deg), 
for pitch horn level (KPIN = 1 or -1) 

DEL3G 

pitch/gimbal coupling tan" 1 ( A> 0 . ) (deg), 
for pitch horn level. 

RFA 

feathering axis radial location, rp ,\/ R 

ZFA 

gimbal undersling, zp\f R 

XFA 

torque offset, .rp\/R 

CONE 

precone angle, (deg) (positive up) 

DROOP 

droop angle, (deg), ^t:> = 9, 

(positive down from precone) 

SWEEP 

sweep angle, dp (deg), at 6-- % = 0, 
(positive aft) 

FDROOP 

feathering axis droop angle, <Va, (deg), 
(positive down from precone) 

FSWEEP 

feathering axis sweep angle, (deg), 

(positive aft) 

MRM 

number of radial stations for integration of inertial 
coefficients (max 50) 

MRB 

number of radial stations in blade mode calculation 
(max 50) 
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4*8 NLHHC 


New Namelist NLHHC 


PSICOLL 

(*•</) 

l'r 

(default 

= 0. 

PSILAT 

(<if' a) 

^ /<»/ 

(default 

= 0. 

P SI LON 

(deg) 

^ itHi 

(default 

= 0. 

THCOLL 

(deg) 

9 r 

(default 

= 0 . 

THLAT 

(deg) 

a, 

(default 

= 0 . 

THLON 

(deg) 

Mfon 

(default 

= 0. 


Variables in this namelist are used only if OPHHC = 1 in namlist NLHHC2. 



4.9 NLHHC2 


New Namelist NLHHC2 


OPHHC 

= 0 

no HHC used ( default ) 


= 1 

use HHC input from NLHHC 


-- 2 

use HHC input from NLHHC 2 

NHHC 


(integer) number of HHC harmonics 
(max 12) 

HHCPO 

(<leg) 

(real) HHC “collective" 
(default = 0.) 

AHHC(NHHC) 

(deg) 

(real) vector of HHC cosine amplitudes 
( default = 12*0.) 

DHHC(NHHC) 

(deg) 

(real) vector of HHC sine amplitudes 
(default = 12*0.) 
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4.10 NLHIRES 


New Namelist NLHIRES 
OPINT 

ITERINT 

MPSIINT 

J FIRST 

JLAST 

MPSIWGP 

FACTINT 

DLSINT 

MR A IN T 

RAEINT(MRAINT-hl) 

ITERNW 

FACTNW 

KNWINT 

OPNEGV 

OPCSNW 


— 0 no HIRES performed (default) 

— 1 HIRES calculations performed 

number of fa, r wake iterations (default — 1) 

number of azimuth steps (max. 720) 

index of first azimuth step ( default = 1) 

index of last azimuth step 

number of azimuths for wake geometry print 

lag factor for ITERINT iterations (default — 1) 

lifting surface correction (default = -1.) 

(see DLS in NLRTR for more detail.) 

number of HIRES radial stations (max 100) 

edges of aerodynamic segments ( re R) 

number of HIRES near wake iterations 
(default = 0) 

lag factor for ITERNW iterations (default = .5) 

extent of HIRES near wake. Near wake extends 
from 0 = 0 to o =KNWINT*\& 
behind reference blade. 

negative tip vortex and, inboard, wake 
= 0 include the tip vortex and inboard, wake in 
the "near ware " . 

= 1 eliminate the Up vortex and inboard wake 
that is in the "near wake", (default) 

= 0 distributed core model in shed NW (default) 

— 1 concentrated eore model in shed NW 
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OPCTNW 

= 0 
= 1 

distributed core model in trailed NW (default) 
concentrated core model in trailed, NW 

MDLSNW 

= 0 
= 1 
= 2 

shed near wake model 

omit shed near wake 

line vortex w/ stepped circulation 

line vortex w/ linear circulation (default) 

MDLTNW 

= 0 
= 1 
= 2 

trailed near wake model 

omit trailed near wake 

line vortex w / stepped circulation 

line vortex w/ linear circulation (default) 

COR ESN W 

< 0 . 

near wake shed core size (re. R) 
use default = .5* (panel length) 
(default — .015) 

CORETNW 

< 0 . 

near wake trailed core size (re. R) 
use default =.5* (panel width) 
(default = .009) 

OPNWMIN 

= 0 
= 1 

minimization of near wake panels option 
supress minimization ( default) 
perform reduction in number of panels 

OPNWC'RC 

= 0 
= 1 

updating of near wake 

update at each azimuth (default) 

update each revolution 

OPDCORE 

= 0 
= 1 

use single core model (default) 
use dual core model 

CO REIN T 


corf size for single core HIRES model 

OPT VC OR 

= 0 
= 1 
= 2 
= 11 
= 12 

single core expansion model 

constant size core ( default) 

step fn. in far wake core size (HIRES only) 

core fn. in far wake core size (HIRES only) 

same as = 1 except also in trim loop 

same as — 2 except also in tri,m. loop 

PHIINC 


wake age at which to increase core size 
to RCORINC by a step function 

RCORINC 


core size after age PHIINC 

RCORO to 
R COR 10 


coefficients of core size expansion function 
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OP COR AC 

= 0 

use constants RCORDCA and RCORDCB for inner 
and, outer core sizes of dual core 


= 1 

use DCORACl and DCORAC2 to calculate inner 
core size, use RCORDCB for outer core size in 
dual core 

RCORDCA 


inner core size (re. R) 

RCORDCB 


outer core size (re. R) 

DC0RAC1 


used, in calculating inner core size 

DC0RAC2 


used, in calculating inner core size 

NDCCODA 


number of terms in dual core expansion (inner) 

NDCCODB 


number of terms in dual core expansion (outer) 

DCCORFAO 


term in inner core expansion function 

DCCORFA (NDCCODA) 


term in inner core expansion function 

DCCORFBO 


term in outer core expansion function 

DCCORFBf NDCCODB ) 


term in outer core expansion function 

GAMACST 


constant used in calculating inner core 
strength 

DR PR GOT 


aero. coll, point offset at root (re R) 
(default = 0.) 

DRPTIP 


aero. coll, point offset at tip (re R) 
(default = 0.) 

DRP 


aero. coll, point offset on blade (re R) 
(default = 100 * 0J 

OPSEGD 

= 0 

segmentation of vortex ( dements 
Off 


= 1 

on (default) 

OPWFCOR 


read, and use TF correction input files 


= 0 

off ( default) 


= 1 

on 

WKMDLl(lS) 


same as WKM()DL(1S) in NLRTR. This used for far 
iterations 1 thru ITERINT-1 

WMDLINT(IS) 


same as WKMODL(IS) in NLRTR. This used, for 
last far wake iteration. 

NLINT 


(do not use) 

NGINT 


(do not use) 
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MRGINT(NGINT) (do not. me) 

MRLINT(NLINT) (no not me) 
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4.11 NLBED 


New Namelist NLP ED 


OPBED 

= 0 


= 1 

HCOR 

<0. 


> 0 . 

ICURV 

= 1 


= 0 

ILEED 

= 1 


= 0 


(integer) use Johnson aerodynamics (original C AMR AD) 
(integer) use Indicia \l aerodynamics 

(real) core size for trailed near wake (re. R) 
use default core size (1/2 )Ar (default) 
use this constant core size 

(integer) use circular arc trailed near wake (default) 
(integer) use straight trailed, near wake. 

(integer) "lead" terms in trailed wake calculation 
( default) 

(integer) do not use "lead" terras (not recommended) 
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4.12 NLSWP 


Namelist 


SWPLO(SO) 


NLSWP 


(deg) (real) sweep of quarter chord line for use in low resolution 
aerodynamic corrections at aerodynamic collocation 
points 

(default = 30*0.) 


202 



4.13 NLWAKE 


Namelist NLWAIv'E 


FACTWN 

factor introducing lag in hound circulation used to 
calculate induced velocity. 

OPYXYY 

integer parameter: 

(= 0 to suppress x and y components of induced 
velocity calculated at the rotors) 

KNW 

extent of near wake, A \ u 

KRW 

extent of rolling up wake, I\ru 

KFW 

extent of far wake and tip vortices, AVu 

KDW 

extent of far wake and tip vortices for points off 
rotor disk. I\r\\ 


(wake age o =K2u’; all K > 1) 

RRU 

initial radial station of wake rollup, i'r { /R 

FRU 

initial tip vortex fraction of r mfU . for rollup, / m - 

PRU 

extent of rollup in wake age, Ori (deg) 

FNW 

tip vortex fraction of I\\/ for near wake. f \\\ 

DYS 

sheet edge test parameter. d r „ 
(< 0 to suppress test) 

DLS 

lifting surface correction parameter. df s 
(< 0 to suppress correction) 

CORE(5) 

vortex core radii, r r j R 

(1) tip vortices 

(2) hurst tip vortices 

(3) tip vortices in far wake off rotor 

(4) trailed lines (< 0 for default = s/2) 

(5) shed lines (< 0 for default = t/2) 

OPCORE(2) 

integer parameter specifying vortex core type: 
values: 

— 0 for distributed vorticity 
= 1 for concentrated vorticity 

elements: 

(1) tip vortices 

(2) inboard wake 
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0PNWS(2) 


LHW 
OPHW 
OP RTS 


integer parameter controlling action when inflow and 
circulation points conincide in near wake (o — 0) 
and sheets are being used: 
values: 

= 0 to use two sheets 
= 1 to use two lines 
= 2 to use single sheet 
elements: 

(1) shed wake 

(2) trailed wake 

number of spirals of far wake for axisymmetric case, Lfjvv 

integer parameter: (= 0 for axisymmetric wake geometry) 

integer parameter: 0 to include rotation matrices 

(i?/\s-,etc.) in influence coefficients) 


WKMODLf 13) 


integer parameter defining wake model: 
values: 

= 0 to omit element 

= 1 for line segment with stepped circulation distribution 
= 2 for line segment with linear circulation distribution 
= 3 for vortex sheet element 
elements: 

(1) tip vortices (stepped line or linear line) 

(2) near wake shed vorticity 

(3) near wake trailed vorticity 

(4) rolling up wake shed vorticity 

(5) rolling up wake trailed vorticity 

(6) far wake shed vorticity 

(7) far wake trailed vorticity 

(8) far wake (off rotor) shed vorticity 

(9) far wake (off rotor) trailed vorticity 

(10) bound vortices (no sheet model) 

(11) axisymmetric wake axial vorticity (no line model) 

(12) axisymmetric wake shed vorticity (no line model) 

(13) axisymmetric wake ring vorticity (no line model) 


MRG 

NG(MRG) 

MRL 

NL(MRL) 


number of circulation points for near wake (<MRA) 

circulation points, identified by aerodynamic segment 
number: no, for i = 1 to MRG (corresponding r, must 
be between r roo/ / R and 1.) 

number of inflow points (<MRA) 

points at which the induced velocity is calculated, 
identified by aerodynamic segment number: 
n / ; for i = 1 to MflL 
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OPWKBP(3) 


VE LB 

DPHIB 

DBV 


QDEBUG 


integer parameter controlling blade position for 
wake analysis: 

( 1 ) =() to suppress inplane motion 

( 2 ) =0 to suppress all harmonics except mean 

(3) = 0 for linear from r = r root to r = 1 . 


core burst propogation rate, \ ) } = 06/ Or 

core burst age increment, A&t, (deg) 

core burst test parameter, d hv 
(< 0 to suppress bursting) 


velocity criterion for debug print: print if 


r-fc/r 


> QDEBUG 


KRWG 


OPRWG 


Prescribed Wake Geometry 

extent of prescribed wake geometry, Kffwa 

(age o = A R\y (f Ar) 

(max 144) 

integer parameter defining prescribed wake geometry model: 
= 1 from A i = f\ A, R > = f 2 A, input A 3 , input K A 

= 2 option # 1 without interference velocity in A 

— 3 from A 1 . I\ > , A' :t , A 4 

Langrebe prescribed wake geometry: 

= 4 from Cr 
= -3 from 

— G from A 

= 7 from A without interference 

Kocurek and Tangier prescribed wake geometry: 

= 8 from Ct 

— 9 from t max 
= 10 from A 

= 11 from A without interference 


FWGT(2) 

FWGSI(2) 

FWGSO(2) 


KWGT(2) 

KWGSI( 2 ) 

KWGSO(2) 


factors / 1 and fi for jirescribed wake geometry: 
tip vortex 
inside sheet edge 
outside sheet edge 

constants A [ , A 2 , A , A 4 for prescribed wake geometry: 
tip vortex 
inside sheet edge 
outside sheet edge 
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KFWG 

OPFWG 


ITERWG 


FACTWG 


RTWG(2) 


Free Wake Geometry 

extent of free wake geometry distortion calculation, AViv r 
(age o — KFWG * Ar): suggest (.4///)MPSI; 

(maximum 144 ) 

integer parameter defining free wake geometry model: 

— 1 Scully free wake geometry 
= 2 same as = 1, without interference velocity 

number of wake geometry iterations (suggest 2 or 3) 
(suggest 4 or 5 for low speed cases) 

factor introducing lag in distortion calculation to 
improve convergence (suggest .5) 

radial station r/R of trailed vorticity: 

(1) inside sheet edge 

(2) outside sheet edge, or trailed line (suggest .4) 


WGM()DL(2) 


COREWG(4) 


MRYBWG 

LDMWG 

NDMWG(MPSI) 


integer parameter defining wake model: 
values: 

= 0 to omit 
= 1 for line segment 
= 2 for sheet element 
elements: 

(1) inboard trailed wake elements 
(1) shed wake elements 

vortex core radii, r r j R 

(1) tip vortices 

(2) burst tip vortices (< 0 for default — unburst value) 

(3) inboard trailed lines (< 0 for default = 
.5*(RTWG(2) - RTWG(l))) 

(4) shed lines (< 0 for default = 0.4 At’) 

number of wake revolutions used below point where 
induced velocity being calculated (suggest 2) 

integer parameter, /da/: general update every 
lr>M At/ 1 increment in boundary age 
(suggest 180 deg/ At’) 

integer parameter, np\f(i'j): 
boundary update every n da/ increment in age. 
function of i'j = j At’, j = 1 to MPSI: 
suggest 90 deg/ At’ for and aft: 
and 45 deg/ A r on shies. 
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DQWG(2) incremental velocity criteria: 

(suggest 0.04A| to().08A|) 

(1) near wake elements defined by |Ar/| >DQWG(1) 

(2) integrate bound vortex line in time over if 

| Ar^ >DQWG(2) 

IPWGDB(2) integer parameters controlling debug level 3 print 
of wake geometry distortion 

(1) IPR: print distortion before general update every 
IPR* Ac : ( = 0 to suppress) 

(1) IPRS: print distortion after each iteration every 
IPRS*Ac ; ( = 0 to suppress: 
last iteration printed in full) 

QW'GDB parameter controlling debug level 3 print: 

induced velocity contribution of wake element printed 
if | Aryl > QWGDB: 

(suggest 0.5A ( to 1.0A | ) 
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4*14 NLBURST 


New Namelist 

OPBURST 

PSITOL 

ZTOL 

CORMULT 

CIRMULT 


NLBURST 


tip vortex bursting model 
= 0 no bursting calculations (default) 
= 1 bursting calculations performed. 

— i ’ tolerance used, to define a. 

"close" interaction (radians). 

< 0 to allow only BVI crossings to 
cause a bursting, (default) 

height tolerance for bursting 
criteria, (re R) 

core size multiplication factor 
applied at a burst event. 

CrMursI = r r%oM *CORMULT 

circulation multiplication factor 
applied at a burst event. 

r (Mu rsi = T rM * CIRMULT 

(not fully functional yet) 
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4.15 NLROLL 


New Name list NLROLL 


OP ROLL U 
OP LOW R 

OPHIWR 
COREL G 
ITER R UP 
ITERFRU 

ITER LG C 

FLGCORG 

NTCOR 

TIPCORE(NTCOR) 
NS COR 

SECCORE(NTCOR) 

NTIPFCT 

TIPFCO 

TIPFC( NTIPFCT) 


— 0 no rollup calculation performed 

— 1 rollup calculations in low resolution only 

— 2 rollup calculations in high resolution only 

— 3 rollup calculations in both low and high, resolution 

= 0 do not print low resolution rollup data. 

— 1 write low resolution rollup data. ( circulation 

distribution , tip and secondary vortex positions. 

tip and. secondary vortex radial circulation distributions) 

— 2 same as = 1 plus large core circulation data. 

= 0 do not print hi resolution rollup data 

— 1 write hi resolution rollup data 


tip core size for the large core calculations 
(ref. radius) typically = 0.3 


number of iterations for convergence between large core 
calculations the wake geometry -trim iteration 

number of iterations between the wake geometry and. tin 
trim iteration when using the rollup results (corresponds 
to ITERF in the trim iteration without rollup ). 


number of iterations between the induced velocity 
calculation and the large core circulation calculation 
(typically = 20) 

lag factor for the large core calculation iteration 
(typically = 0.1, similar to FACTWN in NLWAKE) 

number of tip vortex cores for rollup calculation 
(max 10) 


tip core sizes (ref. radius) 

number of secondary vortex cores for rollup calculation 
(mux 10) 

secondary core sizes (ref. radius) 

number of coefficients in function to phase-in rollup 
calculation results, (max 10) 

coefficient in function to phase -in rollup results 

coefficients for function tip phase-in rollup results 
(for linear phase-in. TIPFC(l) - 1.0/ (age in radians to 
complete phase -in process)) 
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ISPIN 

TA UCO 
TAUC1 
ISECPH 

iruz cor 

NSEARCH 

IRUDZ 

OPROLSS 

CROLSSXY 

CROLSSZ 

ICORYCB 

IFWLGC 

COREXP 


() do not include spinning relation between tip and secondary 
vortex cores. 

1 include spinning relation between tip and secondary 
vortex cores . 

coefficient to phase-in tip /secondary vortex spinning 
coefficient to phase-in tip /secondary vortex spinning 

0 no s eco n d ary p h as e - in ( d efa u 1 1) 

1 secondary phase-in (don’t use) 

z-correction in rollup calculations 

0 don't use (default: leave blank) 

1 use (do not use this . it's not functioning) 

( leave blank) 

(leave blank) 

phase-in model choice 

0 use phase- in function above (default) 

1 use Spreiter- Sachs model (do not use. it's not functioning) 
( leave blank ) 

(leave blank) 

multi- core model choice 

0 use original multi- core model (default) 

1 single core variable, multi- core model 

2 multi core variable, multi- core model 

usage of large core circulation in free wake geometry 

0 use. original circulation (default) 

1 use single weighted large core max circulation 

2 use double weighted, large core max circulation 

core exponent for multi-core model (integer) 

(default = 1) (Scully- type core) 
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4.16 NLCFD 


New Namelist NLCFD 


OPCFD 


OPBVI 


PHICFD 

RDB(G) 


BDB(6) 


OPMOTN 


= 0 do not perform the partial angle calculations 
= 1 compute partial angles (OPMOTN— 0) or partial 
inflow (OPMOTN = l) and output inforrnation to 
the ALPHAP.DAT file 

= 2 same as = 1 except' files CFDAERO.DAT and CAMAERO.DAT 
are read at start of run and used in the aerodynamic 
calculations. 

= 0 if OPCFD= L do not remove tip vortex elements inside 
the BVI box from partial angle or inflow 

— 1 if OPCFD— 1 . remove tip vortex elements inside the 

B VI box from partial angle or inflow, and append a 
table of tip vortex trajectory and strength to the 
file ALPHAP.DAT ' 

( deg) age wise extent of CFD box behind reference blade . 
used, to limit vortex element testing for mclusin 
in the CFD box. 

CFD box dimensions relative to the blade tip (re. R) 

(all are positive numbers ) 

(1) leading edge face 

(2) radial distance outboard of tip (outside edge face ) 

(3) trailing edge face 

(4) radial station of inboard face 

(5) height, of upper face 

(6) height of lower face 

BVI box dimensions relative to hub (re. R) 

( all are positive numbers and. in shaft, '• coordinates ) 

(1) upstream face 

(2) starboard face 

(3) downstream face 

(4) port, face 

( 5) upper face 

(6) lower face 

— 0 if OPCFD — 1. output a partial angle table to the file 

ALPHAP.DAT. and if uPBVI— 1. transform the tip vortex 
segment endpoint coordinates to the M flapped blade “ 
coordinate system before appending the wake table to 
ALPHAP.DAT 

= 1 if OPCFD= 1. output vertical velocity at the blade 

section to ALPHAP.DAT. and if OPBVI— 1. append a wake 
table with untransformed wake coordinates to ALPHAP.DAT 
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4.17 NLMEAS 


New Namelist NLMEAS 


IMODEIN 

= 0 

use regular blade motion calculation 


= 1 

use input blade motion 

IAEROIN 

= 0 

use regular aero calculation 


= 1 

use input aero information 

FLAP FILE 


name of flap file (character* 80) 

LAG FILE 


name of lag file ( character *80 ) 

PITCHFILE 


name of pitch file (character *80) 

CNFILE 


name of cn file (character *80) 

NR A DIN 


number of radials in FL AFFILE. LAG FILE. PITCH FILE 
(max — 50) 

RIN(NRADIN) 


radials in FLAPFILE.LAGFILE. PITCHFILE 

MPSIIN 


number of azimuths FLAPFILE.LAGFILE. PITCHFILE 
(max — 2049) 

NRADINA 


number of radials in CNFILE 
(max — 50) 

RINA(NRADINA) 


radials in CNFILE 

MPSIIN A 


number of azimuths CNFILE 
(max = 2049) 

FLAPO 


fl.ap angle for hinged case (deg) 
(default = 0.) 

LAGO 


lag angle for hinged case (deg) 
(default — 0.) 
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4.18 NLBODY 


Namelist NLBODY 


TITLE (20) 
WEIGHT 


IXX 

IVY 

IZZ 

IXY 

IXZ 

IYZ 

TRATIO 


CONFIG 


ASHAFT(2) 


ACANT(2) 


ATILT 


HMAST 

DPSI21 


CANTHT 

CANTVT 


title for airframe and drive train data (80 characters) 

aircraft gross weight including rotors (lb or kg) 

aircraft moments of inertia including rotors 
(shift- ft_ or k 

t X X 

hfu 

n], 

h, 

hr- 

ratio of rotor 2 rotational speed to rotor 1 
rotational speed, il-i/il \ 

(transmission gear ratio r/,/r/ 2 ) 

integer parameter specifying helicopter eoniguration: 

= 0 for one rotor 

= 1 for single main rotor and tail rotor (rotor 2 is the tail rotor) 
= 2 for tandem main rotors (rotor 2 is the rear rotor) 

= 3 for tilting proprotor aircraft (rotor 2 is the left rotor) 

shaft angle. Op (deg), positive read ward: 

(1) rotor 1 

(2) rotor 2 

shaft cant angle, op (deg), positive to right 
for main rotor: positive upward for tail rotor: 
positive inward in helicopter mode for tilt rotor: 

(1) rotor 1 

(2) rotor 2 

nacelle tilt angle, n /> (deg) for tilting proprotor 
configuration only: (=0. for airplane mode, 

=90. for helicopter mode) 

rotor mast length from pivot to hub (ft or m), 
for tilting proprotor only 

Ac' 2 i (deg); rotor 2 azimuth angle , t 2 
when rotor 1 azimuth angle -Ci = 0.: 
must be 0. if 1 / 1. 

horizontal tail cant angle, Qpr (deg), positive to left 
vertical tail cant angle. & \ t positive to right 
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FSCG 

BLCG 

WLCG 


FSR1 

BLR1 

WLR1 

FSR2 

BLR2 

WLR2 

FSWB 

BLWB 

WLWB 

FSHT 

BLHT 

WLHT 

FSYT 

BL\T 

WLYT 

FSOFF 

BLOFF 

WLOFF 

CONTRLZ(ll) 


location (fuselage station, butt line, and waterline) of 
aircraft components relative to a body fixed reference 
system having an arbitrary orientation and origin; fuselage 
station (FS) positive aft, butt line (BL) positive to right, 
and waterline (\YL) positive up (ft or m) 

aircraft center of gravity location 


rotor 1 hub location (right nacelle pivot location 
for tilting proprotor configuration) 


rotor 2 hub location 


wing- body center of action 


horizontal tail center of action 


vertical tail center of action 


point off rotor disk (for induced velocity calculation 


rotor 1 

control inputs (deg) for all sticks centered ( r' f , = 0. 

< l)0« 

rotor 1 

(2) 0 lr 

rotor 1 

(3) e lH 

rotor 2 

(4) 0« 

rotor 2 

(5) 0„. 

rotor 2 

(6)0 U 

aircraft 

(7) t> f 

aircraft 

(8) ft. 

aircraft 

(9) K 

aircraft 

(10) K 

aircraft 

(11)0. 
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KOCFE 

KC'CFE 

KSCFE 

KPCFE 

PCC’FE 

PSCFE 

PPC’FE 


KFOCFE 

KROCFE 

IvFCCFE 

KRC'CFE 

KFSCFE 

KRSCFE 

KFPCFE 

KRPCFE 

PFC'CFE 

PRCC’FE 

PFPCFE 

PRPCFE 


KFCFE 

KTCFE 

KACFE 

KECFE 

IvRCFE 


description of control system A parameters 

are gains (deg per stick deflection). At/ 1 parameters 
are swashplate azimuth lead angles (deg) 

one rotor, single main rotor and tail rotor 
tilting proprotor configurations: 

Ao, collective stick to collective pitch. 

A'r, lateral cyclic stick to cyclic or differential 
collective pitch. 

A' s . longitudinal cyclic stick to cyclic pitch. 

A ;> , pedal to tail rotor collective or differential 
cyclic pitch. 


lateral cyclic stick to cyclic pitch 
(one rotor, or main rotor and tail rotor configurations) 
longitudinal cyclic stick to cyclic pitch 

pedal to differential cyclic pitch 
(tilting proprotor configuration only) 


tandem main rotor configuration 

A> n . collective stick to front collective pitch 

I\ru , collective stick to rear collective pitch 

AVr, lateral cyclic stick to front cyclic pitch 

I\rc • lateral cyclic stick to rear cyclic pitch 

A longitudinal cyclic stick to front collective pitch 

A'ffs, longitudinal cyclic stick to rear collective pitch 

A p/>, pedal to front cyclic pitch 

Ap/>. pedal to rear cyclic pitch 


A i'fc, lateral cyclic stick to front cyclic pitch 
Ai :, r« lateral cyclic stick to rear cyclic pitch 
Ar pp, pedal to front cyclic pitch 
A t /?/>, pedal to rear cyclic pitch 


aircraft controls (all configurations) 

J >j- collective stick to flaperon 
A/, collective stick to throttle 
A r/ . lateral cyclic stick to ailerons 
AV, longitudinal cyclic stick to elevator 
A,., pedal to rudder 



NEM 

number of airframe modes for which data supplied: 
(max 10). 

QMASS(NEM) 

generalized mass, , including rotors (slug or kg) 

QFREQ(NEM) 

generalized frequency, (Hz) 

QDAMP(NEM) 

structural damping, g s 

QDAMPA(NEM) 

aerodynamic damping, F (fk(Jk - 0(Q^ /(^pV 2 ))/0(q Nk /\ 
(ft 2 or m 2 ) 

QCNTRL(4.NEM) 

control derivatives, F </fr * = d(Q k / (^pY 2 ))/0b 
for bj.br^basby (ft 2 /rad or in 2 /rad) 

DOFSVM(NEM) 

integer vector designating type of mode 
(only required for flutter analysis with OPSYMM ^ 0) 
> 0 for symmetric 

< 0 for antisymmetric 

ZETAR1(3.NEM) 

ZETAR2(3.NEM) 

GAMAR1(3.NEM) 

GAMAR2(3,NEM) 

linear mode shape, at rotor 1 huh (ft /ft or m/m) 

linear mode shape. at rotor 2 hub (ft/ft or m/m) 

linear mode shape, ^ at rotor 1 hub (rad/ft or rad/m] 

linear mode shape, ^ at rotor 2 hub (rad /ft or rad/m] 

KPMCl(NEM) 

KPMSl(NEM) 

KPMC2(NEM) 

KPMS2(NEM) 

piteh/mast-bendine coupling (rad /ft or rad/m): 

K\u k = —dQ\<-/d q Sk for rotor 1 

K\ts k = —dtf\jdq, Sk for rotor 1 

K\n k — -~dti\ r /dq„ k for rotor 2 

K\is k = -d(J\»/dq. Sk for rotor 2 
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Aircraft aerodynamic characteristics: 


wing- body: 


I. HAW 

L n/q 

(ft 2 /rad or m 2 /rad) 
(ft 2 /rad or w 2 /rad) 

LFTFW 


LFTDW 

Lt'r/q 

(ft 2 /rad or m 2 /rad) 

AMAXW 

^ mar 

(deg) 

IWB 

*WB 

(deg) 

DRGOW 

f\YB = D{)/q 

(ft 2 or in 2 ) 

DRGYW 

fverl 

(ft 2 or m 2 ) 

DRGIW 

nelj =d(L/'j) 2 /d(Di/q) 

(ft 2 or m 2 ) 

DRGFW 

D <b } hi 

(ft 2 /rad or m 2 /rad) 

DRGDW 

D »,,hi 

(ft 2 /rad or m 2 /rad) 

MOMOW 

M,/q 

(ft* or in*) 

MOM AW 

MJq 

(ft* /rad or m'* /rad) 
(ft* /rad or m* /rad) 

MOMFW 

M,,/q 

MOMDW 

M^/q 

(ft* /rad or m* (rad) 

SIDEB 

y.hq 

(ft 2 /rad or m 2 f rad) 

SIDEP 

yy„hi 

(ft* /rad or m* /rad) 

SIDER 

yy/q 

(ft* /rad or m*/rad) 

ROLLB 

yjq 

(ft* /rad or in* /rad) 

ROLLP 

y-yjq 

(ft 4 /rad or m 4 /rad) 

ROLLR 

yyjq 

(ft 4 /rad or m 4 /rad) 

ROLLA 

yy,jq 

(ft* /rad or in’* /rad) 

YAWB 

y-.Jq 

(ft* /rad or m* /rad) 

YAWP 

yy,jq 

(ft 4 /rad or m 4 /rad) 

YAWR 

yy-.jq 

(ft 4 /rad or m 4 f rad) 
(ft* /rad or m*/ rad) 

YAWA 

y*,»jq 

Horizontal tail: 

LFTAH 

Ay 

(ft 2 /rad or m 2 /rad) 
(ft 2 / rad or m 2 /rad) 

LFTEH 

L\/q 

AMAXH 

O mar 

(<^g) 

IHT 

1 HT 

Vertical tail: 

(deg) 

LFTAY 

L «/q 

(ft 2 /rad or ta 2 /rad) 
(ft 2 /rad or m 2 /rad) 

LFTRY 

L *,hi 

AMAXY 

C* mar 

(deg) 

IYT 

hr 

(< le g) 
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Airframe interference: 


FETAIL 

LHTAIL 

HYTAIL 


f, = 0(L / <j ) / Of 

horizontal tail length, Iht- for e (ft or m) 
vertical tail height, hyp* for rr, positive up (ft or m) 

OPTINT 

= 0 

integer parameter controlling airframe/tail aerodynamic 
interference: 

to suppress (f = 0. and a = 0.) 



Engine and Drive Train Parameters: 

ENGPOS 

= 0 
= 1 
= •■> 
= 3 

integer parameter specifying drive train configuration: 
one rotor 

asymmetric, engine hv rotor 1 
asymmetric, engine by rotor 2 
symmetric 

IENG 


engine rotational inertia, /y 2 J f , for both engines 
if symmetric configuration (slug-// 2 or kg-m 2 ) 



drive train spring constants (ft-lb/rad or m-N/rad) 

KM ASTI 
KMAST2 
KICS 
KENG 


rotor 1 shaft, A a/, or A y 
rotor 2 shaft, A \/ 2 
interconnect shaft, r/ 2 2 A / or r/ 2 A / 
engine shaft, r p 2 I\p 

GSE 

GSI 


engine shaft structural damping, (} s . (ty degree of freedom) 
interconnect shaft structural damping, g s , 

(i'l degree of freedom) 

KEDAMP 


engine damping factor, k: 

typically — 1.0 for turboshaft engines, or 

10. for induction electric motors 

THRTLC 


dPp;/dtit (dimensional), for both engines 

if symmetric configuration: if the throttle variable, 0 t is only 

used for the governor, just the products: 

KrOPr/OH, = -0P/dr s 
K j OP KjOHi = -OP/Ovs 
must be correct (P = ilpQp = UfQe) 



governor proportional feedback gains (sec): 

KPGOYE 

IvPGOYl 

KPGOY2 


to throttle. I\p = 

to rotor 1 collective, Kp — — 0B/0v s 
to rotor 2 collective, I\p — ~dQ/di\ 

KIGOVE 

KIGOYl 

KIGOY2 


governor integral feedback gains: 
to throttle, I\f — —OOf /di\ s 
to rotor 1 collective, Kj — — d0/0t\ s 
to rotor 2 collective, AT = — 0Q /Oil's 
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T1G0VE 

T1GOY1 

T1GOV2 


T2GOVE 
T2GOV1 
T 2 GOV 2 


governor time lag T\ ~ 2CluJxA&ec\ 

throttle 

rotor 1 

rotor 2 

governor time lag r> = \ / (sec 2 

throttle 

rotor 1 

rotor 2 
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4.19 NLLOAD 


Namelist NLLOAD 




Airframe vibration: 

MYIB 

< 0 

number of stations for airframe vibration calculation 
and print, (max 10). 
to suppress 



airframe location for vibration calculation (ft or m): 

FSYIB(MYIB) 
BLMB(MYIB) 
WL\*IB(MVIB) 
ZETAY( 3.NEM.MMB ) 


fuselage station 
butt line 
water line 

linear mode shape, . at airframe vibration 
stations (ft /ft or m/m) 

MALOAD 

= 0 
< 0 

integer parameter controlling print of motion and 
aerodynamics: 
to suppress 
for only plots 

MHLOAD 

= 0 

integer parameter controlling print of hub and 
control loads: 
to suppress 

MR LOAD 

< o 

number of radial stations for blade section load 
calculation and print (max 20) 
to suppress 

RLOAD(MRLOAD) 


blade radial stations, r/R, for section loads 

MHARML 

< 0 

number of harmonics in loads analysis (max 30): 
for no harmonic analysis: suggest about MPSI/3 

NPOLAR 


integer parameter n for polar plots: 
symbol printed every n-tli step. 

NWKGMP(4) 

= 0 

integer parameter controlling wake geometry printer plot 
to suppress: 



( 1 ) top view 

(2) side view 

(3) back view 

(4) axial convection 

MWKGMP 

< 0 

number of azimuth stations at which wake geometry 
plotted (max 8): 
for no plots 

JWKGMP(MWKGMP) 


azimuth stations at which geometrv plotted 
(c = .JWKGMP * Af) 
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NPLOT(75 


) integer parameter controlling printer- plots of motion 

and aerodynamics: 

values: 

= 0 for no plot 

= 1 for time history 

= 2 for polar plot 

= 3 for noth 

(only time histories available for 1-4 and 68-75) 
elements: 

(1) bending motion 

(2) torsion motion 

(3) maximum circulation 

(4) A off rotor 

(5) o 

(6) M 

(7) A 

(8) O 

(9) G, 

(10) Cm 

(11) 

(12) r 

(13) U„ 

(14) Ur 

(15) 

(16) V 

(17) 0 

(18) 0 

(19) lag 

(20) flap 

(21) <><■//■ lift 

(22) «>//, drag 

(23) Orff, moment 

(24) lift 

(25) Mrff. drag 

(26) M, ff. moment 

(27) A, 

(28) A„ 

(29) A-. 

(30) interference Aj. 

(31) interference X u 

(32) interference A - 

(33) u n 

(34) v« 

(35) it -a 

(36) L/c 

(37) D/c 

(38) M/c 

(39) D,./c 
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(40) 

F,/<- 

(41) 

Fr/C 

(42) 

F. /c = C T a 

(43) 

Male 

(44) 

Fr/e 

(45) 

not used 

(40) 

not used 

(47) 

not used 

(48) 

C'r/(T 

(49) 

C'p, / <7 

(50) 


(51) 

Crjcr 

(52)* 

L 

(53)* 

D 

(54)* 

M 

(55)* 

D r 

(50)* 

F, 

(57)* 

F r 

(58)* 

F- = T 

(59)* 

Ma 

(00)* 

Fr 

(01) 

not use<l 

(02) 

not used 

(03) 

not used 

(04)* 

P 

(05)* 

Pi 

(00)* 

Pint 

(07)* 


(08) 

rotating frame root loads 

(09) 

nonrotating frame hub loads 

(70)* 

rotating frame root loads 

(71)* 

nonrotating frame hub loads 

(72) 

section loads, shaft axes 

(73) 

section loads, principal axes 

(74)* 

section loads, shaft axes 

(75)* 

section loads, principal axes 

^dimensional quantities 


for polar plots, last digit of integer part of (value /increment) is printed, if it is a multiple of 
NPOLAR: the increment is defined as follows: 


.01 plots 27-35 

.10 plots 6. 8-1G, 24-2G. 3G-51 

1.0 plots 5, 7, 17-23, 52-G1 

10. plots 62-6/ 
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KFATIG 


parameter. Iv, in fatigue damage calculation: suggest 3 or 4 

SENDUR(18) 


endurance limit, (dimensional force or moment) 

CMAT(18) 


material constant, C 

EXMAT(18) 


material exponent, M 
rot at ins frame root loads: 


(i) 

inplane shear, f r 


(2) 

axial shear, /,. 


(3) 

vertical shear, f 


(4) 

flap moment, in x 


( 5 ) 

lag moment, m - 


(6) 

control moment, m r 
nonrotatine frame lmh loads: 


(") 

drag force, H 


(8) 

side force, 1’ 


( 9 ) 

thrust, T 


(10) 

roll moment. M r 


(ID 

pitch moment. M y 


(12) 

torque, Q 

section loads (principal axes): 


(13) 

chord shear, 


(14) 

axial shear, /,. 


(15) 

normal shear, / 


(16) 

flapwise moment. m r 


(17) 

edgewise moment, m- 


(18) 

torsion moment, \n\ 

the S-N curve is 

approximated by .V = C/(S/Sfr - l) 2 

use Se < 0. or C < 0. 

to suppress damage fraction calculation: 

use M < 0. to suppress equivalent peak-to-peak load calculation 
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Far field rotational noise: 

MNOISE 

number of microphones (max 10): 
< 0 for no noise analysis 

RANGE(MNOISE) 

microphone range relative hub (ft or m) 

ELVATN(MNOISE) 

microphone elevation relative hub (deg), positive 
above rotor disk 

AZMUTH(MNOISE) 

microphone azimuth relative hul) (deg), defined 
as for rotor azimuth 

MHARMN(3) 

number of harmonics: 

(1) in noise calculation (max 500) 

(2) in aerodynamic load harmonic analysis (suggest MPSI/3) 

(3) in print of noise (< 0 for no print) 

MTIMEN(3) 

number of time steps (< 0 to suppress): 

(1) in period of noise calculation (max 500) 

(2) increment in noise print 

(3) increment in noise plot 

AXS(MRA) 

blade cross section area, A xs /c 2 at aerodynamic 
segments, for thickness calculation (typically 
.685 times thickness ratio) 

OPNOIS(4) 

integer parameter controlling noise calculation: 
= 0 to suppress 

= 1 for impulsive chord wise loading 
= 2 for distributed chordwise loading 

(1) lift noise 

(2) drag noise 

(3) radial force noise 

(4) thickness noise 
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